jpattern.m

来自「利用龙格库塔法（等步长」· M 代码 · 共 10 行

function A=jpattern(t,u)
global N
A=zeros(N);
for i=1:N-1
    A(i,i)=1;
    A(i,i+1)=1;
    A(i+1,i)=1;
end
A(N,N)=1;
spy(A)