| 查看: 1582 | 回复: 5 | ||
欢迎你来了新虫 (小有名气)
|
[求助]
运行结果不随输入变化
|
结果可以运行，但是结果不会随着输入变化，这是为什么呀？ 发自小木虫Android客户端 |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有9人回复
-
时间戳今天，20号变了
已经有5人回复
-
2026国自然函评费到账
已经有9人回复
-
放榜前的不淡定
已经有14人回复
-
建议基金发布提前给出明确的时间点
已经有13人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有20人回复
-
科研孤儿太难了
已经有18人回复
-
看来今天不会放榜了？
已经有11人回复
-
时间戳又变了
已经有19人回复
-
filecode，4个jtjc了
已经有14人回复
raymond137
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1823.8
- 红花: 1
- 帖子: 808
- 在线: 62.2小时
- 虫号: 9463435
- 注册: 2018-07-14
- 性别: GG
- 专业: 电子学与信息系统
2楼2021-05-21 07:06:37
欢迎你来了
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 233.1
- 散金: 15
- 帖子: 131
- 在线: 20.9小时
- 虫号: 21285787
- 注册: 2020-02-28
- 专业: 原子和分子物理
|
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 发自小木虫Android客户端 |
3楼2021-05-21 11:18:56
欢迎你来了
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 233.1
- 散金: 15
- 帖子: 131
- 在线: 20.9小时
- 虫号: 21285787
- 注册: 2020-02-28
- 专业: 原子和分子物理
|
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 发自小木虫Android客户端 |
4楼2021-05-21 11:26:32
raymond137
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1823.8
- 红花: 1
- 帖子: 808
- 在线: 62.2小时
- 虫号: 9463435
- 注册: 2018-07-14
- 性别: GG
- 专业: 电子学与信息系统
5楼2021-05-21 11:53:14
欢迎你来了
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 233.1
- 散金: 15
- 帖子: 131
- 在线: 20.9小时
- 虫号: 21285787
- 注册: 2020-02-28
- 专业: 原子和分子物理
6楼2021-05-21 11:59:59
40