| 查看: 1441 | 回复: 6 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
superdirac木虫 (正式写手)
|
[交流]
【求助】想问个velocity-corrected verlet 算法的问题
|
||
|
v(t)={v(t+dt/2)+v(t-dt/2}/2 - dt*{a(t+dt)-a(t-dt)}/12 + O(dt**4) 这个式子是怎么推出来的?? 谢谢! 假如将 v(t+dt/2),v(t-dt/2),a(t+dt),a(t-dt)作泰勒展开,让后合并同类项: 得出的是 v(t)={v(t+dt/2)+v(t-dt/2}/2 - dt*{a(t+dt)-a(t-dt)}/16 + O(dt**4) 两者一直对不上! 找个高手指点一下! 呵呵!谢谢! [ Last edited by lei0736 on 2009-11-25 at 12:07 ] |
回复此楼» 猜你喜欢
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有3人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有3人回复
superdirac
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 4269.5
- 散金: 33
- 红花: 2
- 帖子: 536
- 在线: 287.5小时
- 虫号: 522969
- 注册: 2008-03-12
- 专业: 理论和计算化学
★ ★
xuefei06(金币+2,VIP+0):3Q！ 3-19 14:59
xuefei06(金币+2,VIP+0):3Q！ 3-19 14:59
|
呵呵! 这个不是 beeman http://en.wikipedia.org/wiki/Beeman's_algorithm |
5楼2009-03-19 11:15:21
老虎大王
木虫 (著名写手)
- 模拟EPI: 17
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
2楼2009-03-18 20:29:23
superdirac
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 4269.5
- 散金: 33
- 红花: 2
- 帖子: 536
- 在线: 287.5小时
- 虫号: 522969
- 注册: 2008-03-12
- 专业: 理论和计算化学
3楼2009-03-18 22:04:05
老虎大王
木虫 (著名写手)
- 模拟EPI: 17
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
4楼2009-03-19 10:17:10
