slove_31.m

抛物型方程古典&c-n格式差分算法 dasdsadsadasdsdsadas

function x=slove_31(a,d)
m=length(d);
w=zeros(1,m);
for i=1:m
    b(i)=a(i,i);
end
aL=zeros(1,m);
aL(1)=a(2,1);
for i=1:m-1
    w(i)=a(i,i+1);
    aL(i+1)=a(i+1,i);
end
d(1)=d(1)/b(1);
w(1)=w(1)/b(1);
for i=2:m-1
    b(i)=b(i)-aL(i)*w(i-1);
    w(i)=w(i)/b(i);
    d(i)=(d(i)-aL(i)*d(i-1))/b(i);
end
d(m)=(d(m)-aL(m)*d(m-1))/(b(m)-aL(m)*w(m-1));
for i=m-1:-1:1
    d(i)=d(i)-w(i)*d(i+1);
end
x=d