| 查看: 1396 | 回復(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è)作家)
金蟲 (著名寫手)
至尊木蟲 (職業(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ā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿武理材料工程348求調(diào)劑 +6 |  ̄^ ̄゜汗 2026-03-19 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 279分求調(diào)劑 一志愿211 +16 | chaojifeixia 2026-03-19 | 18/900 |
|
|
[考研] 考研化學(xué)308分求調(diào)劑 +7 | 你好明天你好 2026-03-23 | 8/400 |
|
|
[考研] 一志愿上海交大生物與醫(yī)藥專碩324分,求調(diào)劑 +5 | jiajunX 2026-03-22 | 5/250 |
|
|
[考研] 350求調(diào)劑 +6 | weudhdk 2026-03-19 | 6/300 |
|
|
[考研] 307求調(diào)劑 +3 | 余意卿 2026-03-21 | 3/150 |
|
|
[考研] 材料與化工085600,總分304,本科有兩篇sci參與,求調(diào)劑 +4 | 幸運的醬醬 2026-03-22 | 5/250 |
|
|
[考研] 275求調(diào)劑 +6 | shansx 2026-03-22 | 8/400 |
|
|
[考研] 生物學(xué)調(diào)劑 +5 | Surekei 2026-03-21 | 5/250 |
|
|
[考研] 318求調(diào)劑 +4 | plum李子 2026-03-21 | 7/350 |
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-22 | 6/300 |
|
|
[考研] 【考研調(diào)劑】化學(xué)專業(yè) 281分,一志愿四川大學(xué),誠心求調(diào)劑 +11 | 吃吃吃才有意義 2026-03-19 | 11/550 |
|
|
[考研] 296求調(diào)劑 +4 | www_q 2026-03-20 | 4/200 |
|
|
[考研] 330求調(diào)劑0854 +3 | assdll 2026-03-21 | 3/150 |
|
|
[考研] 332求調(diào)劑 +3 | 鳳凰院丁真 2026-03-20 | 3/150 |
|
|
[考研] 求調(diào)劑 +6 | Mqqqqqq 2026-03-19 | 6/300 |
|
|
[考研] 299求調(diào)劑 +6 | △小透明* 2026-03-17 | 6/300 |
|
|
[考研] 初始318分求調(diào)劑(有工作經(jīng)驗) +3 | 1911236844 2026-03-17 | 3/150 |
|
|
[考研] 材料,紡織,生物(0856、0710),化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
[考研] 一志愿南京大學(xué),080500材料科學(xué)與工程,調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|