Resolva um problema básico de Sturm–Liouville
Resolva um problema de valor próprio com condições de Dirichlet.
In[1]:=

sol = DSolve[{y''[x] + \[Lambda] y[x] == 0, y[0] == 0, y[\[Pi]] == 0},
y[x], x]
Out[1]=

Crie uma tabela com as primeiras 5 funções próprias.
In[2]:=

eigfuns =
Table[y[x] /.
sol[[1]] //. {\[FormalN] -> i, \[Lambda] -> \[FormalN]^2} /. {C[
1] -> 1}, {i, 5}]
Out[2]=

Faça um gráfico das funções próprias.
In[3]:=

Plot[Evaluate[eigfuns], {x, 0, Pi}]
Out[3]=

Resolva um problema de valores próprios com condições de Neumann.
In[4]:=

sol = DSolve[{y''[x] + \[Lambda] y[x] == 0, y'[0] == 0,
y'[\[Pi]] == 0}, y[x], x]
Out[4]=

Crie uma tabela com as primeiras 5 funções próprias.
In[5]:=

eigfuns =
Table[y[x] /.
sol[[1]] //. {\[FormalN] -> i, \[Lambda] -> \[FormalN]^2} /. {C[
1] -> 1}, {i, 5}]
Out[5]=

Faça um gráfico das funções próprias.
In[6]:=

Plot[Evaluate[eigfuns], {x, 0, Pi}]
Out[6]=
