精选优质文档-倾情为你奉上用显式格式求解二维抛物型偏微分方程2010-05-14 10:41function varargout=liu(varargin)T=1;a=1;h=1/32;dt=1/200;X,T,Z=chfenmethed(h,dt,a,T);mesh(X,T,Z(:,:,3);shading flat;% xlabel(X,FontSize,14);% ylabel(t,FontSize,14);% zlabel(error,FontSize,14);% title(误差图);function X,Y,Z=chfenmethed(h,dt,a,T);%求解下问题%u_t-a*(u_xx+u_yy)=f(x,y,t) 0x,y1,0tT%u(0,y,t)=g0,u(1,y,t)=g1, %u(x,y,0)=d %h离散x y方向的步长%dt离散t方向的步长x=0:h:1;y=x;t=0:dt:T;m=length(x);n=length(t);r=a*dt/h2
