我用k表示个数，为什么出来的个数和k对应不起来？

小范范1989

[Fortran] 纯文本查看 复制代码

查看代码

打印?

01	program test

02	implicit none

03	integer :: k, j

04	real :: E, T1, F1, F, T, V, a

05	real, PARAMETER :: pi = 3.1415926535897932

06	write (*,*) '输入V和a，都是原子单位'

07	read (*,*) V, a

08	open (unit=100,file='E.dat')

09	open (unit=200,file='k.dat')

10	k=int(a*sqrt(2.0*V)/pi)+1   !解的总个数，通过横坐标ε进行判断

11	write (200,*) k             !看一下解的个数，和后面得到的E的个数应该是相等的

12	j=1                         !初始解个数为1个

13	E=0.001                     !开始寻找时E的初始数值

14	F1=-1.0                     !偶宇称时，表达式在E=0时的初始值

15	T1=-1.0                     !奇宇称时，表达式在E=0时的初始值

16	do while (j .LE. k)         !和解的总个数K比较，得到进行的循环的次数

17	F=sqrt(E/(V-E))*tan(0.5*a*sqrt(2.0*E))-1.0  ! 偶宇称时化简后的表达式

18	T=sqrt(E/(V-E))*tan(0.5*a*sqrt(2.0*E)+pi/2.0)-1.0   !奇宇称时化简后的表达式

19	if (F*F1 .LE. 0.0) then                           !寻找满足的E,它使f的符号变化。

20	if(abs(F) .LT. abs(F1) .and. F .GT. -1.0) write (100,*) V, E !看看谁和0更近，去进的时候的E

21	if(abs(F) .GT. abs(F1) .and. F .GT. -1.0) write (100,*) V, E-0.001  !因为0.001是E的步长

22	j=j+1                    !偶宇称得到一个解

23

end if

24	if (T*T1 .LE. 0.0)  then

25	if (abs(T) .LT. abs(T1) .and. T .GT. -1.0) write (100,*) V, E

26	if (abs(T) .GT. abs(T1) .and. T .GT. -1.0) write (100,*) V, E-0.001

27	j=j+1                    !奇宇称时的一个解

28

end if

29	F1=F       !把新的f放到f1中，进行下一步的比较

30	T1=T       !把新的t放到t1中，在进行下一步比较

31	E=E+0.001 ! E的步长

32	end do   !do 循环结束，这样随着E的增大，得到了k个解

33	end program test

teddy2003

是怎么对应不起来啊，程序看样没问题啊。

Saphy

虽然不了解你算的东西，但纯从代码方面来看，问题有两个
1. 当F=F1，T=T1时，不会输出E
2. 当F，T<-1时，也不会输出E
其中2是主要问题，你可以从这方面检查你的计算公式

小范范1989

Saphy 发表于 2014-4-14 12:09
虽然不了解你算的东西，但纯从代码方面来看，问题有两个
1. 当F=F1，T=T1时，不会输出E
2. 当F，T ...

谢谢，我会好好看看的，之前没太注意

小范范1989

teddy2003 发表于 2014-4-13 22:09
是怎么对应不起来啊，程序看样没问题啊。

按我的分子，k和个数应该是相等的，但是算出来后确实不相等的