li40.m

这是机械工业张德丰《MATLAB与外部程序接口编程》随书源码

n=2000;
D=sparse(1:n,1:n,-2*ones(1,n),n,n);
E=sparse(2:n,1:n-1,ones(1,n-1),n,n);
S=E+D+E'  %生成稀疏矩阵S
b=ones(n,1)
tic;  %计时开始
x=S\b;
t1=toc  %计时结束