求助，IMSL库调用失败 - Fortran Coder

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


include 'link_fnl_shared.h'

use lin_sol_gen_int

use rand_gen_int

use error_option_packet

  integer, parameter :: NDATA = 10, N = 20

  real XDATA(NDATA), FDATA(NDATA)

  real XVEC(N), VALUE(N)

  real, parameter :: xmin = -5.0, xmax = 5.0

  real xinc, xp

  integer i



  xinc = (xmax-xmin)/(NDATA-1)

  xp = xmin

  do i=1, NDATA

    XDATA(I) = xp

    FDATA(I) = SIN( XDATA(I) )

        xp = xp+xinc

  end do



  xinc = (xmax-xmin)/(N-1)

  xp = xmin

  do i=1, N

    XVEC(I) = xp

        xp = xp+xinc

  end do

  ! 做插值

  call CSIEZ (NDATA, XDATA, FDATA, N, XVEC, VALUE)

  ! 输出插值结果

  do i=1, N

    write(*,"('(',F5.2,',',F5.2,') error:'F6.3)") &

        XVEC(i), VALUE(i), VALUE(i)-sin( XVEC(i) )

  end do

  stop

end