clear; // *** 連立方程式の定義 *** A = [-2 1 0; 1 -2 1; 0 1 -2]; b = [-1; 0; 0]; // 境界条件 u0 = 1; u4 = 0; // *** ラプラス方程式の解 *** u = A \ b; // *** グラフのプロット *** // グラフ用のx座標 x = linspace(0,1,5); // グラフのプロット plot(x,[u0; u; u4],'-ob'); // グラフの装飾 xlabel("x"); ylabel("u"); zoom_rect([0,0,1.2,1.2]);