| 4 | 1/1 | 返回列表 |
| 查看: 813 | 回復(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 |
木蟲 (正式寫手)

鐵蟲 (初入文壇)
鐵蟲 (初入文壇)
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0703化學(xué)調(diào)劑 ,六級已過,有科研經(jīng)歷 +14 | 曦熙兮 2026-03-15 | 14/700 |
|
|---|---|---|---|---|
|
[考研] 311求調(diào)劑 +3 | 勇敢的小吳 2026-03-20 | 3/150 |
|
|
[考研] 313求調(diào)劑 +4 | 肆叁貳壹22 2026-03-19 | 4/200 |
|
|
[考研] 材料 271求調(diào)劑 +5 | 展信悅_ 2026-03-21 | 5/250 |
|
|
[考研] 298求調(diào)劑 +4 | 上岸6666@ 2026-03-20 | 4/200 |
|
|
[考研] 生物學(xué)一志愿985,分?jǐn)?shù)349求調(diào)劑 +3 | zxts12 2026-03-21 | 3/150 |
|
|
[考研] 一志愿華中科技大學(xué)071000,求調(diào)劑 +3 | 沿岸有貝殼6 2026-03-21 | 3/150 |
|
|
[考研] 306求調(diào)劑 +4 | chuanzhu川燭 2026-03-18 | 4/200 |
|
|
[考研]
|
.6lL 2026-03-18 | 8/400 |
|
|
[考研] 317求調(diào)劑 +5 | 申子申申 2026-03-19 | 9/450 |
|
|
[考研] 一志愿 南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩 +5 | @taotao 2026-03-20 | 5/250 |
|
|
[考博] 招收博士1-2人 +3 | QGZDSYS 2026-03-18 | 3/150 |
|
|
[考研] 085410人工智能專碩317求調(diào)劑(0854都可以) +4 | xbxudjdn 2026-03-18 | 4/200 |
|
|
[考研]
|
不想起名字112 2026-03-19 | 3/150 |
|
|
[考研] 344求調(diào)劑 +6 | knight344 2026-03-16 | 7/350 |
|
|
[考研] 材料專碩326求調(diào)劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
[考博] 26申博 +4 | 八6八68 2026-03-16 | 4/200 |
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
[考研] 機(jī)械專碩325,尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|