24小时热门版块排行榜    

查看: 635  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 依旧是则徐 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 301求调剂 +4 A_JiXing 2026-03-16 4/200 2026-03-17 17:32 by ruiyingmiao
[考研] 085600材料与化工 +4 安全上岸! 2026-03-16 4/200 2026-03-17 14:02 by 勇敢太监王公公
[论文投稿] 有没有大佬发小论文能带我个二作 +3 增锐漏人 2026-03-17 4/200 2026-03-17 09:26 by xs74101122
[考研] 药学383 求调剂 +3 药学chy 2026-03-15 4/200 2026-03-16 20:51 by 元子^0^
[基金申请] 今年的国基金是打分制吗? 50+3 zhanghaozhu 2026-03-14 3/150 2026-03-16 17:07 by 北京莱茵润色
[考研] 304求调剂 +4 ahbd 2026-03-14 4/200 2026-03-16 16:48 by 我的船我的海
[考研] 304求调剂 +3 曼殊2266 2026-03-14 3/150 2026-03-16 16:39 by houyaoxu
[考研] 283求调剂 +10 小楼。 2026-03-12 14/700 2026-03-16 16:08 by 13811244083
[考研] 材料与化工一志愿南昌大学327求调剂推荐 +7 Ncdx123456 2026-03-13 8/400 2026-03-16 12:15 by karry wen
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[考研] 【0703化学调剂】-一志愿华中师范大学-六级475 +5 Becho359 2026-03-11 5/250 2026-03-14 11:35 by 哦哦123
[考研] 材料080500调剂求收留 +3 一颗meteor 2026-03-13 3/150 2026-03-14 10:54 by peike
[考研] 285 求调剂 资源与环境 一志愿北京化工大学 +3 未名考生 2026-03-10 3/150 2026-03-13 23:04 by JourneyLucky
[考研] 材料与化工求调剂一志愿 985 总分 295 +8 dream…… 2026-03-12 8/400 2026-03-13 22:17 by 星空星月
[考研] 26调剂/材料/英一数二/总分289/已过A区线 +6 步川酷紫123 2026-03-13 6/300 2026-03-13 21:59 by 星空星月
[考研] 工科,求调剂 +3 我887 2026-03-11 3/150 2026-03-13 21:39 by JourneyLucky
[考研] 求调剂 +5 一定有学上- 2026-03-12 5/250 2026-03-13 18:31 by ms629
[考研] 求b区学校调剂 +3 周56 2026-03-11 3/150 2026-03-13 16:20 by JourneyLucky
[考研] 333求调剂 +3 152697 2026-03-12 4/200 2026-03-13 07:08 by Iveryant
[考博] 2026年博士申请 +3 QwQwQW10 2026-03-11 3/150 2026-03-12 17:58 by gxch43
信息提示
请填处理意见