用imsl求非线性方程组时运行出错，显示结果如下

fcode

第一个参数 f 应该是一个函数，而不是数组。

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

program kakusann
  Include 'link_f90_static.h'
  USE NEQNF_INT
  USE UMACH_INT
  implicit none
  integer,parameter::n=2
  integer::nout
  real::x(n),xguess(n),fnorm
  external fcn
  xguess=3.0
  call umach(2,nout)
  CALL NEQNF (Fcn, X, xguess=xguess, fnorm=fnorm)
  WRITE (NOUT,*) x , FNORM
end program kakusann

SUBROUTINE FCN (X, F, N)
  INTEGER N
  REAL X(N), F(N)
  F(1)=55*x(1)**3+2*x(1)
  F(2)=4*x(2)**3+2*x(2)
END SUBROUTINE FCN