打靶法c语言程序,打靶法求解一维Schrodinger方程程序示例

流程图

# include

# include

# defineh 10

floatksqu(float

x,float e)

{float

k=0.1,m=0.01,ksq,v;

v=k*(x-1)*(x-1)*(x-2)*(x-2);

ksq=2*m*(e-v)/(h*h);

returnksq;

}

intmain()

{FILE *fp;

floatx,e,de=1;

floatdx,xmin,xmax,k1,k2,k3,k4,l1,l2,l3,l4,k,psi,psi1;

fp=fopen("1.txt","w+");

xmin=0;xmax=20;dx=0.05;

for(e=0;e<=6000;e=e+1)

{psi=0.0001;psi1=0;dx=0.1;

for(x=xmin;x<=xmax;x=x+dx)

{k=-ksqu(x,e);

k1=k*psi;

l1=psi1;

k2=k*(psi+dx*l1/2);

l2=psi1+dx*k1/2;

k3=k*(psi+dx*l2/2);

l3=psi1+dx*k2/2;

k4=k*(psi+l3*dx);

l4=psi1+dx*k3;

psi=psi+dx*(l1+2*l2+2*l3+l4)/6;

psi1=psi1+dx*(