function u = heat1d(n)
% 熱伝導シミュレーション（１次元、定常）
%  n: 節点の数

d = 5; % 棒の長さ
f = 20; % 発熱量
h = d/(n - 1); % 分割の幅
A = zeros(n); b = zeros(n,1); % 初期化
A(1,1) = 1; b(1) = 0; % u(1) = 0
A(n,n) = 1; b(n) = 0; % u(n) = 0
for i=2:n-1
    A(i,i) = 2;
    A(i,i-1) = -1;
    A(i,i+1) = -1;
    b(i) = f*h^2;
end
u = A\b; % Ax=bを解く（Macの場合、\はoption+\）