应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 【求助】急需一个程序
91/1返回列表
查看: 1446  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

飞扬123

新虫 (初入文坛)

[交流] 【求助】急需一个程序

急需：一个求解一元多次方程的Fortran程序，谢谢！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼 2010-11-30 09:57:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)
自已写，或找现成的软件做……
赞一下 回复此楼
2楼2010-11-30 10:05:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gsview

木虫 (小有名气)
★
余泽成(金币+1):谢谢参与应助，欢迎常来程序语言版！ 2010-12-24 16:15:28
从书上给你找的！






        DIMENSION A(4,4),B(4),X(4),JS(4)
        DOUBLE PRECISION A,B,X
        DATA A/0.2368,0.1968,0.1582,1.1161,0.2471,0.2071,1.1675,0.1254
     *       ,0.2568,1.2168,0.1768,0.1397,1.2671,0.2271,0.1871,0.1490/
        DATA B/1.8471,1.7471,1.6471,1.5471/
        N=4
        CALL AGAUS(A,B,N,X,L,JS)
        IF (L.NE.0) THEN
          WRITE(*,10) (I,X(I),I=1,4)
        END IF
10        FORMAT(1X,'X(',I2,' )=',D15.6)
        END




        SUBROUTINE AGAUS(A,B,N,X,L,JS)
        DIMENSION A(N,N),X(N),B(N),JS(N)
        DOUBLE PRECISION A,B,X,T
        L=1
        DO 50 K=1,N-1
          D=0.0
          DO 210 I=K,N
          DO 210 J=K,N
            IF (ABS(A(I,J)).GT.D) THEN
              D=ABS(A(I,J))
              JS(K)=J
              IS=I
            END IF
210          CONTINUE
          IF (D+1.0.EQ.1.0) THEN
            L=0
          ELSE
            IF (JS(K).NE.K) THEN
              DO 220 I=1,N
                T=A(I,K)
                A(I,K)=A(I,JS(K))
                A(I,JS(K))=T
220              CONTINUE
            END IF
            IF (IS.NE.K) THEN
              DO 230 J=K,N
                T=A(K,J)
                A(K,J)=A(IS,J)
                A(IS,J)=T
230              CONTINUE
              T=B(K)
              B(K)=B(IS)
              B(IS)=T
            END IF
          END IF
          IF (L.EQ.0) THEN
            WRITE(*,100)
            RETURN
          END IF
          DO 10 J=K+1,N
            A(K,J)=A(K,J)/A(K,K)
10          CONTINUE
          B(K)=B(K)/A(K,K)
          DO 30 I=K+1,N
            DO 20 J=K+1,N
              A(I,J)=A(I,J)-A(I,K)*A(K,J)
20            CONTINUE
            B(I)=B(I)-A(I,K)*B(K)
30          CONTINUE
50        CONTINUE
        IF (ABS(A(N,N))+1.0.EQ.1.0) THEN
          L=0
          WRITE(*,100)
          RETURN
        END IF
        X(N)=B(N)/A(N,N)
        DO 70 I=N-1,1,-1
          T=0.0
          DO 60 J=I+1,N
            T=T+A(I,J)*X(J)
60          CONTINUE
          X(I)=B(I)-T
70        CONTINUE
100        FORMAT(1X,' FAIL ')
        JS(N)=N
        DO 150 K=N,1,-1
          IF (JS(K).NE.K) THEN
            T=X(K)
            X(K)=X(JS(K))
            X(JS(K))=T
          END IF
150        CONTINUE
        RETURN
        END
赞一下(1人) 回复此楼
3楼2010-11-30 16:41:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)
★
nono2009(金币+1):呵呵 2010-12-02 07:46:04
人家求的是一元多次方程，不是线性方程组，呵呵……

一元多次方程一般都是通过迭代的方法求解的，效率之间的差别会很大，但查一下书也应该都会有的。不管怎么样，这种事情还是应该自己动手……
赞一下(1人) 回复此楼
4楼2010-11-30 21:52:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gsview

木虫 (小有名气)
不好意思，看错了！呵呵！
赞一下 回复此楼
5楼2010-12-01 13:13:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

553宿舍

新虫 (初入文坛)
呵呵，挺有意思的！
回复此楼
6楼2010-12-02 16:44:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangxn06

金虫 (正式写手)
彭国伦的书里面有吧，数值分析教材中也应该有具体的方法，用程序实现起来也不难
赞一下 回复此楼
7楼2010-12-17 18:47:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

klinedancer

铜虫 (小有名气)
★
余泽成(金币+1):谢谢参与应助，欢迎常来程序语言版！ 2010-12-24 16:15:45
一元3次以上的方程没有公式解，常用的方法有二分法，牛顿迭代法，还是要看你主要采用什么算法，程序不是问题，算法固定后后容易实现！
赞一下(1人) 回复此楼
8楼2010-12-17 20:41:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luoxitc

新虫 (初入文坛)
★
余泽成(金币+1):谢谢参与应助，欢迎常来程序语言版！ 2010-12-24 16:16:25
就用牛顿迭代法
一般书上有现成的
赞一下(1人) 回复此楼
9楼2010-12-22 14:12:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 飞扬123 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +31 医学老男孩 2026-08-13 69/3450 2026-08-16 00:16 by mo_mo_mo_mo
[基金申请] 有时候，自然基金真的不能太认真 （我的申报经验） +9 majunge000 2026-08-11 11/550 2026-08-15 23:20 by botar
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +7 Ldrop2023 2026-08-13 7/350 2026-08-15 15:57 by mzhh2000
[基金申请] 各位道友，我要去昆明玩几天，回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 22/1100 2026-08-14 23:45 by Noways
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +23 Lanmanbaby 2026-08-09 37/1850 2026-08-14 13:40 by Equinoxhua
[硕博家园] 请教兼职经验 +3 是阿文鸭 2026-08-09 3/150 2026-08-14 12:07 by HER12025
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] FileCode能看出啥？ +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[基金申请] 我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人 +10 家与远方 2026-08-10 15/750 2026-08-14 02:08 by 绵羊哥哥
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[基金申请] 结合人工智能，周易传统文化，filecode打分制来了，3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见，大家看看几个A，几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分？ +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 确定了，国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
信息提示
关闭
请填处理意见
关闭