Fortran Coder

查看: 19897|回复: 19
打印 上一主题 下一主题

[数学库] 用imsl求非线性方程组时运行出错,显示结果如下

[复制链接]

8

帖子

1

主题

0

精华

入门

F 币
39 元
贡献
11 点
楼主
发表于 2018-5-27 22:51:50 | 显示全部楼层
我也碰到同样的问题。

8

帖子

1

主题

0

精华

入门

F 币
39 元
贡献
11 点
沙发
发表于 2018-5-27 22:52:11 | 显示全部楼层
我在用 compaq visual fortran 6.6

8

帖子

1

主题

0

精华

入门

F 币
39 元
贡献
11 点
板凳
发表于 2018-5-27 22:52:33 | 显示全部楼层
Error: This is an actual argument keyword name, and not a dummy argument name.   [XGUESS]
call neqnf(fcn, x, XGUESS=xguess, FNORM=fnorm)

8

帖子

1

主题

0

精华

入门

F 币
39 元
贡献
11 点
地板
发表于 2018-5-27 23:16:08 | 显示全部楼层
program hello

    use imsl

        implicit none

        integer n
        parameter (n = 3)

        integer k, nout
        real fnorm, x(n), xguess(n)
        external fcn

data xguess/4.0, 4.0, 4.0/

call umach (2,nout)
call neqnf(fcn, x, XGUESS=xguess, FNORM=fnorm)

write (nout, 99999) (x(k),k=1,n), fnorm
99999 format (' the solution to the system is', /, ' x = (', 3f5.1, &
')', /, ' with fnorm = ', f5.4, //)


        end program hello

        subroutine fcn(x,f,n)
        integer n
        real x(n), f(n)
        real exp,sin
        intrinsic exp, sin

        f(1) = x(1)+exp(x(1)-1.0)+(x(2)+x(3))*(x(2)+x(3)) -27.0
        f(2) = exp(x(2)-2.0)/x(1)+x(3)*x(3)-10.0
        f(3) = x(3)+sin(x(2)-2.0)+x(2)*x(2)-7.0

        return
        end subroutine fcn
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

捐赠本站|Archiver|关于我们 About Us|小黑屋|Fcode ( 京ICP备18005632-2号 )

GMT+8, 2024-5-3 12:10

Powered by Tencent X3.4

© 2013-2024 Tencent

快速回复 返回顶部 返回列表