| 查看: 694 | 回复: 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 |
回复此楼» 猜你喜欢
面上本子正文33页，违规吗？会被低分嘛？
已经有7人回复
-
上海大学实验技术岗位非升即走
已经有8人回复
-
今年审到国自然15份，谈谈感受
已经有14人回复
-
考博自荐
已经有6人回复
-
青C资助名额大幅增加！
已经有16人回复
-
重磅！青年科学基金项目（C类）资助增幅预计超过50%
已经有10人回复
-
我在等一个没有答案的答案
已经有3人回复
-
半夜喝咖啡
已经有5人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有7人回复
徐谊
木虫 (著名写手)
- 应助: 25 (小学生)
- 金币: 3693.4
- 散金: 2438
- 红花: 9
- 帖子: 1184
- 在线: 603.9小时
- 虫号: 1214236
- 注册: 2011-02-26
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2014-12-14 22:17:37
3楼2014-12-15 09:31:30
snoopyzhao
至尊木虫 (职业作家)
- 程序强帖: 16
- 应助: 157 (高中生)
- 贵宾: 0.02
- 金币: 18844.7
- 红花: 29
- 帖子: 3803
- 在线: 1422.4小时
- 虫号: 183750
- 注册: 2006-02-13
- 专业: 污染生态化学
4楼2014-12-21 03:55:28
snoopyzhao
至尊木虫 (职业作家)
- 程序强帖: 16
- 应助: 157 (高中生)
- 贵宾: 0.02
- 金币: 18844.7
- 红花: 29
- 帖子: 3803
- 在线: 1422.4小时
- 虫号: 183750
- 注册: 2006-02-13
- 专业: 污染生态化学
5楼2014-12-21 03:57:21
