追赶法求解三对角线性方程组

作业五：编写用追赶法求解三对角线性方程组的标准程序，并求下列方程组的解。matlab1233212417259xxxfunction[x,flag]=zgf(A,b)%A为方程组的系数矩阵%b为方程组的右端项%x为方程组的解%flag为指标向量，flag=‘failure’表示失败，flag=‘OK’表示成功[n,m]=size(A);n=length(b);[L,U,flag]=LU_decom(A);fori=1:nifi==1y(i)=b(i)/L(i,i);elsey(i)=(b(i)-L(i,i-1)*y(i-1))/L(i,i);endendx(n)=y(n)/U(n,n);fork=n-1:-1:1y(k)=y(k)-U(k,k+1)*x(k+1);x(k)=y(k)/U(k,k);end%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%function[L,U,flag]=LU_decom(A)[n,m]=size(A);ifn~=merror('TherowsandcolumnsofmatrixAmustbeequal!');return;endL=eye(n);U=zeros(n);flag='OK';fork=1:nforj=k:nz=0;forq=1:k-1z=z+L(k,q)*U(q,j);endU(k,j)=A(k,j)-z;endifabs(U(k,k))epsflag='failure';return;endfori=k+1:nz=0;forq=1:k-1z=z+L(i,q)*U(q,k);endL(i,k)=(A(i,k)-z)/U(k,k);endend流程图开始赋值LU=A分解由Ly=b求y由Ux=y求出x输出x运行结果