调用二分法子程序解方程出错，求助

Karmi

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

查看代码

打印?

01	Program test

02	Implicit none

03	Integer :: i

04	Real :: Maact1,Maact2

05	Real,dimension(3) ::fact

06	fact(1)=26.

07	fact(2)=30.

08	fact(3)=32.

09	Do i=1,3                                      !输出的结果Maact1和Maact2结果一致，都为1.985840

10	Call equation(fact(i),Maact1)               !但是把9-12行注释掉的话，

11	Write(*,*) Maact1                           !输出的Maact2结果为2.1938，是正确的

12

End Do

13	Call equation(31.57,Maact2)

14	Write(*,*) Maact2

15	End Program test

16

17	Subroutine equation(f,Ma)              !已知函数值，根据普朗特-迈耶函数求马赫数,1<=Ma<=3

18	Implicit none

19	Real,intent(in) :: f                   !输入普朗特-迈耶函数值

20	Real,intent(out) :: Ma                 !输出马赫数

21	Real :: eps=0.001                      !计算精度

22	Real :: Ma1=1.0                        !区间下限

23	Real :: Ma2=3.0                        !区间上限

24	Real :: f_Ma1,f_Ma2,f_Ma               !对应函数值

25	Real :: gamma=1.4

26	Real,External :: P_Ma                  !普朗特迈耶函数

27

!二分法解方程

28	Ma=(Ma1+Ma2)/2

29	f_Ma=P_Ma(Ma)-f

30	f_Ma1=P_Ma(Ma1)-f

31	f_Ma2=P_Ma(Ma2)-f

32	Do while((abs(f_ma)>eps).and.((Ma2-Ma1)>eps))

33	If(f_Ma*f_Ma1<0.) then

34

Ma2=Ma

35	Ma=(Ma1+Ma2)/2

36	f_Ma2=f_Ma

37

Else

38

Ma1=Ma

39	Ma=(Ma1+Ma2)/2

40	f_Ma1=f_Ma

41

End If

42	f_Ma=P_Ma(Ma)-f

43

End Do

44

45	!f_Ma1=P_Ma(Ma1)-f

46	!f_Ma2=P_Ma(Ma2)-f

47

!Do

48	!  Ma=(Ma1+Ma2)/2

49	!  f_Ma=P_Ma(Ma)-f

50	!  If ((abs(f_Ma)<=eps).or.abs(Ma1-Ma2)<=eps) exit

51	!  If (f_Ma1*f_Ma<0) then

52	!    Ma2=Ma

53	!    f_Ma2=f_Ma

54

!  Else

55	!    Ma1=Ma

56	!    f_Ma1=f_Ma

57

!  End if

58

!End do

59	End subroutine equation

60

61	Real Function P_Ma(Ma)

62	Implicit None

63	Real,Intent(in) :: Ma

64	Real :: gamma=1.4

65	P_Ma=sqrt((gamma+1)/(gamma-1))*atand(sqrt((gamma-1)/(gamma+1)*(Ma**2-1)))-atand(sqrt(Ma**2-1))

66	End FUnction P_Ma

当我调用子程序的时候，输入如果是一个具体的数字的话，就能得到正确结果，但如果输入一个赋值的变量名的话就结果错误。
二分法也换了种写法，但还是同样的问题，是我的调用有什么问题吗？

li913

在定义变量时进行赋值，变量具有save属性，会保存前一次的计算结果。所以，把22,23行改一下就是。

Karmi

li913 发表于 2020-5-19 12:10
在定义变量时进行赋值，变量具有save属性，会保存前一次的计算结果。所以，把22,23行改一下就是。 ...

原来是这样，非常感谢，那么我把Ma1，Ma2作为参数输入来处理也是可以的了。