版塊導(dǎo)航: 正在加載中...

登錄注冊

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

24小時熱門版塊排行榜

返回列表

當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖

caoweiyue

新蟲 (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 127
散金: 30
帖子: 89
在線: 15.5小時
蟲號: 1403468
注冊: 2011-09-15
專業(yè): 水力學(xué)與水信息學(xué)

[求助] 剛開始學(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)有5人回復(fù)
081700 調(diào)劑 267分已經(jīng)有8人回復(fù)
請教下大家 2026年國家基金申請是雙盲審嗎？已經(jīng)有5人回復(fù)
276求調(diào)劑。有半年電池和半年高分子實習(xí)經(jīng)歷已經(jīng)有10人回復(fù)
一志愿南航材料專317分求調(diào)劑已經(jīng)有4人回復(fù)
求調(diào)劑已經(jīng)有4人回復(fù)
生物學(xué)學(xué)碩求調(diào)劑已經(jīng)有5人回復(fù)
284求調(diào)劑已經(jīng)有10人回復(fù)
一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑已經(jīng)有4人回復(fù)
07化學(xué)280分求調(diào)劑已經(jīng)有4人回復(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ù)

1樓 2012-05-16 18:39:43

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

jwwang_lq

鐵蟲 (初入文壇)

應(yīng)助: 2 (幼兒園)
金幣: 17
帖子: 7
在線: 1.6小時
蟲號: 1719996
注冊: 2012-03-27
性別: GG
專業(yè): 生物大分子結(jié)構(gòu)與功能

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

write語句里應(yīng)該輸出b(i)，而不是x(i)

贊一下

回復(fù)此樓

老實人

3樓2012-05-17 08:51:20

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

查看全部 4 個回答

Jasonluo

至尊木蟲 (知名作家)

中央書籍處書記

應(yīng)助: 61 (初中生)
貴賓: 0.01
金幣: 22628.1
散金: 166
紅花: 14
帖子: 9258
在線: 680.9小時
蟲號: 518969
注冊: 2008-03-05
性別: GG
專業(yè): 自然地理學(xué)

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

樓主中間那個害羞的表情是啥意思啊？呵呵

贊一下

回復(fù)此樓

朝為田舍郎，暮登天子堂。將相本無種，學(xué)子當(dāng)自強！

2樓2012-05-16 20:59:37

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

daizhijuan

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 12.5
帖子: 1
在線: 27分鐘
蟲號: 2184351
注冊: 2012-12-13
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

【答案】應(yīng)助回帖

write(*,"('x',i1,'=',f6.4)" i,x(i)
少了右邊的括號
real:: A(n,n),b(n),x(n)
子程序中的這句之前最好也給n一個值：比如 parameter（n=4）

贊一下

回復(fù)此樓

4樓2012-12-13 16:52:02

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

查看全部 4 個回答

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

[基金申請] 請教下大家 2026年國家基金申請是雙盲審嗎？ +3	lishucheng1 2026-03-22	5/250	2026-03-24 08:22 by gltch
[考研] 求調(diào)劑 +7	十三加油 2026-03-21	7/350	2026-03-23 23:48 by 熱情沙漠
[考研] 0854電子信息求調(diào)劑 324 +3	Promise-jyl 2026-03-23	3/150	2026-03-23 13:43 by wangkm
[考研] 333求調(diào)劑 +6	87639 2026-03-21	10/500	2026-03-23 10:41 by Iveryant
[考研] 求調(diào)劑材料學(xué)碩080500，總分289分 5+3	@taotao 2026-03-19	21/1050	2026-03-23 10:17 by 冠c哥
[考研] 293求調(diào)劑 +3	濤濤Wjt 2026-03-22	5/250	2026-03-22 22:21 by jiangpengfei
[考研] 0854電子信息求調(diào)劑 +3	α____ 2026-03-22	3/150	2026-03-22 21:28 by zhq0425
[考研] 354求調(diào)劑 +7	Tyoumou 2026-03-18	10/500	2026-03-22 11:11 by 人來盛
[考研] 資源與環(huán)境調(diào)劑申請(333分) +5	holy J 2026-03-21	5/250	2026-03-21 22:42 by Catalysis25
[考研] 求調(diào)劑 +3	13341 2026-03-20	3/150	2026-03-21 18:28 by 學(xué)員8dgXkO
[考研] 297求調(diào)劑 +11	戲精丹丹丹 2026-03-17	12/600	2026-03-21 17:47 by ColorlessPI
[考研] 336求調(diào)劑 +5	rmc8866 2026-03-21	5/250	2026-03-21 17:24 by 學(xué)員8dgXkO
[考研] 332求調(diào)劑 +3	鳳凰院丁真 2026-03-20	3/150	2026-03-21 10:27 by luoyongfeng
[考研] 化學(xué)求調(diào)劑 +4	臨澤境llllll 2026-03-17	5/250	2026-03-21 02:23 by JourneyLucky
[考研] 一志愿南昌大學(xué)，327分，材料與化工085600 +9	Ncdx123456 2026-03-19	9/450	2026-03-20 23:41 by lovewei0727
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑，總分315（英一） +5	sbdksD 2026-03-19	5/250	2026-03-20 22:10 by luoyongfeng
[考研] 329求調(diào)劑 +9	想上學(xué)吖吖 2026-03-19	9/450	2026-03-20 22:01 by luoyongfeng
[考研] 求調(diào)劑一志愿南京航空航天大學(xué)289分 +3	@taotao 2026-03-19	3/150	2026-03-20 21:34 by JourneyLucky
[考研] 295復(fù)試調(diào)劑 +8	簡木ChuFront 2026-03-19	8/400	2026-03-20 20:44 by zhukairuo
[考研] 生物學(xué)調(diào)劑招人�。�！ +3	山海天嵐 2026-03-17	4/200	2026-03-19 21:34 by 怎么釋懷

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時熱門版塊排行榜

caoweiyue

[求助] 剛開始學(xué)Fortran 幫忙看一下這個程序 為什么算不出結(jié)果

» 猜你喜歡

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

jwwang_lq

【答案】應(yīng)助回帖

Jasonluo

【答案】應(yīng)助回帖

daizhijuan

【答案】應(yīng)助回帖

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

[求助] 剛開始學(xué)Fortran 幫忙看一下這個程序為什么算不出結(jié)果

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