24小时热门版块排行榜    

查看: 1390  |  回复: 5
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

欢迎你来了

新虫 (小有名气)

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

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

运行结果不随输入变化


发自小木虫Android客户端
回复此楼
已阅   关注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 个回答

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 320材料与化工,求调剂 +12 鹤遨予卿 2026-03-04 14/700 2026-03-10 22:37 by peike
[考研] 一志愿中科大材料与化工。353分能调剂到哪? +4 MayUxw1 2026-03-08 6/300 2026-03-10 17:09 by MayUxw1
[考研] 一志愿天津大学,英一数二305分求调剂,四六级已过 +7 小小番的茄 2026-03-09 7/350 2026-03-10 16:48 by ztnimte
[考研] 328,0703考生求调剂,一志愿为东北师范大学 +3 观素律 2026-03-09 4/200 2026-03-10 15:58 by houyaoxu
[考研] 304求调剂(085602一志愿985) +8 化工人999 2026-03-09 8/400 2026-03-10 15:21 by houyaoxu
[考研] 085701环境工程295求调剂,一志愿A区211 +8 等春来, 2026-03-06 9/450 2026-03-10 10:38 by 沙漠之狐994
[考研] 家人们 调剂不迷路 看这里 +8 likeihood 2026-03-09 13/650 2026-03-10 08:09 by likeihood
[考研] 307求调剂 +3 辛仁豆腐 2026-03-08 5/250 2026-03-09 14:09 by macy2011
[考研] 【求调剂】293分环境工程求调剂材料/化工,服从调剂,抗压能力强! +13 xiiiia 2026-03-04 14/700 2026-03-09 14:06 by macy2011
[考研] 中科大材料299求调剂 +10 DAIjiayo 2026-03-05 16/800 2026-03-09 10:45 by 斩魂滴兔子!
[考研] 347求调剂 +4 浮云满足 2026-03-07 4/200 2026-03-08 16:46 by 星空星月
[考研] 求调剂 +4 呼呼?~+123456 2026-03-06 4/200 2026-03-06 23:11 by L135790
[考研] 一志愿中石油,化学工程,301求调剂 +4 豚馒枳枳 2026-03-04 4/200 2026-03-06 22:15 by 星空星月
[考研] 08工科求调剂 +3 隆LLL 2026-03-06 4/200 2026-03-06 15:26 by 隆LLL
[考研] 材料调剂 +5 ws 上岸鸭 2026-03-05 5/250 2026-03-05 20:06 by houyaoxu
[考研] 304求调剂 +4 曼殊2266 2026-03-05 4/200 2026-03-05 17:10 by zhukairuo
[考研] 一志愿西交化工专硕288专业课93求调剂 +6 好运好运接接 2026-03-04 6/300 2026-03-05 09:36 by houyaoxu
[考研] 286 +6 ksncj 2026-03-04 6/300 2026-03-04 20:49 by lature00
[考研] 0703化学求调剂,总分320分,一志愿中南民族大学 +3 19286712440 2026-03-04 3/150 2026-03-04 18:30 by houyaoxu
[考研] 一志愿西交大材料学硕 346 求调剂 +3 zju51 2026-03-04 3/150 2026-03-04 13:01 by zhukairuo
信息提示
请填处理意见