| 6 | 1/1 | 返回列表 |
| 查看: 1398 | 回復: 5 | ||
歡迎你來了新蟲 (小有名氣)
|
[求助]
運行結(jié)果不隨輸入變化
|
|
結(jié)果可以運行,但是結(jié)果不會隨著輸入變化,這是為什么呀? 發(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客戶端 |
金蟲 (正式寫手)
新蟲 (小有名氣)
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[基金申請]
NSFC申報書里申請人簡歷中代表性論著還需要在申報書最后的附件里面再上傳一遍嗎
20+3
|
NSFC2026我來了 2026-03-10 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 材料與化工085600調(diào)劑求老師收留 +5 | jiaanl 2026-03-11 | 5/250 |
|
|
[考研] 271求調(diào)劑 +4 | 生如夏花… 2026-03-11 | 4/200 |
|
|
[考研] 341求調(diào)劑 +3 | 搗蛋豬豬 2026-03-11 | 3/150 |
|
|
[考研] 大連大學化學專業(yè)研究生調(diào)劑 +3 | 琪久. 2026-03-10 | 8/400 |
|
|
[考研] 考研材料與化工,求調(diào)劑 +7 | 戲精丹丹丹 2026-03-09 | 7/350 |
|
|
[考研] 哈工大材料324求調(diào)劑 +6 | 閆旭東 2026-03-10 | 8/400 |
|
|
[考研] 0856材料與化工353分求調(diào)劑 +11 | NIFFFfff 2026-03-09 | 11/550 |
|
|
[考研] 求調(diào)劑材料專碩293 +6 | 段_(:з」∠)_ 2026-03-10 | 6/300 |
|
|
[考研] 2026考研求調(diào)劑-材料類-本科211一志愿985-初試301分 +10 | 蟲友233 2026-03-07 | 10/500 |
|
|
[考研] 328,0703考生求調(diào)劑,一志愿為東北師范大學 +3 | 觀素律 2026-03-09 | 4/200 |
|
|
[考研] 求調(diào)劑! +3 | 朔朔話 2026-03-09 | 3/150 |
|
|
[考研] 304求調(diào)劑(085602一志愿985) +8 | 化工人999 2026-03-09 | 8/400 |
|
|
[考研] 復試調(diào)劑 +6 | 呼呼?~+123456 2026-03-08 | 8/400 |
|
|
[考博] 找博導 +4 | 小呆呆熊 2026-03-07 | 4/200 |
|
|
[考研] 一志愿山東大學,總分327,英語二79,有論文,有競賽,已過四六級 +3 | 木木目目1 2026-03-09 | 3/150 |
|
|
[考研] 334求調(diào)劑 +5 | Trying] 2026-03-07 | 5/250 |
|
|
[考研] 一志愿中國石油大學(華東) 本科齊魯工業(yè)大學 求調(diào)劑 +3 | snw石 2026-03-07 | 3/150 |
|
|
[考研] 第一志愿上海大學,專業(yè)化學工程與技術(shù),總分288,求調(diào)劑 +3 | 1829197082 2026-03-07 | 3/150 |
|
|
[考研]
|
zbcm_zbcm 2026-03-05 | 6/300 |
|