clear;

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

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

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