参考文献:科学计算中的偏微分方程有限差分法
B=[4 -2 0 0 0;-1 4 -1 0 0;0 -1 4 -1 0;0 0 -1 4 -1;0 0 0 -2 4]; I=eye(5); Z=zeros(5); A=[B -2*I Z Z Z;-I B -I Z Z;Z -I B -I Z;Z Z -I B -I;Z Z Z -2*I B];
为了保持形式上的对称,试着将A改写为
A=[1.0/2*B -I Z Z Z;-I B -I Z Z;Z -I B -I Z;Z Z -I B -I;Z Z Z -1*I 1.0/2*B];