50 likes | 187 Views
Set the problem up, initialize wave function, Q matrix. do j=0,jmax-1 x=dx*j psi(j)=(ar/(dsqrt(sigma)*sspi))* * dexp(-(x-x0)**2/(2.0d0*sigma**2))* * (ar*dcos(k0*x)+ai*dsin(k0*x)) write(6,100) x,real(psi(j)),imag(psi(j)) enddo
E N D
Set the problem up, initialize wave function, Q matrix do j=0,jmax-1 x=dx*j psi(j)=(ar/(dsqrt(sigma)*sspi))* * dexp(-(x-x0)**2/(2.0d0*sigma**2))* * (ar*dcos(k0*x)+ai*dsin(k0*x)) write(6,100) x,real(psi(j)),imag(psi(j)) enddo do j=0,jmax-1 a(j)=-ai*alpha/8.0d0 b(j)=ar*0.5d0+ai*alpha/4.0d0 c(j)=-ai*alpha/8.0d0 enddo a(0)=0.0d0 ; c(0)=0.0d0
Then integrate in time, output wave function at end do nt=1,ntmax call tridag(a,b,c,psi,chi,jmax) psi=chi-psi enddo do j=0,jmax-1 x=dx*j write(6,100) x,real(psi(j)),imag(psi(j)) enddo 100 format(f8.4,2f12.6) stop end