lecture5 - Specification part Execution part do m=1,nmax

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

View Full Document Right Arrow Icon
Project 1 results nmax=100 E N (0) 2nd order perturbation theory Exact N=1 0.0493 -0.3793 -0.4012 N=2 0.1974 0.1869 0.1866 N=3 0.4441 0.2737 0.3048 N=4 0.7896 0.7531 0.7525 N=5 1.2337 1.0901 1.0902
Background image of page 1

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

View Full DocumentRight Arrow Icon
Dependence on nmax Nmax= 10 Nmax= 20 Nmax= 50 Nmax= 100 Nmax= 200 N=1 -0.3858 -0.4009 -0.4012 -0.4012 -0.4012 N=2 0.1869 0.1866 0.1866 0.1866 0.1866 N=3 0.3069 0.3048 0.3048 0.3048 0.3048 N=4 0.7538 0.7525 0.7525 0.7525 0.7525 N=5 1.0944 1.0904 10902 1.0902 1.0902
Background image of page 2
Wave function n=1
Background image of page 3

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

View Full DocumentRight Arrow Icon
Wave function n=2
Background image of page 4
Wave function n=3
Background image of page 5

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

View Full DocumentRight Arrow Icon
Wave function n=4
Background image of page 6
Wave function n=5
Background image of page 7

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

View Full DocumentRight Arrow Icon
implicit none INTEGER, PARAMETER :: nmax=100 INTEGER, PARAMETER :: nbin=1000 REAL*8, PARAMETER :: L=10.0d0 REAL*8 :: pi,fac REAL*8 , dimension(nmax) :: w,fv1,fv2,wpert,e0 REAL*8 , dimension(nmax,nmax) :: hamil,z REAL*8 :: psi,x,dx INTEGER :: matz,ierr,m,n,nb pi=4.0d0*datan(1.0d0) hamil=0.0d0 wpert=0.0d0 dx=L/nbin
Background image of page 8
Background image of page 9

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

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

Unformatted text preview: Specification part Execution part do m=1,nmax e0(m)=(pi*m)**2/(2.0d0*L**2) enddo do m=1,nmax hamil(m,m)=e0(m) do n=1,nmax fac=-(3.0d0/(2.0d0*L))*dsqrt(2.0d0*pi/L)* * (dcos(pi*(m-n)/2.0d0)*dexp(-(n-m)**2*pi**2/(2.0d0*L**3))-* dcos(pi*(m+n)/2.0d0)*dexp(-(n+m)**2*pi**2/(2.0d0*L**3))) hamil(m,n)=hamil(m,n)+fac if(m.eq.n) wpert(m)=wpert(m)+hamil(m,m) if(m.ne.n) wpert(m)=wpert(m)+fac**2/(e0(m)-e0(n)) enddo write(6,*) m,hamil(m,m) enddo matz=1 call rs(nmax,nmax,hamil,w,matz,z,fv1,fv2,ierr) do n=1,nmax write(6,*) w(n),wpert(n),e0(n) enddo fac=dsqrt(2.0d0/L) do m=1,5 do nb=1,nbin psi=0.0d0 do n=1,nmax x=dx*nb psi=psi+fac*z(n,m)*dsin(pi*x*n/L) enddo write(6,*) x,psi enddo write(6,*) enddo stop...
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 / 10

lecture5 - Specification part Execution part do m=1,nmax

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

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