| 查看: 945 | 回复: 6 | |||
[求助]
二阶常微分方程的编写思路 已有2人参与
|
|
求高手解决一下方程的编写思路 最好详细一些。或者给出相似的计算例子也行 y··（t）+2y·（t)+y=0 初始条件y（0）=4 y·（0）=2 |
回复此楼» 猜你喜欢
268求调剂
已经有9人回复
-
301求调剂
已经有4人回复
-
0703化学336分求调剂
已经有5人回复
-
301求调剂
已经有4人回复
-
材料专硕306英一数二
已经有9人回复
-
304求调剂
已经有8人回复
-
材料工程专硕调剂
已经有5人回复
-
344求调剂
已经有4人回复
-
326求调剂
已经有6人回复
-
296求调剂
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于二阶偏微分方程求解
已经有21人回复
- 师兄师姐们，帮小弟看看这个 Laplace变换后的二阶常微分方程怎么解啊
已经有23人回复
- 有限差分法 第二类边界条件
已经有12人回复
- 非线性二阶微分方程组求解，matlab
已经有12人回复
- 求高手帮忙编写 常微分方程组中 高阶导数用其他低阶导数表示的 简单数学软件代码
已经有8人回复
- 微分方程求解（通解中含hypergeom函数）
已经有17人回复
- 中心差分求解二维热传导的matlab程序
已经有12人回复
- Matlab数值求解非线性常微分方程
已经有5人回复
- 非线性二阶常微分方程组+第三类边界条件
已经有17人回复
- 一个二阶常微分方程的求解方法？
已经有5人回复
- 如何用 数值解法 求 二阶微分方程的本征值
已经有18人回复
- 微分方程组求参数问题，求高人指教，感谢
已经有11人回复
- Matlab用龙格库塔法求二阶非线性方程
已经有13人回复
- 二阶非线性微分方程求解
已经有14人回复
- 如何求解描述振动的二阶微分方程
已经有21人回复
- 请问如何求解二元一阶微分方程组
已经有24人回复
- matlab-常微分方程参数估计
已经有12人回复
- 求助三元函数二阶偏微分C语言写法
已经有10人回复
- 【求助】求助一个一阶常微分方程的初值问题
已经有16人回复
- 【求助】请教二阶微分方程的特解（含初始条件）
已经有8人回复
- 【求助】有限元法求解二阶偏微分方程组
已经有6人回复
- 【求助】一个常微分方程的通解
已经有6人回复
angocn
至尊木虫 (文坛精英)
- 应助: 206 (大学生)
- 金币: 17572.5
- 散金: 100
- 红花: 33
- 帖子: 11967
- 在线: 451.8小时
- 虫号: 2595191
- 注册: 2013-08-13
- 性别: GG
- 专业: 光学
2楼2014-12-09 13:23:20
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68127.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
3楼2014-12-09 21:56:02
4楼2014-12-10 10:09:25
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68127.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
5楼2014-12-10 13:52:56
|
版主，按照你的帮助，我这样写的程序，但是子程序里的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,jPARAM = 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 |
6楼2014-12-14 11:45:56
| 祝福。 |
7楼2014-12-14 11:47:24
