h=L/N; e=ones(N-1,1); a=spdiags([-e 2*e -e],-1:1,N-1,N-1)/(h*h);
ずっと以前は、疎行列関係の命令を知らずに次のように作っていた。
h=L/N; e=ones(N-1,1); a=(2*eye(N-1,N-1)-diag(e,1)-diag(e,-1))/(h*h);