clc clear all syms alpha n=2; %n=input('n=');%未知数个数 dt=0.01; x0=zeros(n,1); d=zeros(n); for i=1:n for j=1:n if(i==j) d(i,j)=1; end end end