| 5 | 1/1 | 返回列表 |
| 查看: 815 | 回復(fù): 3 | ||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||
小丫鴨鐵蟲 (初入文壇)
|
[求助]
Fortran中的一個高斯程序 已有1人參與
|
|
|
各個大神門 有沒有見過下面這個程序的 可以給解答下這是求什么的嗎 ?編寫的什么公式??? 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é)070300 學(xué)碩336求調(diào)劑 +3 | vv迷 2026-03-21 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +10 | 楓橋ZL 2026-03-18 | 12/600 |
|
|
[考研] 材料學(xué)碩333求調(diào)劑 +3 | 北道巷 2026-03-18 | 3/150 |
|
|
[考研] 336求調(diào)劑 +5 | rmc8866 2026-03-21 | 5/250 |
|
|
[考研] 0805材料320求調(diào)劑 +3 | 深海物語 2026-03-20 | 3/150 |
|
|
[考研] 一志愿山大07化學(xué) 332分 四六級已過 本科山東雙非 求調(diào)劑! +3 | 不想理你 2026-03-16 | 3/150 |
|
|
[考研] 310求調(diào)劑 +3 | baibai1314 2026-03-16 | 3/150 |
|
|
[考研]
|
.6lL 2026-03-18 | 8/400 |
|
|
[考研] 288求調(diào)劑 +16 | 于海海海海 2026-03-19 | 16/800 |
|
|
[考研] 319求調(diào)劑 +3 | 小力氣珂珂 2026-03-20 | 3/150 |
|
|
[考研] 招收調(diào)劑碩士 +4 | lidianxing 2026-03-19 | 12/600 |
|
|
[考研]
|
不想起名字112 2026-03-19 | 3/150 |
|
|
[考研] 328求調(diào)劑,英語六級551,有科研經(jīng)歷 +4 | 生物工程調(diào)劑 2026-03-16 | 12/600 |
|
|
[考研] 材料,紡織,生物(0856、0710),化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
[考研] 293求調(diào)劑 +11 | zjl的號 2026-03-16 | 16/800 |
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6 | 太想進步了0608 2026-03-16 | 6/300 |
|