24小时热门版块排行榜    

查看: 816  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 欢迎你来了 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 年年想打听,年年打听不到。。 +11 hdzw9071 2024-06-03 12/600 2024-06-03 14:57 by gy116024
[论文投稿] 投稿的时候编辑发来邮件说作者机构信息没填写怎么办 5+3 飞飞0421 2024-06-02 7/350 2024-06-03 13:39 by topedit
[硕博家园] 实验室太吵闹,无法安静学习,怎么办? +11 utahh 2024-05-31 22/1100 2024-06-03 11:18 by utahh
[论文投稿] 选期刊 5+3 jfdhj 2024-05-29 5/250 2024-06-03 07:47 by 黑大环境队长
[基金申请] 2024杰青和万人领军什么时候会评 +6 墨香琴韵 2024-06-02 6/300 2024-06-03 07:24 by llhljsy
[硕博家园] 又想换工作 +15 brightmj 2024-05-27 26/1300 2024-06-03 07:01 by zyqchem
[论文投稿] 求助大神,Fe和Al离子对MOF都有淬灭,当两种离子共存时,怎么区分两种离子? 10+3 maoxiao 2024-06-02 3/150 2024-06-03 01:43 by nono2009
[基金申请] +5 河马の史诗 2024-06-02 5/250 2024-06-02 19:05 by 328838485
[基金申请] 讨论贴:好本子一定会被打招呼的淘汰吗 +40 ioe3937 2024-05-27 51/2550 2024-06-02 13:52 by Howard28
[基金申请] 九部门发文:不得将专利授权数量作为人才评价、项目评审、职称评定、高校评价等的条件 +17 sjtu2012 2024-05-28 21/1050 2024-06-02 13:43 by 欢乐颂叶蓁
[考博] 导师不让硕转博,让我去国外读博,能理解吗? +12 萧山幽谷 2024-05-29 20/1000 2024-06-02 12:01 by yuan0806
[教师之家] 在大地上我们只过一生---看完我的阿勒泰上头了好几天,完结那天晚上几乎失眠 +11 瞬息宇宙 2024-05-27 13/650 2024-06-01 22:15 by otani
[论文投稿] 求Sci期刊推荐 10+4 甄小鱼 2024-05-30 7/350 2024-06-01 10:41 by bobvan
[基金申请] 博后特别资助状态变化 +24 随梦而飞2017 2024-05-30 35/1750 2024-06-01 10:10 by 青岛阳仔
[考博] 申博求助?本硕双非一篇三区和四区去985工科非天坑专业有没有希望? +4 LYS1200 2024-05-29 6/300 2024-06-01 08:50 by 学术渣渣神
[教师之家] 中年 (金币+3) +18 459582015 2024-05-28 19/950 2024-06-01 00:41 by 沈婉婷.Girl
[硕博家园] 求问endnote上Wiley引用格式在哪下载 +3 Elon. 2024-05-30 3/150 2024-05-31 10:13 by crazy peng
[文学芳草园] 对对子啊 +5 天若孤独 2024-05-29 7/350 2024-05-31 09:00 by wjykycg
[基金申请] 信息学部函评结束了吗? +6 ducan21 2024-05-28 7/350 2024-05-29 12:10 by WORLD0256
[论文投稿] 核心初审被拒,理由是“选题的意义不明确,文章写得不像是科技论文”,怎么改 5+3 工藤雷花樱 2024-05-27 8/400 2024-05-29 10:09 by topedit
信息提示
请填处理意见