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]=
