我在半隐式欧拉的Add global damping位置加了一句：
m->velocity = (1 – 0.00005)*m->velocity;
就可以停下来了。因为感觉阻尼的存在就是和速度对抗，让速度降下来（不过老师讲要根据弹簧两个端点的相对速度计算，这么做感觉有点偷懒了。。。真是迷）
显示欧拉，因为一开始速度是NULL，就在Mass的构造函数给velocity设了初值（0，0）（顺手把last_position也给了（0，0）的初值）然后使用这一帧的速度计算当前位置，结果绳子直接起飞了。。。。。真是不知道怎么写才对了。。。。
verlet可以用当前质点位置的前一帧、后一帧的二阶泰勒展开解出来，不过里面带了一个加速度项；那求加速度不是又得在物理上分析受力了吗？感觉没有啥思路实现。。。。
感觉这次实验写的好迷啊，哪都不清楚，求点拨
This post has received 1 vote up.
