| 5 | 1/1 | 返回列表 |
| 查看: 1399 | 回復(fù): 5 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
歡迎你來(lái)了新蟲 (小有名氣)
|
[求助]
運(yùn)行結(jié)果不隨輸入變化
|
||
|
結(jié)果可以運(yùn)行,但是結(jié)果不會(huì)隨著輸入變化,這是為什么呀? 發(fā)自小木蟲Android客戶端 |
新蟲 (小有名氣)
|
module contants implicit none real::a1, a2 real::ed = 9.91, edm = 14.5 real::b2 = 2.3 real::kb = 1.380649e-23 real::a_i = 0.08 real::a_j = 4.33d-5 real(kind=8)::h = 1.06e-34 real(kind=8)::w = 2.37e3 integer::b1 = 3390end module contants real function gos(x) use contants implicit none real:: x real(kind=4):: m, m1, m2 real(kind=4):: em1 real::a_1, b_1 integer::j, v integer::tv, t0, trot, tread *,x,tv,t0,trot,t a_1 = a1 - (b1/trot) - (b1/t0) a1 = -3*a_i*kb*t/(2*ed) a2 = -3*a_j*kb*t/(2*ed) if (v >= 0 .and. v < 9) then m1 = 0 m2 = 9 b_1 = b1 m = 1 em1 = (m1 + 0.5)*h*w gos= exp(x*em1 + a_1*m1)*((1 - exp((m2 - m1)*(x*b_1*kb + a_1)))/(1 - exp(x*b_1*kb + a_1))) else if (v >= 9 .and. v < 31) then m1 = 9 m2 = 31 b_1 = 0.75*b1 m = 2 em1 = (m1 + 0.5)*h*w gos= exp(x*em1 + a_1*m1)*((1 - exp((m2 - m1)*(x*b_1*kb + a_1)))/(1 - exp(x*b_1*kb + a_1))) else if (v >= 31 .and. v < 55) then m1 = 31 m2 = 55 b_1 = 0.45*b1 m = 3 em1 = (m1 + 0.5)*h*w gos= exp(x*em1 + a_1*m1)*((1 - exp((m2 - m1)*(x*b_1*kb + a_1)))/(1 - exp(x*b_1*kb + a_1))) end if! fl=q+d+c write(*,*) gos return end function program piok !use contants implicit none real::y,x real,external::gos y=gos(x) ! print*,gos end 發(fā)自小木蟲Android客戶端 |
金蟲 (正式寫手)
新蟲 (小有名氣)
|
module contants implicit none real::a1, a2 real::ed = 9.91, edm = 14.5 real::b2 = 2.3 real::kb = 1.380649e-23 real::a_i = 0.08 real::a_j = 4.33d-5 real(kind=8)::h = 1.06e-34 real(kind=8)::w = 2.37e3 integer::b1 = 3390 end module contants real function gos(x) use contants implicit none real:: x real(kind=4):: m, m1, m2 real(kind=4):: em1 real::a_1, b_1 integer::j, v integer::tv, t0, trot, t read *,x,tv,t0,trot,t a_1 = a1 - (b1/trot) - (b1/t0) a1 = -3*a_i*kb*t/(2*ed) a2 = -3*a_j*kb*t/(2*ed) if (v >= 0 .and. v < 9) then m1 = 0 m2 = 9 b_1 = b1 m = 1 em1 = (m1 + 0.5)*h*w gos= exp(x*em1 + a_1*m1)*((1 - exp((m2 - m1)*(x*b_1*kb + a_1)))/(1 - exp(x*b_1*kb + a_1))) else if (v >= 9 .and. v < 31) then m1 = 9 m2 = 31 b_1 = 0.75*b1 m = 2 em1 = (m1 + 0.5)*h*w gos= exp(x*em1 + a_1*m1)*((1 - exp((m2 - m1)*(x*b_1*kb + a_1)))/(1 - exp(x*b_1*kb + a_1))) else if (v >= 31 .and. v < 55) then m1 = 31 m2 = 55 b_1 = 0.45*b1 m = 3 em1 = (m1 + 0.5)*h*w gos= exp(x*em1 + a_1*m1)*((1 - exp((m2 - m1)*(x*b_1*kb + a_1)))/(1 - exp(x*b_1*kb + a_1))) end if ! fl=q+d+c write(*,*) gos return end function program piok !use contants implicit none real::y,x real,external::gos y=gos(x) ! print*,gos end 發(fā)自小木蟲Android客戶端 |
金蟲 (正式寫手)
|
我是做C++的。你這個(gè)是什么語(yǔ)言代碼?Fortran?建議你用debug模式,加幾個(gè)斷點(diǎn),分步跟蹤。 發(fā)自小木蟲Android客戶端 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 308 085701 四六級(jí)已過求調(diào)劑 +5 | 溫喬喬喬喬 2026-03-12 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 307求調(diào)劑 +3 | 超級(jí)伊昂大王 2026-03-12 | 3/150 |
|
|
[考研] 2026考研求調(diào)劑-材料類-本科211一志愿985-初試301分 +10 | 蟲友233 2026-03-07 | 11/550 |
|
|
[考研] 求調(diào)劑 +5 | 18880831720 2026-03-11 | 5/250 |
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +4 | 邱gl 2026-03-11 | 5/250 |
|
|
[考研] 0856化工原理 +5 | z2839474511 2026-03-10 | 5/250 |
|
|
[考研] 0857環(huán)境調(diào)劑 +5 | 熠熠_11 2026-03-10 | 5/250 |
|
|
[考研] 307求調(diào)劑 +6 | 超級(jí)伊昂大王 2026-03-10 | 6/300 |
|
|
[考研] 306求調(diào)劑 +3 | 唐薏薏 2026-03-09 | 3/150 |
|
|
[考研] 0817學(xué)碩華東區(qū)求調(diào)劑 +3 | 30660438 2026-03-08 | 3/150 |
|
|
[考研] 311求調(diào)劑 +3 | 牛乳糖的卡卡 2026-03-10 | 3/150 |
|
|
[考研] 328,0703考生求調(diào)劑,一志愿為東北師范大學(xué) +3 | 觀素律 2026-03-09 | 4/200 |
|
|
[考研] 0703求調(diào)劑 +10 | c297914 2026-03-06 | 11/550 |
|
|
[考研] 材料調(diào)劑 +4 | xxxcm 2026-03-08 | 7/350 |
|
|
[考研] 286求調(diào)劑 +12 | Faune 2026-03-06 | 14/700 |
|
|
[考研]
|
Sixuan wang 2026-03-06 | 7/350 |
|
|
[考研] 0307化學(xué)求調(diào)劑 +6 | 0ok0k 2026-03-07 | 6/300 |
|
|
[考研] 誠(chéng)求調(diào)劑,323分 有專利/科研/實(shí)習(xí)經(jīng)歷 +6 | 橙子cyx 2026-03-06 | 8/400 |
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +4 | 邱gl 2026-03-06 | 4/200 |
|
|
[考研] 08工科求調(diào)劑 +3 | 隆LLL 2026-03-06 | 4/200 |
|