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

风和

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

风和

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

program main
include 'link_f90_static.h'
use test

use numerical_libraries

call ttt


pause
end program



module test
include 'link_f90_static.h'
use numerical_libraries
contains

real function f(x)  ！样例，实际主计算程序比这个多的多
  implicit none
  real x
  f=X**2+5*X+3
!  write(*,*)x,f
  return
end function


subroutine ttt
  implicit none
  integer, parameter :: ITMAX = 100   !最大计算次数
  integer, parameter :: NROOT = 1     !初始值数量
  real, parameter :: EPS = 1.0E-6     !定义两个解之间最短距离，任何两个解的差值一定不小于EPS,ABS(Xi-Xj)>=eps
  real, parameter :: ERRABS = 1.0E-6  !当ABS(F(Xi))<ERRABS时，就视Xi为F(Xi)的解
  real, parameter :: ERRREL = 1.0E-6  !当新Xi和就Xi值，很接近时，就视Xi为解，不再做新的逼近
  real, parameter :: ETA =0.01 !0.0005 !求解步长
  integer INFO(NROOT)
  real,external :: f
  real :: X(NROOT)     !答案
  real :: XGUESS(NROOT) =0 !(/ -1.0000, 1.0000/)
  
  call ZREAL (f, ERRABS, ERRREL, EPS, ETA, &
              NROOT, ITMAX, XGUESS, X, INFO)
  write(*,*) X !,f(X)

end subroutine

end module

风和

该问题已解决