hothch.m

都是大家需要的程序,大家要下的快下,大家要多多支持

a=0;b=1;m=10;T=0.5;N=50;af=1;
f=inline('sin(pi*x)','x');
h=(b-a)/m;
k=T/N;
lmd=af^2*k/h^2;
x=linspace(a,b,m+1);
x=x(2:m);
i=1:m-1;
u=f(i.*h);
for j=1:N
    t=j*k;
    u=trisys(-lmd*ones(m-2,1),1+2*lmd*ones(m-1,1),-lmd*ones(m-2,1),u);
end
%u'
true=exp(-pi^2*T).*sin(pi*x);
error=abs(u-true);
re=[x'     u'       true'        error']