Wolfram Language

Ecuaciones diferenciales parciales

Solución de un problema de Dirichlet para una ecuación de Laplace

Especifique la ecuación de Laplace en 2D.

In[1]:=
Click for copyable input
leqn = Laplacian[u[x, y], {x, y}] == 0;

Establezca una condición de Dirichlet en un rectángulo.

In[2]:=
Click for copyable input
\[CapitalOmega] = Rectangle[{0, 0}, {1, 2}];
In[3]:=
Click for copyable input
dcond = DirichletCondition[ u[x, y] == Piecewise[{{UnitTriangle[2 x - 1], y == 0 || y == 2}}, 0], True];

Resuelva el problema de Dirichlet.

In[4]:=
Click for copyable input
sol = DSolveValue[{leqn, dcond}, u[x, y], {x, y} \[Element] \[CapitalOmega]] // FullSimplify
Out[4]=

Extraiga los primeros 300 términos de la suma de Inactive.

In[5]:=
Click for copyable input
asol = sol /. {\[Infinity] -> 300} // Activate;

Visualice la solución en el rectángulo.

In[6]:=
Click for copyable input
Plot3D[asol // Evaluate, {x, y} \[Element] \[CapitalOmega], PlotRange -> All, PlotTheme -> "Business"]
Out[6]=

Ejemplos relacionados

de en fr ja ko pt-br ru zh