Matlab Barrera de potencial con matlab

[DBG]

¡Buenas! estoy intentando resolver el sistema cuántico de la barrera de potencial empleando Matlab, en concreto quiero emplear el método de Runge-Kutta. El problema es que no sé como hallar las condiciones iniciales del sistema (o sea, ψ(0) y dψ(0)).
¿Alguien puede ayudarme? no se puede obtener la ecuación de onda de forma analítica.
Gracias

[jorgext]

Re: Barrera de potencial con matlab

Para resolver la ecuación de Schrödinger es conveniente que revises el método de Numerov.

[DBG]

Re: Barrera de potencial con matlab

Al final he empleado Runge-Kutta, aunque no es ni de lejos tan preciso, ya que con el método de Numerov no me manejo muy bien; lo que he hecho al final es tomar como punto inicial de integración el primer punto de la tercera región, donde la función de onda siempre puede obtenerse analíticamente, y es ψ(x)= C exp(-ikx), se escoge C=1 y después se normaliza la función de onda numéricamente. A partir de aquí, se integra la ecuación de onda desde el punto a al punto inicial (es decir, al revés en el espacio), y desde el punto a hasta el final.