24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 641  |  回复: 4

依旧是则徐

铜虫 (小有名气)

[求助] 公共参量设定(子程序里的i,j 怎样设置才能使用上面主程序传入的i,j 值呢) 已有1人参与

各位,我这样写的程序,但是子程序里的i,j 怎样设置才能使用上面主程序传入的i,j 值呢。
!求解Y''(T)=-(i+j*cos(2*T))*Y
program main
   use IMSL
   implicit none
   integer, parameter :: MXPARM = 50
   integer, parameter :: N = 2
   integer IDO, ISTEP, NOUT
   real PARAM(MXPARM), T, TEND, Y(N)
   real, parameter :: TOL = 1E-8
   external FCN
   real :: i,j

  do i=0,1,0.001
      do j=0,1.5,0.1
           write(*,"(I2,I2)" i,j

       PARAM = 0
        WRITE (*,"(' Time', 9X,'Y')"
        Y(1) = 0 ! Initial condition
        Y(2) = 2.0
        T = 0.0    ! Y(T)=4.0中的时间T值
        IDO = 1
          do ISTEP=0,100,10
            TEND = ISTEP
                call IVPRK (IDO, N, FCN, T, TEND, TOL, PARAM, Y)
            write (*,'(F5.1,F12.4)') T, Y
          end do

       call IVPRK (3, N, FCN, T, TEND, TOL, PARAM, Y)
          end do
   end do
     
  stop
   end program
! 计算微分方程式
! YPRIME=Y'=f(t,y)
subroutine FCN (N, T, Y, YPRIME)
   implicit none
   integer N
   real T, Y(N), YPRIME(N)
   real i,j
  !Y''(T)=-(i+j*cos(2*T))*Y
   YPRIME(2) =-(i+j*cos(2*T))*Y(1)
   return
end subroutine
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

徐谊

木虫 (著名写手)

【答案】应助回帖


感谢参与,应助指数 +1
ben_ladeng: 金币+1, 谢谢指教 2014-12-17 10:57:33
i j定义为全局变量就行了啊

[ 发自手机版 http://muchong.com/3g ]
2楼2014-12-14 22:17:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

依旧是则徐

铜虫 (小有名气)

引用回帖:
2楼: Originally posted by 徐谊 at 2014-12-14 22:17:37
i j定义为全局变量就行了啊

就是需要设成全局变量 我在主程序和子程序的real i,j 后面都加上common i,j
对吗。请指教!
3楼2014-12-15 09:31:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

Fortran 中可以把 subroutine 当一个参数传给另外一个 subroutine 吗?实在想不起来有这样的写法了,只记得 function 可以作为参数……
4楼2014-12-21 03:55:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

另外,i, j 这类的变量名,一般还是用来表示 integer 变量吧,否则如果不开 implicit none 的话,容易出错……
5楼2014-12-21 03:57:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 依旧是则徐 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 334分 一志愿武理 材料求调剂 +7 李李不服输 2026-03-26 7/350 2026-03-29 16:25 by 学员8dgXkO
[考研] 299求调剂 +9 15188958825 2026-03-25 9/450 2026-03-29 14:16 by wxiongid
[考研] 297求调剂 +11 田洪有 2026-03-26 11/550 2026-03-29 13:14 by peike
[考研] 349求调剂 +6 李木子啊哈哈 2026-03-25 6/300 2026-03-29 12:47 by 无际的草原
[考研] 求调剂,一志愿 南京航空航天大学 ,080500材料科学与工程学硕,总分289分 +7 @taotao 2026-03-29 7/350 2026-03-29 12:03 by longlotian
[考研] 一志愿华东师范大学有机化学专业,初试351分,复试被刷求调剂! +3 真名有冰 2026-03-29 4/200 2026-03-29 08:47 by qingfeng258
[考研] 321求调剂 +7 璞玉~~ 2026-03-25 8/400 2026-03-29 06:41 by 544594351
[考研] 一志愿哈尔滨工业大学材料与化工方向336分 +9 辰沐5211314 2026-03-26 9/450 2026-03-29 01:12 by 我是小康
[考研] 0703 化学 求调剂,一志愿山东大学 342 分 +4 Shern—- 2026-03-28 4/200 2026-03-29 00:47 by 544594351
[考研] 312,生物学求调剂 +3 小译同学abc 2026-03-28 3/150 2026-03-28 15:32 by 落睿可思
[考研] 0703本科郑州大学求调剂 +3 nhj_ 2026-03-25 3/150 2026-03-28 13:24 by Iveryant
[考研] 272求调剂 +7 脚滑的守法公民 2026-03-27 7/350 2026-03-27 17:23 by laoshidan
[考研] 266分求材料化工冶金矿业等专业的调剂 +4 哇呼哼呼哼 2026-03-26 4/200 2026-03-27 17:02 by zhyzzh
[考研] 307求调剂 +8 超级伊昂大王 2026-03-24 9/450 2026-03-27 15:34 by 超级伊昂大王
[考研] 085601 材料工程 313分 求调剂 +5 Ong3 2026-03-27 5/250 2026-03-27 12:24 by goldfish51
[考研] 271求调剂 +6 生如夏花… 2026-03-22 6/300 2026-03-26 16:48 by 张凯十八号
[考研] 材料专硕 335 分求调剂 +4 拒绝冷暴力 2026-03-25 4/200 2026-03-25 18:45 by haxia
[考研] 347求调剂 +4 L when 2026-03-25 4/200 2026-03-25 13:37 by cocolv
[考研] 生物学学硕求调剂 +7 小羊睡着了? 2026-03-23 10/500 2026-03-25 02:24 by 清风拂扬。 m
[考研] 材料专硕331求调剂 +4 鲜当牛 2026-03-24 4/200 2026-03-24 15:58 by JourneyLucky
信息提示
请填处理意见