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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 【求助】急需一个程序
91/1返回列表
查看: 1125  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 复试调剂，一志愿郑州大学材料与化工289分 +14 硕星赴 2026-04-08 14/700 2026-04-08 23:59 by GouQ
[考研] 22408 270分 +4 sanjin020722 2026-04-08 5/250 2026-04-08 23:57 by GouQ
[考研] 296求调剂 +3 汪！？！ 2026-04-08 3/150 2026-04-08 21:46 by 土木硕士招生
[考研] 材料专硕(0856) 339分求调剂 +16 哈哈哈鹅哈哈哈 2026-04-05 16/800 2026-04-08 16:02 by luoyongfeng
[考研] 考研调剂 +11 冰冰，，， 2026-04-07 11/550 2026-04-08 15:08 by Delta2012
[考研] 求考研材料调剂 +3 材化李可 2026-04-07 3/150 2026-04-08 00:21 by JourneyLucky
[考研] 315求调剂 +17 小羊小羊_ 2026-04-02 18/900 2026-04-07 22:01 by lijunpoly
[考研] 生物工程求调剂 +13 喜欢还是不甘心 2026-04-05 13/650 2026-04-07 16:55 by Ecowxq666！
[考研] 081700学硕，323分，一志愿中国海洋大学求调剂学校 +19 披星河 2026-04-04 19/950 2026-04-07 15:00 by 上岸快快
[考研] 软工学硕299求调剂 +6 useryy 2026-04-07 6/300 2026-04-07 09:50 by vgtyfty
[考研] 考研调剂 +3 Wwwwwww哇 2026-04-06 3/150 2026-04-06 20:55 by lbsjt
[考研] 081200-11408-276学硕求调剂 +5 崔wj 2026-04-05 5/250 2026-04-06 15:40 by lin-da
[考研] 0703求调剂383分 +9 W55j 2026-04-03 9/450 2026-04-06 06:50 by houyaoxu
[考研] 材料调剂 +6 一样YWY 2026-04-05 6/300 2026-04-05 20:30 by 南航~万老师
[考研] 0854求调剂 +4 assdll 2026-04-04 4/200 2026-04-05 09:44 by zhq0425
[考研] 材料与化工306分找调剂 +23 沧海轻舟e 2026-04-02 27/1350 2026-04-04 21:52 by laoshidan
[考研] 297求调剂 +11 ljy20040718！ 2026-04-03 13/650 2026-04-04 09:23 by 来看流星雨10
[考研] 315求调剂 +6 顺理成张 2026-04-03 8/400 2026-04-03 14:04 by 百灵童888
[考研] 数一英一285求调剂 +7 AZMK 2026-04-03 9/450 2026-04-03 13:03 by ms629
[考研] 298求B区调剂 +4 zzz，，r 2026-04-02 5/250 2026-04-02 12:17 by 土木硕士招生
信息提示
关闭
请填处理意见
关闭