24小时热门版块排行榜    

查看: 1389  |  回复: 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿211,0860总分286食品工程专业求调剂 +8 林林Winnie 2026-03-05 8/400 2026-03-10 18:44 by 清风月
[考研] 2026考研求调剂-材料类-本科211一志愿985-初试301分 +10 虫友233 2026-03-07 10/500 2026-03-10 17:10 by Demonsssss
[考研] 一志愿东北大学080500,英一数二343分求调剂 +8 Seele_v5 2026-03-04 8/400 2026-03-10 14:51 by ruiyingmiao
[考研] 接受26届调剂生 +21 猪猪猪毛 2026-03-06 21/1050 2026-03-10 13:39 by Jack_333
[考研] 家人们 调剂不迷路 看这里 +8 likeihood 2026-03-09 13/650 2026-03-10 08:09 by likeihood
[考研] 320材料与化工,求调剂 +11 鹤遨予卿 2026-03-04 13/650 2026-03-09 19:39 by 30660438
[考研] 313求调剂 +4 Yyt杨1 2026-03-07 5/250 2026-03-09 13:48 by macy2011
[考研] 一志愿苏大材料工程专硕293求调剂 +5 钢铁大炮 2026-03-04 5/250 2026-03-09 09:49 by EBSD
[考研] 290 材料与化工求调剂 +7 Nebulala 2026-03-08 7/350 2026-03-09 05:21 by houyaoxu
[考博] 26年博士申请 +4 科研狗111 2026-03-07 4/200 2026-03-08 21:56 by 0611517sll
[教师之家] 交大前校长王树国:现在最先进的科技并不在大学实验室,而是在企业研究院 +4 zju2000 2026-03-08 6/300 2026-03-08 19:15 by zju2000
[考研] 085701环境工程专业,初试305,均过国家A区线 +7 卡卡来了@ 2026-03-07 8/400 2026-03-08 06:40 by 刘兵
[考研] 材料化工求调剂 +7 博斯特525 2026-03-06 8/400 2026-03-08 03:51 by lfhuang
[考研] 322分 085600求调剂,有互联网+国金及主持省级大创经历 +4 熊境喆 2026-03-05 4/200 2026-03-07 22:04 by 星空星月
[考研] 求调剂,不管什么专业,我是可塑造的人才一枚,希望遇到知己老师捞捞我 +4 13102137290 2026-03-06 5/250 2026-03-07 21:21 by lissomchan
[考研] 诚求调剂,323分 有专利/科研/实习经历 +6 橙子cyx 2026-03-06 8/400 2026-03-07 19:06 by 新篇章DFSS
[考研] 267化工调剂求助 +7 聪少OZ 2026-03-04 7/350 2026-03-07 17:33 by 无际的草原
[考研] 考研282分求调剂,接受跨专业 +4 刘淄博 2026-03-04 7/350 2026-03-05 22:56 by 刘淄博
[考研] 274求调剂 +9 一个学习者 2026-03-04 9/450 2026-03-05 20:56 by 黑衣馒头人
[考研] 材料085600 303求调剂 +7 1bygone 2026-03-04 7/350 2026-03-05 09:39 by oxidpl
信息提示
请填处理意见