一个数值求解 2个函数2个自变量偏微分方程的例子

dlpg070

一个数值求解 2个函数2个自变量偏微分方程的例子

希望对倪鹏举的问题有启发

1. Clear["Global`*"]
   
2. 
   
3. Print["2个函数2个自变量 NDVolve求解的简单例子  "]
   
4. 
   
5. s=NDSolve[{D[y1[x1,x2],x1]==D[y1[x1,x2],x2,x2],y1[0,x2]==0,y1[x1,0]==Sin[x1],y1[x1,5]==0,
   
6. D[y2[x1,x2],x1]==D[y2[x1,x2],x2,x2],y2[0,x2]==0,y2[x1,0]==(Sin[x1])^2,y2[x1,5]==0},{y1,y2},(*{t,x}\[Element]\[CapitalOmega]*){x1,0,10},{x2,0,5},MaxStepSize->0.05]
   
7. Print["显示y1[x1,x2]:"]
   
8. Plot3D[Evaluate[y1[x1,x2]/.s],{x1,0,10},{x2,0,5}, PlotRange->All]
   
9. Print["显示y2[x1,x2]:"]
   
10. Plot3D[Evaluate[y2[x1,x2]/.s],{x1,0,10},{x2,0,5}, PlotRange->All]
    

复制代码