 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 829 | 回復(fù): 3 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||
caoweiyue新蟲 (小有名氣)
|
[求助]
剛開始學(xué)Fortran 幫忙看一下這個程序 為什么算不出結(jié)果
|
|
|
下面的程序算的是全選主元的高斯消去法， 我是照徐士良的算法集編的，但是結(jié)果都是0 program main implicit none integer::i,L integer,parameter::n=4 real::A(n,n),b(n),x(n) 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 / call agaus(A,b,n,L) if(L==0) then write(*,*) "fail" else do i=1,n write(*,"('x',i1,'=',f6.4)"  i,x(i)end do end if stop end program subroutine agaus(A,b,n,L) implicit none integer:: i,j,k,n,is,js,L real:: D, T real:: A(n,n),b(n),x(n) L=1 do k=1,n-1 D=0.0 !選出了最大的，并保存其位置 do i=k,n do j=k,n if( abs(A(i,j))>D ) then D=A(i,j) is=i js=j end if end do end do if(D+1.0==1.0) then L=0 return end if !列交換 if(js/=k) then do i=1,n T=A(i,k) A(i,k)=A(i,js) A(i,js)=T end do end if !行交換 if(is/=k) then do j=k,n T=A(k,j) A(k,j)=A(is,j) A(is,j)=A(k,j) end do T=b(k) b(k)=b(is) b(is)=b(k) end if !化成三角矩陣 do i=k+1,n do j=k,n A(i,j)=A(i,j)-A(i,k)/A(k,k)*A(k,j) end do b(i)=b(i)-A(i,k)/A(k,k)*b(k) end do end do if(A(n,n)+1.0==1.0) then L=0 return end if !回代 x(n)=b(n)/A(n,n) do i=n-1,1,-1 T=0.0 do j=i+1,n T=T+A(i,j)*x(j) end do x(i)=(b(i)-T)/A(i,i) end do return end subroutine agaus |
回復(fù)此樓» 猜你喜歡
0854電子信息求調(diào)劑
已經(jīng)有7人回復(fù)
-
0854人工智能方向招收調(diào)劑
已經(jīng)有3人回復(fù)
-
調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿北化315 求調(diào)劑
已經(jīng)有3人回復(fù)
-
有機(jī)合成求助
已經(jīng)有5人回復(fù)
-
材料學(xué)碩333求調(diào)劑
已經(jīng)有3人回復(fù)
-
289求調(diào)劑
已經(jīng)有7人回復(fù)
-
收08調(diào)劑生
已經(jīng)有6人回復(fù)
-
材料專碩找調(diào)劑
已經(jīng)有5人回復(fù)
-
085602 289分求調(diào)劑
已經(jīng)有5人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 配置一臺1萬左右的計算機(jī)跑Fortran程序以及高斯的計算
已經(jīng)有8人回復(fù)
- 求高人幫忙看一下我這個Fortran程序的錯誤出在哪里？
已經(jīng)有12人回復(fù)
- 物理模擬是用c語言好還是fortran好？
已經(jīng)有30人回復(fù)
- 我編的Simpson積分法fortran程序給不出結(jié)果，大俠們看看哪里出了問題？
已經(jīng)有4人回復(fù)
- Fortran 95／2003科學(xué)計算與工程 對應(yīng)的各章程序源代碼
已經(jīng)有16人回復(fù)
- 有關(guān)fortran的一次而問題，希望大家能幫幫忙，謝謝
已經(jīng)有4人回復(fù)
- 寫了一個fortran90的小程序，編譯通不過，請大俠幫忙
已經(jīng)有59人回復(fù)
- 【求助】用fortran怎么描述下面這樣的情況，我是新手，請大家?guī)兔?/a>
已經(jīng)有8人回復(fù)
- 【求助】初學(xué)fortran，
已經(jīng)有15人回復(fù)
- 【求助】請問一下fortran的可視化編程
已經(jīng)有7人回復(fù)
- 【原創(chuàng)】elk中的狀態(tài)方程擬合工具eos
已經(jīng)有128人回復(fù)
- 【求助】fortran計算時出現(xiàn)NaN
已經(jīng)有11人回復(fù)
daizhijuan
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 12.5
- 帖子: 1
- 在線: 27分鐘
- 蟲號: 2184351
- 注冊: 2012-12-13
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
4樓2012-12-13 16:52:02
Jasonluo
至尊木蟲 (知名作家)
中央書籍處書記
- 應(yīng)助: 61 (初中生)
- 貴賓: 0.01
- 金幣: 22628.1
- 散金: 166
- 紅花: 14
- 帖子: 9258
- 在線: 680.9小時
- 蟲號: 518969
- 注冊: 2008-03-05
- 性別: GG
- 專業(yè): 自然地理學(xué)
2樓2012-05-16 20:59:37
jwwang_lq
鐵蟲 (初入文壇)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 17
- 帖子: 7
- 在線: 1.6小時
- 蟲號: 1719996
- 注冊: 2012-03-27
- 性別: GG
- 專業(yè): 生物大分子結(jié)構(gòu)與功能
3樓2012-05-17 08:51:20
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研]
 材料與化工考研調(diào)劑
+6
|
孅華 2026-03-22 | 6/300 |
|
|---|---|---|---|---|
|
|
[考研] 【雙一流院校新能源、環(huán)境材料，材料加工與模擬招收大量調(diào)劑】 +4 | Higraduate 2026-03-22 | 7/350 |
|
|
[考研] 求材料，環(huán)境專業(yè)調(diào)劑 +3 | 18567500178 2026-03-18 | 3/150 |
|
|
|
[考研] 0703化學(xué)求調(diào)劑 +4 | 奶油草莓. 2026-03-22 | 5/250 |
|
|
|
[考研] 263求調(diào)劑 +6 | yqdszhdap－ 2026-03-22 | 9/450 |
|
|
|
[考研] 291求調(diào)劑 +5 | 孅華 2026-03-22 | 5/250 |
|
|
|
[考研] 323求調(diào)劑 +6 | 洼小桶 2026-03-18 | 6/300 |
|
|
|
[考研] 307求調(diào)劑 +11 | 冷笙123 2026-03-17 | 11/550 |
|
|
|
[考研] 求調(diào)劑一志愿海大，0703化學(xué)學(xué)碩304分，有大創(chuàng)項目，四級已過 +6 | 幸運(yùn)哩哩 2026-03-22 | 10/500 |
|
|
|
[考研] 308求調(diào)劑 +3 | 墨墨漠 2026-03-21 | 3/150 |
|
|
|
[考研] 275求調(diào)劑 +6 | shansx 2026-03-22 | 8/400 |
|
|
|
[考研] 考研調(diào)劑 +4 | 來好運(yùn)來來來 2026-03-21 | 4/200 |
|
|
|
[考研] 085600材料與化工306 +4 | z1z2z3879 2026-03-21 | 4/200 |
|
|
|
[考研] 化學(xué)調(diào)劑 +5 | yzysaa 2026-03-21 | 5/250 |
|
|
|
[考研] 材料 271求調(diào)劑 +5 | 展信悅_ 2026-03-21 | 5/250 |
|
|
|
[基金申請]
學(xué)校已經(jīng)提交到NSFC，還能修改嗎？
 40+4
|
babangida 2026-03-19 | 9/450 |
|
|
|
[考研] 材料與化工（0856）304求 B區(qū) 調(diào)劑 +3 | 邱gl 2026-03-21 | 3/150 |
|
|
|
[考研] 一志愿華南師大 070300（化學(xué)）304分求調(diào)劑 +3 | 0703武芊慧雪304 2026-03-18 | 3/150 |
|
|
|
[考研] 295材料求調(diào)劑，一志愿武漢理工085601專碩 +5 | Charlieyq 2026-03-19 | 5/250 |
|
|
|
[考研]
320求調(diào)劑0856
+3
|
不想起名字112 2026-03-19 | 3/150 |
|