版塊導航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

飛揚123

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 301.6
帖子: 3
在線: 46分鐘
蟲號: 1151311

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

急需：一個求解一元多次方程的Fortran程序，謝謝！

回復此樓

» 猜你喜歡

0703化學336分求調(diào)劑已經(jīng)有6人回復
304求調(diào)劑已經(jīng)有9人回復
293求調(diào)劑已經(jīng)有13人回復
281求調(diào)劑（0805）已經(jīng)有8人回復
環(huán)境領(lǐng)域全國重點實驗室招收博士1-2名已經(jīng)有3人回復
材料專碩306英一數(shù)二已經(jīng)有10人回復
301求調(diào)劑已經(jīng)有6人回復
一志愿天津大學化學工藝專業(yè)（081702）315分求調(diào)劑已經(jīng)有7人回復
302求調(diào)劑已經(jīng)有6人回復
26博士申請已經(jīng)有3人回復

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

需要用matlab擬合數(shù)據(jù)的外行人急需幫助已經(jīng)有10人回復
MATLAB 編程求解普通實驗數(shù)據(jù)時間序列的最大李雅普諾夫指數(shù)！��！已經(jīng)有12人回復
兩本經(jīng)典的入門教程C Primer plus（第五版）和C++ Primer plus（第五版）已經(jīng)有58人回復
求助~~!!請懂SQL SERVER的進. 已經(jīng)有13人回復
【求助】急需最新全國行政區(qū)劃數(shù)據(jù)和常見漢字及其拼音（也可以是聲母的第一個字母）數(shù) 已經(jīng)有4人回復
【求助】用matlab編程，需要將數(shù)據(jù)自檢分類（高手進來指點下）已經(jīng)有5人回復
【求助】購買程序降溫的烘箱求助��！已經(jīng)有11人回復
【討論】哪位有VB6.0英文版啊，急需【已完成】已經(jīng)有3人回復
【分享】Origin 8.x 視頻教程（75個）已經(jīng)有177人回復

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

1樓 2010-11-30 09:57:33

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

snoopyzhao

至尊木蟲 (職業(yè)作家)

自已寫，或找現(xiàn)成的軟件做……

贊一下

回復此樓

2樓2010-11-30 10:05:25

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

gsview

木蟲 (小有名氣)

應(yīng)助: 4 (幼兒園)
金幣: 2446.9
帖子: 146
在線: 234.2小時
蟲號: 443276

★
余澤成(金幣+1):謝謝參與應(yīng)助，歡迎常來程序語言版！ 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

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

snoopyzhao

至尊木蟲 (職業(yè)作家)

★
nono2009(金幣+1):呵呵 2010-12-02 07:46:04

人家求的是一元多次方程，不是線性方程組，呵呵……

一元多次方程一般都是通過迭代的方法求解的，效率之間的差別會很大，但查一下書也應(yīng)該都會有的。不管怎么樣，這種事情還是應(yīng)該自己動手……

贊一下(1人)

回復此樓

4樓2010-11-30 21:52:48

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

gsview

木蟲 (小有名氣)

應(yīng)助: 4 (幼兒園)
金幣: 2446.9
帖子: 146
在線: 234.2小時
蟲號: 443276

不好意思，看錯了！呵呵！

贊一下

回復此樓

5樓2010-12-01 13:13:20

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

553宿舍

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 74.5
帖子: 25
在線: 9.8小時
蟲號: 1002314

呵呵，挺有意思的！

回復此樓

6樓2010-12-02 16:44:06

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wangxn06

金蟲 (正式寫手)

應(yīng)助: 2 (幼兒園)
金幣: 576.7
帖子: 509
在線: 250.5小時
蟲號: 1156835

彭國倫的書里面有吧，數(shù)值分析教材中也應(yīng)該有具體的方法，用程序?qū)崿F(xiàn)起來也不難

贊一下

回復此樓

7樓2010-12-17 18:47:11

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

klinedancer

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 911.7
帖子: 300
在線: 98.7小時
蟲號: 1066832

★
余澤成(金幣+1):謝謝參與應(yīng)助，歡迎常來程序語言版！ 2010-12-24 16:15:45

一元3次以上的方程沒有公式解，常用的方法有二分法，牛頓迭代法，還是要看你主要采用什么算法，程序不是問題，算法固定后后容易實現(xiàn)！

贊一下(1人)

回復此樓

8樓2010-12-17 20:41:58

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

luoxitc

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 1103.9
帖子: 37
在線: 7.6小時
蟲號: 635873

★
余澤成(金幣+1):謝謝參與應(yīng)助，歡迎常來程序語言版！ 2010-12-24 16:16:25

就用牛頓迭代法
一般書上有現(xiàn)成的

贊一下(1人)

回復此樓

9樓2010-12-22 14:12:46

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主飛揚123 的主題更新

返回列表

普通表情龍兔虎貓高級回復 (可上傳附件)

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 293求調(diào)劑 +8	zjl的號 2026-03-16	13/650	2026-03-18 04:32 by anny19840123
[考研] 302求調(diào)劑 +6	呼呼呼。。。。 2026-03-17	6/300	2026-03-17 23:56 by 星空星月
[考研] 化學工程321分求調(diào)劑 +11	大米飯！ 2026-03-15	14/700	2026-03-17 17:11 by ruiyingmiao
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6	薛云鵬 2026-03-15	6/300	2026-03-17 11:05 by 學員h26Tkc
[碩博家園] 深圳大學碩士招生（2026秋，傳感器方向，僅錄取第一志愿） +4	xujiaoszu 2026-03-11	9/450	2026-03-17 10:29 by xujiaoszu
[基金申請] 國自科面上基金字體 +6	iwuli 2026-03-12	7/350	2026-03-16 21:18 by sculhf
[考研] 一志愿211 0703方向310分求調(diào)劑 +3	努力奮斗112 2026-03-15	3/150	2026-03-16 16:44 by houyaoxu
[考研] 304求調(diào)劑 +3	曼殊2266 2026-03-14	3/150	2026-03-16 16:39 by houyaoxu
[考研] 0703 物理化學調(diào)劑 +3	我可以上岸的對?/a> 2026-03-13	5/250	2026-03-16 10:50 by 我可以上岸的對?/a>
[教師之家] 焦慮 +7	水冰月月野兔 2026-03-13	9/450	2026-03-16 10:00 by Quakerbird
[考研] 326求調(diào)劑 +3	mlpqaz03 2026-03-15	3/150	2026-03-16 07:33 by Iveryant
[考研] 0856專碩279求調(diào)劑 +5	加油加油！? 2026-03-15	5/250	2026-03-15 11:58 by 2020015
[考研] 265求調(diào)劑 +4	威化餅07 2026-03-12	4/200	2026-03-14 17:23 by userper
[考研] 一志愿哈工大材料324分求調(diào)劑 +5	閆旭東 2026-03-14	5/250	2026-03-14 14:53 by 木瓜膏
[考研] 材料與化工求調(diào)劑一志愿 985 總分 295 +8	dream…… 2026-03-12	8/400	2026-03-13 22:17 by 星空星月
[考研] 求材料調(diào)劑 085600英一數(shù)二總分302 前三科235 精通機器學習一志愿哈工大 +4	林yaxin 2026-03-12	4/200	2026-03-13 22:04 by 星空星月
[考研] 311求調(diào)劑 +3	冬十三 2026-03-13	3/150	2026-03-13 20:41 by JourneyLucky
[考研] 工科材料085601 279求調(diào)劑 +8	困于星晨 2026-03-12	10/500	2026-03-13 15:42 by ms629
[考研] 314求調(diào)劑 +7	無懈可擊的巨人 2026-03-12	7/350	2026-03-13 15:40 by JourneyLucky
[考研] 求調(diào)劑 +3	程雨杭 2026-03-12	3/150	2026-03-13 15:06 by JourneyLucky

24小時熱門版塊排行榜

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

» 猜你喜歡

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

» 搶金幣啦！回帖就可以得到: 查看全部散金貼

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

» 搶金幣啦！回帖就可以得到:

查看全部散金貼