Nugroho's blog.: Even Wave Equation on Infinite-depth Potential Well (with Some Variation Potential in Between) using Shooting Method,

## Wednesday, March 30, 2016

### Even Wave Equation on Infinite-depth Potential Well (with Some Variation Potential in Between) using Shooting Method,

It's from previous code, I use it to simulate even function.

`from pylab import *figure(3)n   = 39psi0= zeros(n) #psi = psi0     # x   = linspace(0,1,n)V   = zeros(n)#V   = 39*pow(x,2)for i in arange(n):    if i<n/2.:        V[i]    = 73.    else:        V[i]    = 0.       psi0[0] = 1.psi0[1] = psi0[0] #for odd function, use another value, `
`                  #but psi0[0] must be zerot   = 0dx  = 1./nE   = 1.dE  = .1err = .05while t< 771:    #k   = 2*dx*dx*(E-V)    for i in range (1,n-1):        k   = 2*dx*dx*(E-V[i])        psi[i+1]  = 2*psi0[i]-psi0[i-1]-k*psi0[i]    psi0    = psi    if abs(psi[n-1])<err:        print E        plot(x,psi)    t   += 1    E   += dExlabel('x')ylabel('psi')title('Fungsi Gelombang')grid(True)savefig("shooting.png")figure(2)plot(x,V)xlabel('x')ylabel('V')title('Potensial')grid(True)show()`
.
Here's some result with different potential V