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

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

24小時(shí)熱門版塊排行榜

返回列表

小丫鴨

鐵蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 224
帖子: 16
在線: 6.1小時(shí)
蟲(chóng)號(hào): 2805936
注冊(cè): 2013-11-16
專業(yè): 化工熱力學(xué)和基礎(chǔ)數(shù)據(jù)

[求助] Fortran中的一個(gè)高斯程序已有1人參與

各個(gè)大神門  有沒(méi)有見(jiàn)過(guò)下面這個(gè)程序的可以給解答下這是求什么的嗎？編寫的什么公式？？？
subroutine gauss(n,ns,BB)
implicit none
integer i,j,k,n1,nt,n,ns,ip,i1
double precision BB(2,3),x,AA(3,4)

n1=n+1
nt=n+ns
do i=1,2
do j=1,3
aa(i+1,j+1)=bb(i,j)
enddo
enddo
if (n.gt.1) then
   do 555 i=2,n
      ip=i
      i1=ip
      x=abs(aa(i1,i1))
      do 444 j=i+1,n+1
   if (abs(aa(j,i1)).ge.x) then
         x=abs(aa(j,i1))
         ip=j
      endif
444       continue

if (ip.ne.I1) then
      do 333 j=i1,nt+1
x=aa(i1,j)
aa(i1,j)=aa(ip,j)
         aa(ip,j)=x
333       continue
endif

do 222 j=i+1,n+1
   if (abs(aa(i1,i1)).lt.1.d-10) then
         aa(i1,i1)=1.0
      endif
x=aa(j,i1)/aa(i1,i1)
   do 111 k=i+1,nt+1
aa(j,k)=aa(j,k)-x*aa(i1,k)
111    continue
222 continue
555    continue
endif

   do 50 ip=1,n
   i=n1-ip+1
   do 40 k=n1+1,nt+1
      if (abs(aa(i,i)).lt.1.d-10) then
   aa(i,i)=1.0
      endif
   aa(i,k)=aa(i,k)/aa(i,i)

if (i.ne.2) then
      i1=i-1
      do 30 j=2,i1
aa(j,k)=aa(j,k)-aa(i,k)*aa(j,i)
30    continue
endif
40    continue
50 continue

   do i=1,n
do j=1,n+1
bb(i,j)=aa(i+1,j+1)
enddo
enddo
   end

回復(fù)此樓

» 猜你喜歡

哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有7人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有266人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請(qǐng)問(wèn)還有沒(méi)有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有31人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有68人回復(fù)
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(fù)

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

注冊(cè)化工工程師須知已經(jīng)有23人回復(fù)
使用Multiwfn繪制原子軌道圖形、研究原子殼層結(jié)構(gòu)及相對(duì)論效應(yīng)的影響已經(jīng)有17人回復(fù)

1樓 2015-11-03 15:22:21

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

nagami

木蟲(chóng) (正式寫手)

計(jì)算強(qiáng)帖: 1
應(yīng)助: 64 (初中生)
金幣: 604.9
散金: 281
紅花: 33
帖子: 675
在線: 1131.6小時(shí)
蟲(chóng)號(hào): 1912637
注冊(cè): 2012-07-27
性別: GG
專業(yè): 拓?fù)鋵W(xué)

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

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

從上往下看，有
1）列選主元；
2）判斷主元是否在對(duì)角線上，然后行交換
3）一系列約化
整體的編號(hào)很亂，猜測(cè)是列主元Gauss消去，
但這個(gè)BB（2，3）就看不明白了，就算是增廣矩陣，為什么要指定維數(shù)
覺(jué)得F77可讀性太差。。。F90和C++哪個(gè)都比它強(qiáng)

贊一下(1人)

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

小丫鴨

女靠衣裝；男靠金裝

2樓2015-11-04 15:50:01

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

小丫鴨

鐵蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 224
帖子: 16
在線: 6.1小時(shí)
蟲(chóng)號(hào): 2805936
注冊(cè): 2013-11-16
專業(yè): 化工熱力學(xué)和基礎(chǔ)數(shù)據(jù)

引用回帖:

2樓: Originally posted by nagami at 2015-11-04 15:50:01
從上往下看，有
1）列選主元；
2）判斷主元是否在對(duì)角線上，然后行交換
3）一系列約化
整體的編號(hào)很亂，猜測(cè)是列主元Gauss消去，
但這個(gè)BB（2，3）就看不明白了，就算是增廣矩陣，為什么要指定維數(shù)
覺(jué)得F77可 ...

非常感謝你的回復(fù) 問(wèn)過(guò)老師老師也不確定說(shuō)是求逆還是不太懂

贊一下

回復(fù)此樓

3樓2015-11-04 16:56:51

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

小丫鴨

鐵蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 224
帖子: 16
在線: 6.1小時(shí)
蟲(chóng)號(hào): 2805936
注冊(cè): 2013-11-16
專業(yè): 化工熱力學(xué)和基礎(chǔ)數(shù)據(jù)

送紅花一朵

引用回帖:

送朵小紅花給你哈哈

贊一下

回復(fù)此樓

4樓2015-11-04 16:57:49

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主小丫鴨的主題更新

返回列表

普通表情龍兔虎貓

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

[考研] 280求調(diào)劑 +11	咕嚕曉曉 2026-03-18	12/600	2026-03-21 22:40 by ACS Nano——
[考研] 廣西大學(xué)材料導(dǎo)師推薦 +3	夏夏夏小正 2026-03-17	5/250	2026-03-21 22:20 by 金昊ML
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +5	張vvvv 2026-03-15	7/350	2026-03-21 19:23 by ColorlessPI
[考研] 311求調(diào)劑 +3	勇敢的小吳 2026-03-20	3/150	2026-03-21 17:40 by ColorlessPI
[考研] 266求調(diào)劑 +3	哇呼哼呼哼 2026-03-20	3/150	2026-03-21 16:46 by barlinike
[考研] 070300化學(xué)319求調(diào)劑 +7	錦鯉0909 2026-03-17	7/350	2026-03-21 03:46 by JourneyLucky
[考研] 材料 336 求調(diào)劑 +3	An@. 2026-03-18	4/200	2026-03-21 01:39 by JourneyLucky
[考研] 材料專碩英一數(shù)二306 +7	z1z2z3879 2026-03-18	7/350	2026-03-20 23:48 by JourneyLucky
[考研] 北科281學(xué)碩材料求調(diào)劑 +5	tcxiaoxx 2026-03-20	5/250	2026-03-20 21:35 by laoshidan
[考研] 求調(diào)劑一志愿南京航空航天大學(xué)289分 +3	@taotao 2026-03-19	3/150	2026-03-20 21:34 by JourneyLucky
[考研] 一志愿華中農(nóng)業(yè)071010，總分320求調(diào)劑 +3	困困困困坤坤 2026-03-20	3/150	2026-03-20 20:38 by 學(xué)員8dgXkO
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3	葵梓衛(wèi)隊(duì) 2026-03-20	3/150	2026-03-20 19:28 by zhukairuo
[考研] 廣西大學(xué)家禽遺傳育種課題組2026年碩士招生（接收計(jì)算機(jī)專業(yè)調(diào)劑） +3	123阿標(biāo) 2026-03-17	3/150	2026-03-20 15:58 by 飛行琦
[考研] 085410人工智能專碩317求調(diào)劑（0854都可以） +4	xbxudjdn 2026-03-18	4/200	2026-03-20 09:07 by 不168
[考研] 本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑 +4	我不是一根蔥 2026-03-18	4/200	2026-03-19 09:11 by 浮云166
[考研] 290求調(diào)劑 +3	p asserby. 2026-03-15	4/200	2026-03-17 16:35 by wangkm
[考研] 有沒(méi)有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3	TqlXswl 2026-03-16	7/350	2026-03-17 15:23 by TqlXswl
[考研] 11408 一志愿西電，277分求調(diào)劑 +3	zhouzhen654 2026-03-16	3/150	2026-03-17 07:03 by laoshidan
[考研] 321求調(diào)劑 +5	大米飯！ 2026-03-15	5/250	2026-03-16 16:33 by houyaoxu
[考研] 085601材料工程315分求調(diào)劑 +3	yang_0104 2026-03-15	3/150	2026-03-15 10:58 by peike