 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 820 | 回復(fù): 3 | ||
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ù)此樓» 猜你喜歡
求調(diào)劑
已經(jīng)有4人回復(fù)
-
求調(diào)劑
已經(jīng)有4人回復(fù)
-
277材料科學(xué)與工程080500求調(diào)劑
已經(jīng)有4人回復(fù)
-
學(xué)校已經(jīng)提交到NSFC，還能修改嗎？
已經(jīng)有9人回復(fù)
-
286分人工智能專業(yè)請求調(diào)劑愿意跨考！
已經(jīng)有7人回復(fù)
-
286求調(diào)劑
已經(jīng)有6人回復(fù)
-
280求調(diào)劑
已經(jīng)有10人回復(fù)
-
材料學(xué)碩301分求調(diào)劑
已經(jīng)有4人回復(fù)
-
278求調(diào)劑
已經(jīng)有8人回復(fù)
-
求助
已經(jīng)有5人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 配置一臺1萬左右的計算機跑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ù)
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
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
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 caoweiyue 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085601調(diào)劑 358分 +3 | zzzzggh 2026-03-20 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 材料 271求調(diào)劑 +3 | 展信悅_ 2026-03-21 | 3/150 |
|
|
|
[考研] 070300化學(xué)319求調(diào)劑 +7 | 錦鯉0909 2026-03-17 | 7/350 |
|
|
|
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +4 | 生物工程調(diào)劑 2026-03-17 | 8/400 |
|
|
|
[考研] 一志愿武理材料305分求調(diào)劑 +6 | 想上岸的鯉魚 2026-03-18 | 7/350 |
|
|
|
[考研] 274求調(diào)劑 +10 | S.H1 2026-03-18 | 10/500 |
|
|
|
[考研] 295求調(diào)劑 +4 | 一志愿京區(qū)211 2026-03-18 | 6/300 |
|
|
|
[考研] 294求調(diào)劑材料與化工專碩 +15 | 陌の森林 2026-03-18 | 15/750 |
|
|
|
[考研] 一志愿西安交通大學(xué) 學(xué)碩 354求調(diào)劑211或者雙一流 +3 | 我想要讀研究生 2026-03-20 | 3/150 |
|
|
|
[考研] 353求調(diào)劑 +3 | 拉鉤不許變 2026-03-20 | 3/150 |
|
|
|
[考研] 085410人工智能專碩317求調(diào)劑（0854都可以） +4 | xbxudjdn 2026-03-18 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +3 | 暗涌afhb 2026-03-16 | 3/150 |
|
|
|
[考研] 生物學(xué)調(diào)劑招人�。�！ +3 | 山海天嵐 2026-03-17 | 4/200 |
|
|
|
[考研] 286求調(diào)劑 +6 | lemonzzn 2026-03-16 | 10/500 |
|
|
|
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +4 | 生物工程調(diào)劑 2026-03-16 | 12/600 |
|
|
|
[考研] 344求調(diào)劑 +6 | knight344 2026-03-16 | 7/350 |
|
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
|
[考研] 材料專碩326求調(diào)劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
|
[考研] 288求調(diào)劑 +4 | 奇點0314 2026-03-14 | 4/200 |
|