****************************************************************Р%追赶法Рfunction [x]=zhuiganfa(a,b,c,d)Р%对角线下方的元素,个数比A少一个Р% %对角线元素Р%对角线上方的元素,个数比A少一个Р%d为方程常数项Р%用追赶法解三对角矩阵方程Рr=size(a);Рm=r(2);Рr=size(b);Рn=r(2);Рif size(a)~=size(c)|m~=n-1|size(b)~=size(d)Р error('变量不匹配,检查变量输入情况!');РendР%%Р%LU分解Рu(1)=b(1);Рfor i=2:nР l(i-1)=a(i-1)/u(i-1);Р u(i)=b(i)-l(i-1)*c(i-1);Р v(i-1)=(b(i)-u(i))/l(i-1); РendР%追赶法实现Р%%Р%求解Ly=d,"追"的过程Рy(1)=d(1);Рfor i=2:nР y(i)=d(i)-l(i-1)*y(i-1);РendР%%Р%求解Ux=y,"赶"的过程Рx(n)=y(n)/u(n);Рfor i=n-1:-1:1Р x(i)=y(i)/u(i);Р x(i)=(y(i)-c(i)*x(i+1))/u(i);РendР%********************************************************************Р(5.4)精确解函数Р%t时刻,u的取值;Рfunction [ f]=uexact(x,y,t)Рf=sin(x*pi)*cos(y*pi)*exp(-pi*pi/8*t);Р%********************************************************************