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

流程图

# 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*(k1+2*k2+2*k3+k4)/6;}

if(psi<=0.001&&psi>=-0.001&&psi1>=-0.00001&&psi1<=0.00001)

{psi=0.0001;psi1=0;

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*(k1+2*k2+2*k3+k4)/6;

fprintf(fp,"%f %f

%f\n",x,psi,e);}

}

fclose(fp);

}

运行结果

输入：

please input the number of xmax

20.

please input the number of Emax

4000.

please input the number of K

please input the number of date group

N1,N2

1000 4000

图像：

能量本征值1的图像：

能量本征值2的图像：

00087

3040.000000

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

傅里叶谱方法-傅里叶谱方法求解一维 KdV 方程及其Matlab程序实现
3.3 傅里叶谱方法求解复杂偏微分方程 (组) 3.3.1 一维 KdV 方程背景介绍科特韦赫-德弗里斯方程(英语:Korteweg-De Vries equation),一般简称KdV方程,是1 ...
问题: 一维Schrodinger方程的双线性Strichartz估计的最佳性
一维的Schrodinger方程的双线性Strichartz估计, 有一个看似简单但目前仍没答案的问题, 表述如下: 假设$f,g \in L^2(\mathbf{R})$, 且$\widehat{f ...
深度学习求解一维burgers方程和Galerkin求解泊松方程
一维burgers方程介绍 { u t + u u x − ( 0.01 / π ) u x x
隐方程求解一维抛物型方程（热传导方程）
记录一下数值计算课大作业设时间t和空间x 离散:xi_{i}i=ih , tn_{n}n=nΔ\DeltaΔ*t (i = 0,1,-I , n = 0,1,-) 设 U(xi_{i}i,tn ...
numerov算法matlab,Numerov算法求解一维薛定谔方程研究.doc
Numerov算法求解一维薛定谔方程研究摘要本文主要讲述的是运用MATLAB对一维定态薛定谔方程求解中遇到问题的分析.问题指的是在运用文献[1]所提供的程序,解一维定态薛定谔方程时出现了波函数与 ...
matlab向量的模_基于MATLAB使用矩阵方法求解一维定态薛定谔方程
摘要:此文介绍了一种使用MATLAB求解一维定态薛定谔方程的方法.利用充分格式进行离散化,得出相应的矩阵方程,用MATLAB求解本征值和本征函数.此方法简单可靠,可以处理各种时间无关的束缚态问题.所用 ...
求解欧拉方程的c语言,用有限体积方法求解欧拉方程
<用有限体积方法求解欧拉方程>由会员分享,可在线阅读,更多相关<用有限体积方法求解欧拉方程(12页珍藏版)>请在人人文库网上搜索. 1.有限体积法求解二维可压缩Euler方程计 ...
c语言实验报告5数组,c语言实验报告五一维数组.doc
c语言实验报告五一维数组 <C程序设计>课程实验报告学院: 班级: 姓名: 学号: 实验设备:计算机1台2011年3月1日实验项目名称一维数组掌握一维数组的定义.初始化,以及一维数组 ...
c语言实现二分法_C语言实现二分法求解方程在区间内的根
C语言实现二分法求解方程在区间内的根. 设有非线性方程: 其中, 为上连续函数且设 (不妨设方程在内仅有一个实根),求上述方程实根的二分法过程,就是将含根区间[a,b]逐步分半,检查函数值符号的变 ...

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

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

最新文章

热门文章