| 查看: 1395 | 回復(fù): 13 | ||
| 【懸賞金幣】回答本帖問題,作者夕陽西下將贈送您 20 個金幣 | ||
夕陽西下金蟲 (著名寫手)
|
[求助]
程序無法運行怎么回事
|
|
|
program L implicit none integer,parameter:: nx=100,nt=100,p=2,te=5,E=2 integer:: i,j real v real:: ne(x,t) real:: n0=5E+2,m=2.1,q=0.6E-19,k=0.514E-9,R=50.0,me=1.2E-3 v=n0*k x=0.0 t=0.0 do i=1,nx x=x+i*0.02 do j=1,nt t=t+j*0.01 ne(0,0)=0.6E-19 ne(x+i*0.02,t+(j+1)*0.01)=-25*te/(me*v)*ne(x+(i+1)*0.02,t+j*0.01)+(1+50*te/(me*v))*ne(x+i*0.02,t+j*0.01)-25*te/(me*v)*ne(x+(i-1)*0.02,t+j*0.01) open (unit=10,file="1.txt" ![]() write (10,*) ne(x,t) end do end do stop end |
木蟲 (著名寫手)

至尊木蟲 (職業(yè)作家)
|
你這個程序中有太多的錯誤,先根據(jù)編譯器的提示進行修改…… 舉個例子: real:: ne(x,t) 中的 x,t 都不是常數(shù),那么,你的 ne 的定義就是錯誤的 ne(0,0)=0.6E-19 默認情況下數(shù)組的下標(biāo)是從 1,1 開始的…… ne(x+i*0.02,t+(j+1)*0.01)= 下標(biāo)不能是實數(shù),只能是整數(shù)…… 等等…… |
金蟲 (著名寫手)
至尊木蟲 (職業(yè)作家)
金蟲 (著名寫手)
|
脗茅路魯脛茫擄茂脦脪驢麓驢麓脮芒賂枚魯脤脨貌隆攏 fin.txt 脢盲脠毛路陸魯脤脧碌脢媒 脫婁賂脙脭玫脙麓脢盲脠毛脛脴攏驢脦脪脮芒脩霉脢盲脠毛脭玫脙麓虜?shù)撀睹斆劽様n驢脳脺脢脟脤謾脢戮麓鉚脦貿(mào)A(n,n) -5 11 -5 -1 0 0 0 0 0 0 -5 11 -5 0 0 0 0 0 0 -5 11 -5 -1 0 0 0 0 0 -5 11 -5 0 0 0 0 0 0 B(n,m) 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 module mat_eq contains subroutine solve(A,B,X,N,M) implicit real*8(a-z) integer::i,k,N,M integer::id_max !脰梅脭陋脣脴鹵錨潞脜 real*8::A(N,N),B(N,M),X(N,N) real*8::Aup(N,N),Bup(N,M) !Ab脦陋脭枚鹿茫戮脴脮貿(mào) [AB] real*8::AB(N,N+M) real*8::vtemp1(N+M),vtemp2(N+M) real*8::vtmp(N),xtmp(N) AB(1:N,1:N)=A AB(:,N+1:N+M)=B ! 脮芒露脦脢脟 脕脨脰梅脭陋脧沒脠樓路簍碌脛潞脣脨脛虜驢路脰 do k=1,N-1 elmax=dabs(Ab(k,k)) id_max=k !脮芒露脦脦陋虜茅脮脪脰梅脭陋脣脴 !脮芒露脦魯脤脨貌碌脛脰梅脪陋脛驢碌脛虜?shù)撁娒嚸幝伱嬞T魯脰碌脳卯麓貿(mào)脭陋脣脴賂酶elmax攏盧露酶脢脟脦陋脕脣脮脪魯枚脳卯麓貿(mào)脭陋脣脴露脭脫婁碌脛鹵錨潞脜 do i=k+1,n if (dabs(Ab(i,k))>elmax) then elmax=Ab(i,k) id_max=i end if end do !脰脕麓脣攏盧脪脩戮顱脥錨魯脡虜茅脮脪脳卯麓貿(mào)脭陋脣脴攏盧虜茅脮脪脥錨魯脡脪脭潞貿(mào)脫毛 碌脷k脨脨陸祿祿祿 !陸祿祿祿脕陸脨脨脭陋脣脴攏盧脝盲脣沒虜?shù)擕u盲 vtemp1=Ab(k, ![]() vtemp2=Ab(id_max, ![]() Ab(k, =vtemp2Ab(id_max, =vtemp1 !脪脭脡脧脪祿麓貿(mào)露脦脢脟脦陋陸祿祿祿脕陸脨脨脭陋脣脴攏盧陸祿祿祿脥錨魯脡脪脭潞貿(mào)錄麓擄麓脮脮脧沒脭陋路簍陸酶脨脨 do i=k+1,27 temp=Ab(i,k)/Ab(k,k) Ab(i, =Ab(i, -temp*Ab(k,![]() end do end do !----------------------------- ! 戮顱鹿媒脡脧脪祿虜陸攏盧Ab脪脩戮顱祿爐脦陋脠莽脧脗脨脦脢陸碌脛戮脴脮貿(mào) ! | * * * * # #| ! [A b]= | 0 * * * # #| ! | 0 0 * * # #| ! | 0 0 0 * # #| ! Aup(:, =AB(1:N,1:N)do i=1,m !碌梅脫脙脫脙脡脧脠媒陸脟路陸魯脤脳茅碌脛祿脴麓酶路陸路簍 vtmp=AB(:,N+i) call uptri(Aup,vtmp,xtmp,n) !擄脩錄脝脣茫陸謾鹿沒賂魯脰碌賂酶X X(:,i)=xtmp end do end subroutine solve subroutine uptri(A,b,x,N) !---------------------------------subroutine comment ! Version : V1.0 ! Coded by : syz ! Date : 2010-4-8 !----------------------------------------------------- ! Purpose : 脡脧脠媒陸脟路陸魯脤脳茅碌脛祿脴麓酶路陸路簍 ! Ax=b !----------------------------------------------------- ! Input parameters : ! 1. A(N,N)脧碌脢媒戮脴脮貿(mào) ! 2. b(N)脫脪脧貌脕驢 ! 3. N路陸魯脤脦盧脢媒 ! Output parameters : ! 1. x 路陸魯脤碌脛賂霉 ! 2. ! Common parameters : ! !---------------------------------------------------- implicit real*8(a-z) integer::i,j,N real*8::A(N,N),b(N),x(N) x(N)=b(N)/A(N,N) !祿脴麓酶虜驢路脰 do i=n-1,1,-1 x(i)=b(i) do j=i+1,N x(i)=x(i)-a(i,j)*x(j) end do x(i)=x(i)/A(i,i) end do end subroutine uptri subroutine driver(N,M) !---------------------------------subroutine comment ! Version : V1.0 ! Coded by : syz ! Date : 2010-4-8 !----------------------------------------------------- ! Purpose : 脟媒露爐魯脤脨貌 !----------------------------------------------------- ! Input parameters : ! 1. N 脙貓脢枚 A(N,N) ! 2. M 脙貓脢枚路陸魯脤 X(N,M)攏盧B(N,M) ! Output parameters : ! 1. ! 2. ! P.S : ! N,M 麓脫脦脛錄鎂脰脨露脕脠隆 ! !---------------------------------------------------- implicit real*8(a-z) integer::N,M integer::i,j real*8::A(N,N),B(N,M),X(N,M) read(11,*) A(N,N) read(11,*) B(N,M) call solve(A,B,X,N,M) write(12,*) X(n,m) end subroutine driver end module mat_eq program main !--------------------------------------program comment ! Version : V1.0 ! Coded by : syz ! Date : 2010-4-8 !----------------------------------------------------- ! Purpose : 脧沒脠樓路簍陸芒戮脴脮貿(mào)路陸魯脤 ! !----------------------------------------------------- ! In put data files : ! 1. fin.txt 脢盲脠毛路陸魯脤脧碌脢媒 ! 2. ! Output data files : ! 1. fout.txt 錄脝脣茫陸謾鹿沒 ! 2. !----------------------------------------------------- ! Post Script : ! 1. 脨貓脪陋脳錄鹵賂脢盲脠毛脢媒戮脻 ! 2. 脫脡脫脷脟媒露爐潞爐脢媒碌梅脫脙路陸路簍潞爐脢媒 !----------------------------------------------------- use mat_eq integer::N,M open(unit=11,file='fin.txt') open(unit=12,file='fout.txt') call driver(N,M) end program main |
新蟲 (初入文壇)
金蟲 (著名寫手)
至尊木蟲 (職業(yè)作家)
金蟲 (著名寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境,總分308求調(diào)劑 +6 | 墨墨漠 2026-03-23 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 考研化學(xué)308分求調(diào)劑 +5 | 你好明天你好 2026-03-23 | 6/300 |
|
|
[考研] 276求調(diào)劑 +3 | YNRYG 2026-03-21 | 4/200 |
|
|
[考研] 323求調(diào)劑 +6 | 洼小桶 2026-03-18 | 6/300 |
|
|
[考研] 308求調(diào)劑 +3 | 墨墨漠 2026-03-21 | 3/150 |
|
|
[考研] 305分求調(diào)劑(食品工程) +4 | Sxy112 2026-03-21 | 6/300 |
|
|
[考研] 384求調(diào)劑 +3 | 子系博 2026-03-22 | 4/200 |
|
|
[基金申請] 山東省面上項目限額評審 +4 | 石瑞0426 2026-03-19 | 4/200 |
|
|
[考研] 材料工程專碩 348分求調(diào)劑 +3 | 冬辭. 2026-03-17 | 5/250 |
|
|
[考研] 一志愿深大,0703化學(xué),總分302,求調(diào)劑 +4 | 七月-七七 2026-03-21 | 4/200 |
|
|
[考研] 材料與化工(0856)304求 B區(qū) 調(diào)劑 +3 | 邱gl 2026-03-21 | 3/150 |
|
|
[考研] 265求調(diào)劑 +12 | 梁梁校校 2026-03-19 | 14/700 |
|
|
[考研] 材料工程(專)一志愿985 初試335求調(diào)劑 +3 | hiloiy 2026-03-17 | 4/200 |
|
|
[考研] 求調(diào)劑一志愿南京航空航天大學(xué)289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
[考研]
|
簡木ChuFront 2026-03-19 | 8/400 |
|
|
[考研] 材料學(xué)求調(diào)劑 +4 | Stella_Yao 2026-03-20 | 4/200 |
|
|
[考研] 廣西大學(xué)家禽遺傳育種課題組2026年碩士招生(接收計算機專業(yè)調(diào)劑) +3 | 123阿標(biāo) 2026-03-17 | 3/150 |
|
|
[考研] 本科鄭州大學(xué)物理學(xué)院,一志愿華科070200學(xué)碩,346求調(diào)劑 +4 | 我不是一根蔥 2026-03-18 | 4/200 |
|
|
[考研] 材料,紡織,生物(0856、0710),化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|