lecture8 - psi=chi-psi enddo do j=0,jmax-1 x=dx*j...

Info iconThis preview shows pages 1–5. Sign up to view the full content.

View Full Document Right Arrow Icon
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 Set the problem up, initialize wave function, Q matrix
Background image of page 1

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon
do nt=1,ntmax call tridag(a,b,c,psi,chi,jmax)
Background image of page 2
Background image of page 3

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon
Background image of page 4
Background image of page 5
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: 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 Then integrate in time, output wave function at end QuickTime and a TIFF (LZW) decompressor are needed to see this picture. QuickTime and a TIFF (LZW) decompressor are needed to see this picture....
View Full Document

This note was uploaded on 08/08/2011 for the course PHZ 5156 taught by Professor Johnson,m during the Fall '08 term at University of Central Florida.

Page1 / 5

lecture8 - psi=chi-psi enddo do j=0,jmax-1 x=dx*j...

This preview shows document pages 1 - 5. Sign up to view the full document.

View Full Document Right Arrow Icon
Ask a homework question - tutors are online