Язык Wolfram Language™

Создание галереи собственных функций лапласиана в шаре

Определим трехмерный лапласовский оператор.

In[1]:=

\[ScriptCapitalL] = -Laplacian[u[x, y, z], {x, y, z}];

Укажем однородные граничные условия Дирихле.

In[2]:=

\[ScriptCapitalB] = DirichletCondition[u[x, y, z] == 0, True];

Найдем 16 наименьших собственных значений и собственных функций в шаре.

In[3]:=

\[CapitalOmega] = Ball[{0, 0, 0}, 2];

In[4]:=

{vals, funs} = DEigensystem[{\[ScriptCapitalL], \[ScriptCapitalB]}, u[x, y, z], {x, y, z} \[Element] \[CapitalOmega], 16];

Собственные значения заданы через функцию BesselJZero.

In[5]:=

vals[[1]] // TraditionalForm

Out[5]//TraditionalForm=

Создадим галерею собственных функций.

In[6]:=

Grid[Partition[ ParallelTable[ DensityPlot3D[ funs[[i]] // N // Evaluate, {x, y, z} \[Element] \[CapitalOmega], Boxed -> False, Axes -> False, ColorFunction -> Hue, Method -> {"ShrinkWrap" -> True}, ImageSize -> 125], {i, 16}], 4]]

Out[6]=