In[]:=
Clear["Global`*"]
In[]:=
L=30;
In[]:=
k=2
Out[]=
2
In[]:=
ℒ=-Laplacian[u[x],{x}]+(x^2)+*u[x]2
2
Sin
2*k*π*x
L
Out[]=
1
2
2
x
2
Sin
2πx
15
′′
u
In[]:=
sol=NDEigensystem[{ℒ,u[-30]==0,u[30]==0},u[x],{x,-30,30},10]
Out[]=
{0.177782,0.179474,0.181851,0.183384,0.191888,0.194827,0.197089,0.441937,0.445448,0.455733},InterpolatingFunction[x],InterpolatingFunction[x],InterpolatingFunction[x],InterpolatingFunction[x],InterpolatingFunction[x],InterpolatingFunction[x],InterpolatingFunction[x],InterpolatingFunction[x],InterpolatingFunction[x],InterpolatingFunction[x]
In[]:=
sol[[1]]
Out[]=
{0.177782,0.179474,0.181851,0.183384,0.191888,0.194827,0.197089,0.441937,0.445448,0.455733}
In[]:=
sol[[2]]
Out[]=
InterpolatingFunction[x],InterpolatingFunction[x],InterpolatingFunction[x],InterpolatingFunction[x],InterpolatingFunction[x],InterpolatingFunction[x],InterpolatingFunction[x],InterpolatingFunction[x],InterpolatingFunction[x],InterpolatingFunction[x]
In[]:=
Plot[Evaluate[sol[[2]],{x,-30,30}],PlotRangeAll]
Out[]=