【正文】
初值 ************** for i=1:41 u(i,1)=0 u(i,41)=0 end for j=2:40 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