24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 642  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 依旧是则徐 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 332求调剂 +12 032500 2026-03-25 12/600 2026-03-30 12:16 by 学员8dgXkO
[有机交流] 考研调剂 +5 watb 2026-03-26 5/250 2026-03-30 11:51 by yjolah
[考研] 070300一志愿211,312分求调剂院校 +4 小黄鸭宝 2026-03-30 4/200 2026-03-30 11:49 by 棒棒球手
[考研] 297求调剂 +13 田洪有 2026-03-26 13/650 2026-03-30 11:33 by 替代品000
[考研] 求化学调剂 +11 wulanna 2026-03-28 11/550 2026-03-30 10:59 by 探123
[考研] 285求调剂 +5 AZMK 2026-03-29 8/400 2026-03-30 10:41 by AZMK
[考研] 375求调剂 +6 雨夏整夜 2026-03-29 6/300 2026-03-30 10:21 by herarysara
[考研] 322求调剂:一志愿湖南大学 材料与化工(085600),已过六级。 +4 XX小邓 2026-03-29 4/200 2026-03-29 17:34 by 无际的草原
[考研] 329求调剂 +10 钮恩雪 2026-03-25 10/500 2026-03-29 13:32 by peike
[考研] 279求调剂 +4 蝶舞轻绕 2026-03-29 4/200 2026-03-29 09:45 by laoshidan
[考研] 调剂求院校招收 +6 鹤鲸鸽 2026-03-28 6/300 2026-03-29 08:15 by fmesaito
[考研] 调剂考研 +3 王杰一 2026-03-29 3/150 2026-03-29 08:09 by fmesaito
[考研] 材料求调剂 一志愿哈工大总分298分,前三科223分 +5 dongfang59 2026-03-27 5/250 2026-03-28 04:53 by wxiongid
[考研] 315分求调剂 +7 26考研上岸版26 2026-03-26 7/350 2026-03-28 04:05 by fmesaito
[考研] 265求调剂11408 +3 刘小鹿lu 2026-03-27 3/150 2026-03-27 20:53 by nihaoar
[考研] 298调剂 +3 jiyingjie123 2026-03-27 3/150 2026-03-27 11:57 by wxiongid
[考研] 341求调剂 +7 青柠檬1 2026-03-26 7/350 2026-03-27 00:19 by wxiongid
[考研] 材料科学与工程 317求调剂 +4 JKSOIID 2026-03-26 4/200 2026-03-26 15:58 by 不吃魚的貓
[考研] 0854人工智能方向招收调剂 +4 章小鱼567 2026-03-24 4/200 2026-03-25 13:29 by 2177681040
[考研] 材料考研调剂生 +3 黄粱一梦千年 2026-03-24 3/150 2026-03-24 17:00 by barlinike
信息提示
请填处理意见