24小时热门版块排行榜    

查看: 798  |  回复: 5

欢迎你来了

新虫 (小有名气)

[求助] 运行结果不随输入变化

结果可以运行,但是结果不会随着输入变化,这是为什么呀?

运行结果不随输入变化


发自小木虫Android客户端
回复此楼
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

raymond137

金虫 (正式写手)

2楼2021-05-21 07:06:37
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

欢迎你来了

新虫 (小有名气)

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
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

欢迎你来了

新虫 (小有名气)

引用回帖:
2楼: Originally posted by raymond137 at 2021-05-21 07:06:37
贴出源代码

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
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

raymond137

金虫 (正式写手)

引用回帖:
4楼: Originally posted by 欢迎你来了 at 2021-05-21 11:26:32
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 ...

我是做C++的。你这个是什么语言代码?Fortran?建议你用debug模式,加几个断点,分步跟踪。

发自小木虫Android客户端
5楼2021-05-21 11:53:14
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

欢迎你来了

新虫 (小有名气)

6楼2021-05-21 11:59:59
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 欢迎你来了 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 博三一直没文章怎么办 +24 133456 2024-05-17 38/1900 2024-05-21 07:26 by 133456
[论文投稿] word转成pdf之后公式里面的字体变了,正文字体没变。 +5 1255037206 2024-05-20 7/350 2024-05-21 00:00 by prepg
[教师之家] 谁要做个人网站,可以交流啊! +12 hjc404 2024-05-14 36/1800 2024-05-20 21:11 by 枯井中抽水
[考博] 双非博士还是985科研助理,然后再读博 +6 lxdatj123 2024-05-18 12/600 2024-05-20 21:06 by chgan2
[基金申请] 青年拔尖出炉 +6 babu2015 2024-05-14 7/350 2024-05-20 20:20 by hema6g13
[硕博家园] 又见奇葩学生,可咋办 +24 appleapple2 2024-05-17 29/1450 2024-05-20 17:59 by appleapple2
[基金申请] 有知道工材e01情况的吗? +13 xiaopang8958 2024-05-17 19/950 2024-05-20 09:52 by 蜕变123~
[教师之家] 优秀毕业论文 指导教师,普通老师有希望吗? +6 河西夜郎 2024-05-17 7/350 2024-05-20 09:36 by 晓目崇
[教师之家] “直接受聘正高专业技术职务”怎么理解 +8 ZHONGWU_U 2024-05-17 10/500 2024-05-19 18:29 by Quakerbird
[教师之家] 执念 +7 459582015 2024-05-16 8/400 2024-05-19 18:22 by Quakerbird
[论文投稿] Angew拒稿重投经历大修返回又送审了 48+3 YOONAZHA 2024-05-18 5/250 2024-05-19 11:28 by nono2009
[考博] 想被211以上高校课题组接收 +8 风起沧澜 2024-05-16 10/500 2024-05-19 06:52 by Guo369
[论文投稿] 一个审稿人已返回,也没邀请审稿人,也没人审稿中,编辑要闹啥? 3+3 林师妹 2024-05-14 12/600 2024-05-17 22:33 by 林师妹
[留学DIY] 求助:医生想润 +7 safrina 2024-05-15 7/350 2024-05-17 18:50 by 卡皮吧啦123
[论文投稿] 投稿成功后又想撤回 +5 otani 2024-05-16 5/250 2024-05-17 16:02 by topedit
[硕博家园] 五氯化铌怎么溶解啊 +3 南南枝枝 2024-05-17 5/250 2024-05-17 11:37 by ad_fish
[硕博家园] 领养小孩 +15 Yyds55 2024-05-14 15/750 2024-05-17 09:53 by 田园1983
[考博] 24应届生寻博导 +3 8139173 2024-05-15 4/200 2024-05-16 21:09 by 不容易456
[基金申请] 粗大,科研之友提示有一篇10多年前的文章被阅读 +9 babu2015 2024-05-16 11/550 2024-05-16 14:53 by TsolmonZ
[教师之家] 宋凯 +10 yexuqing 2024-05-14 10/500 2024-05-16 09:21 by shl2112501
信息提示
请填处理意见