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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 【求助】急需一个程序
91/1返回列表
查看: 1094  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 304求调剂 +3 曼殊2266 2026-03-18 3/150 2026-03-19 14:42 by peike
[考研] 求调剂，一志愿:南京航空航天大学大学 ，080500材料科学与工程学硕，总分289分 +3 @taotao 2026-03-19 3/150 2026-03-19 14:07 by peike
[考研] 294求调剂材料与化工专硕 +11 陌の森林 2026-03-18 11/550 2026-03-19 13:22 by houyaoxu
[考研] 287求调剂 +3 晨昏线与星海 2026-03-19 4/200 2026-03-19 12:32 by peike
[考研] 085601材料工程专硕求调剂 +8 慕寒mio 2026-03-16 8/400 2026-03-19 12:30 by peike
[考研] 能源材料化学课题组招收硕士研究生8-10名 +4 脱颖而出 2026-03-16 11/550 2026-03-19 08:11 by 脱颖而出
[考研] 085410人工智能专硕317求调剂（0854都可以） +3 xbxudjdn 2026-03-18 3/150 2026-03-18 22:14 by zhq0425
[考研] 311求调剂 +4 冬十三 2026-03-18 4/200 2026-03-18 21:47 by 尽舜尧1
[考研] 085601专硕，总分342求调剂，地区不限 +5 share_joy 2026-03-16 5/250 2026-03-18 14:48 by haxia
[考研] 收复试调剂生 +4 雨后秋荷 2026-03-18 4/200 2026-03-18 14:16 by elevennnne
[考研] 070300化学319求调剂 +6 锦鲤0909 2026-03-17 6/300 2026-03-18 13:22 by Iveryant
[考研] 288求调剂，一志愿华南理工大学071005 +4 ioodiiij 2026-03-17 4/200 2026-03-18 12:36 by Linda Hu
[考研] 工科材料085601 279求调剂 +6 困于星晨 2026-03-17 6/300 2026-03-18 10:21 by kkcoco25
[考研] 265求调剂 +3 梁梁校校 2026-03-17 3/150 2026-03-18 09:12 by zhukairuo
[考研] 有没有道铁/土木的想调剂南林，给自己招师弟中～ +3 TqlXswl 2026-03-16 7/350 2026-03-17 15:23 by TqlXswl
[考研] 211本，11408一志愿中科院277分，曾在中科院自动化所实习 +6 Losir 2026-03-12 7/350 2026-03-17 12:09 by danranxie
[考研] 318求调剂 +3 Yanyali 2026-03-15 3/150 2026-03-16 16:41 by houyaoxu
[考研] 中科院材料273求调剂 +4 yzydy 2026-03-15 4/200 2026-03-16 15:59 by Gaodh_82
[考研] 326求调剂 +3 mlpqaz03 2026-03-15 3/150 2026-03-16 07:33 by Iveryant
[考研] 321求调剂（食品/专硕） +3 mxcz321 2026-03-12 6/300 2026-03-13 08:45 by xc321
信息提示
关闭
请填处理意见
关闭