Solving second order ODE eigenvalue problem using RK4

Hi, I'm trying to solve the following ODE using RK4:

y''+(lambda-x^2) y = 0 on [-L,L], subject to y(-L)=y(L)=0 (taking L=5 is fine).

I know the solution is lambda=2n+1, n some positive integer (some residual is fine), with y=exp(-x^2/2) (which obviously is very small at L=\pm 5)

I'm terribley bad at coding and I don't think my code will be any use, so I havent included it.

Please help. Thanks in advance.

[Jan]

You forgot to ask a question. How can we help you?

Of course your code would be useful, because if you show it, we can help you to fix the problems in it.

Answer 1

[RAKESH KUMAR]

i think you should look at this web-link: https://www.cs.ox.ac.uk/pseudospectra/eigtool/

Eigtool