24小时热门版块排行榜    

查看: 1093  |  回复: 8

飞扬123

新虫 (初入文坛)


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

急需:一个求解一元多次方程的Fortran程序,谢谢!
回复此楼

» 猜你喜欢

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

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

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
自已写,或找现成的软件做……
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
3楼2010-11-30 16:41:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nono2009(金币+1):呵呵 2010-12-02 07:46:04
人家求的是一元多次方程,不是线性方程组,呵呵……

一元多次方程一般都是通过迭代的方法求解的,效率之间的差别会很大,但查一下书也应该都会有的。不管怎么样,这种事情还是应该自己动手……
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次以上的方程没有公式解,常用的方法有二分法,牛顿迭代法,还是要看你主要采用什么算法,程序不是问题,算法固定后后容易实现!
8楼2010-12-17 20:41:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luoxitc

新虫 (初入文坛)



余泽成(金币+1):谢谢参与应助,欢迎常来程序语言版! 2010-12-24 16:16:25
就用牛顿迭代法
一般书上有现成的
9楼2010-12-22 14:12:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 飞扬123 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 288求调剂 +6 于海海海海 2026-03-19 6/300 2026-03-19 12:33 by peike
[考研] 招收调剂硕士 +3 lidianxing 2026-03-19 7/350 2026-03-19 12:31 by 山东小li
[考研] 一志愿吉林大学材料学硕321求调剂 +6 Ymlll 2026-03-18 9/450 2026-03-19 10:28 by 星空星月
[考研] 一志愿985,本科211,0817化学工程与技术319求调剂 +10 Liwangman 2026-03-15 10/500 2026-03-19 10:25 by 无际的草原
[教师之家] 焦虑 +9 水冰月月野兔 2026-03-13 13/650 2026-03-19 09:50 by otani
[考研] 本科郑州大学物理学院,一志愿华科070200学硕,346求调剂 +4 我不是一根葱 2026-03-18 4/200 2026-03-19 09:11 by 浮云166
[考研] 304求调剂 +6 司空. 2026-03-18 6/300 2026-03-18 23:03 by 星空星月
[考研] 材料与化工一志愿南昌大学327求调剂推荐 +8 Ncdx123456 2026-03-13 9/450 2026-03-18 14:40 by haxia
[考研] 298-一志愿中国农业大学-求调剂 +7 手机用户 2026-03-17 7/350 2026-03-18 14:34 by vgtyfty
[考研] 331求调剂(0703有机化学 +7 ZY-05 2026-03-13 8/400 2026-03-18 14:13 by 007_lilei
[考研] 0703化学调剂 +4 pupcoco 2026-03-17 7/350 2026-03-18 12:14 by djl2006
[考研] 307求调剂 +3 冷笙123 2026-03-17 3/150 2026-03-18 09:55 by macy2011
[考研] 268求调剂 +8 一定有学上- 2026-03-14 9/450 2026-03-17 17:47 by laoshidan
[考博] 26申博 +4 八6八68 2026-03-16 4/200 2026-03-17 13:00 by 轻松不少随
[考研] [导师推荐]西南科技大学国防/材料导师推荐 +3 尖角小荷 2026-03-16 6/300 2026-03-16 23:21 by 尖角小荷
[考研] 321求调剂 +5 大米饭! 2026-03-15 5/250 2026-03-16 16:33 by houyaoxu
[考研] 070300化学学硕求调剂 +6 太想进步了0608 2026-03-16 6/300 2026-03-16 16:13 by kykm678
[考研] 085600材料与化工 求调剂 +13 enenenhui 2026-03-13 14/700 2026-03-16 15:19 by 了了了了。。
[考研] 招收0805(材料)调剂 +3 18595523086 2026-03-13 3/150 2026-03-14 00:33 by 123%、
[考研] 295求调剂 +3 小匕仔汁 2026-03-12 3/150 2026-03-13 15:17 by vgtyfty
信息提示
请填处理意见