【正文】
u(41,j)=(1/8)*[sin(pi*j*h)] u(40,j)=(1/8)*[sin(pi*j*h)] end %************近似解 ************** for i=39:1:2 for j=2:40 u(i,j)=u(i+1,j1)+u(i+1,j+1)u(i+2,j) end end %**************精確解 *************** for i=1:40 for j=1:40 x1=i*h。 t1=j*h。 u1(i,j)=1/16*sin(pi*(x1t1))+1/16*sin(pi*(x1+t1)) end end