Résolvez un problème de Dirichlet pour l'équation de Helmholtz
Spécifiez une équation de Helmholtz en 2D.
In[1]:=
heqn = {Laplacian[u[x, y], {x, y}] + 5 u[x, y] == 0};Définissez les conditions de Dirichlet pour l'équation dans un rectangle.
In[2]:=
bc = {u[x, 0] == UnitTriangle[x - 2]/2, u[x, 2] == 0, u[0, y] == 0,
u[4, y] == 0};Résolvez le problème de Dirichlet en utilisant DSolveValue.
In[3]:=
(sol = DSolveValue[{heqn, bc}, u[x, y], {x, y}]) // TraditionalFormOut[3]//TraditionalForm=
Extrayez les 30 premiers termes de la somme Inactive.
In[4]:=
fsol = sol /. \[Infinity] -> 30 // Activate;Visualisez la solution approximative.
In[5]:=
Plot3D[fsol // Evaluate, {x, 0, 4}, {y, 0, 2}, PlotRange -> All,
PlotTheme -> "Scientific"]Out[5]=
