流程图
# 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;