clear;

// xの刻み幅
dx = 1 / 4;
// グラフ用のx座標
x = [0:dx:1];

// *** 連立方程式の定義 ***
A = [1 -1  0  0  0;
     1 -2  1  0  0;
     0  1 -2  1  0;
     0  0  1 -2  1;
     0  0  0  0  1];
b = [dx;
     0;
     0;
     0;
     0];

// *** ラプラス方程式の解 ***
u = A \ b;

// *** グラフのプロット ***
// グラフのプロット
plot(x,u,'-ob');
// グラフの装飾
xlabel("x");
ylabel("u");
zoom_rect([0,0,1.2,1.2]);