微分方程数值解法实验报告样稿

_2_0,inf)<errРР flag=0;Р endР u_2_0=u_2;Р m_2=m_2+1;%迭代次数值РendРuu=zeros(n+1);Рfor mm=1:n-1Р for nn=1:n-1Р uu(nn+1,mm+1)=u_2((mm-1)*(n-1)+nn,1);Р endРendР%Gauss-Seidel迭代差分解图像Рx=[0:1/n:1];Рy=[0:1/n:1];Рfigure(3)Рmesh(x,y,uu); title('Gauss-Seidel迭代差分解图像')РРРSOR迭代法Рfunction[u_3,m_3]=SOR_Solve(A,b,err,w,n)РD=diag(A);РD=diag(D);РL=-tril(A,-1);РR2=-triu(A,1);РB1=(D-w*L)\((1-w)*D+w*R2);Рg=w*(D-w*L)\b;РРm_3=0;Рu_3_0=zeros(size(b,1),1);%初始迭代值Рflag=1;Рwhile flagР u_3=B1*u_3_0+g;Р if norm(u_3-u_3_0,inf)<errР flag=0;Р endР u_3_0=u_3;Р m_3=m_3+1;%迭代次数值РendРuu=zeros(n+1);Рfor mm=1:n-1Р for nn=1:n-1Р uu(nn+1,mm+1)=u_3((mm-1)*(n-1)+nn,1);Р endРendР%SOR迭代差分解图像Рx=[0:1/n:1];Рy=[0:1/n:1];Рfigure(4)Рmesh(x,y,uu); title('SOR迭代差分解图像') РРPR方法РРclc;