求助，关于imsl5.0函数使用

风和

想用imsl解非线性方程组的方法求计算程序的解（就是逐步逼近），用的imsl里的是zreal函数，但是这个函数非要用external，外部函数，我这整套程序都是在算这个结果，总不至于全部都放外边，求大神支招。

风和

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

查看代码

打印?

01	program main

02	include 'link_f90_static.h'

03

use test

04

05	use numerical_libraries

06

07

call ttt

08

09

10

pause

11	end program

12

13

14

15	module test

16	include 'link_f90_static.h'

17	use numerical_libraries

18

contains

19

20	real function f(x)  ！样例，实际主计算程序比这个多的多

21	implicit none

22

real x

23	f=X**2+5*X+3

24	!  write(*,*)x,f

25

return

26	end function

27

28

29	subroutine ttt

30	implicit none

31	integer, parameter :: ITMAX = 100   !最大计算次数

32	integer, parameter :: NROOT = 1     !初始值数量

33	real, parameter :: EPS = 1.0E-6     !定义两个解之间最短距离，任何两个解的差值一定不小于EPS,ABS(Xi-Xj)>=eps

34	real, parameter :: ERRABS = 1.0E-6  !当ABS(F(Xi))<ERRABS时，就视Xi为F(Xi)的解

35	real, parameter :: ERRREL = 1.0E-6  !当新Xi和就Xi值，很接近时，就视Xi为解，不再做新的逼近

36	real, parameter :: ETA =0.01 !0.0005 !求解步长

37	integer INFO(NROOT)

38	real,external :: f

39	real :: X(NROOT)     !答案

40	real :: XGUESS(NROOT) =0 !(/ -1.0000, 1.0000/)

41

42	call ZREAL (f, ERRABS, ERRREL, EPS, ETA, &

43	NROOT, ITMAX, XGUESS, X, INFO)

44	write(*,*) X !,f(X)

45

46	end subroutine

47

48	end module

风和

该问题已解决