에어리 방정식의 스튀름-리우빌 문제 풀이
에어리 방정식의 고유값 문제를 해결합니다.
In[1]:=
![Click for copyable input](assets.ko/solve-a-sturm-liouville-problem-for-the-airy-equat/In_35.png)
eqns = {-y''[x] + x y[x] == \[Lambda] y[x], y[0] == 0, y[1] == 0};
In[2]:=
![Click for copyable input](assets.ko/solve-a-sturm-liouville-problem-for-the-airy-equat/In_36.png)
sol = DSolveValue[eqns, y[x], x, Assumptions -> 1 < \[Lambda] < 300];
문제의 고유값 목록을 만듭니다.
In[3]:=
![Click for copyable input](assets.ko/solve-a-sturm-liouville-problem-for-the-airy-equat/In_37.png)
eigvals = {ToRules[sol[[1, 1, 2]]]}
Out[3]=
![](assets.ko/solve-a-sturm-liouville-problem-for-the-airy-equat/O_21.png)
이 문제의 고유 함수를 구합니다.
In[4]:=
![Click for copyable input](assets.ko/solve-a-sturm-liouville-problem-for-the-airy-equat/In_38.png)
eigfuns = (sol[[1, 1, 1]] /. {C[1] -> 1})
Out[4]=
![](assets.ko/solve-a-sturm-liouville-problem-for-the-airy-equat/O_22.png)
1<λ<300의 범위에서 고유 함수를 플롯합니다.
In[5]:=
![Click for copyable input](assets.ko/solve-a-sturm-liouville-problem-for-the-airy-equat/In_39.png)
Plot[Evaluate[eigfuns /. eigvals], {x, 0, 1}]
Out[5]=
![](assets.ko/solve-a-sturm-liouville-problem-for-the-airy-equat/O_23.png)