list for(iter=verlet_vars.verlet()[i].begin();iter!=verlet_vars.verlet()[i].end();iter++) //call the force() function to calculate奇怪的是,如果我在内部循环被调用之前对迭代器做了任何事情,就像这样:
iter=verlet
我完成了之前为Verlet方法所做的代数,没有力--这导致了如下所示的代码,但是当我忽略外力时,"+(2*F/D)“项就没有了。我的问题是,是否有人能发现我在Verlet方法中添加的一个力项的缺陷。% verlet.m% oscillator.% stepsize h, for a second-order ODE
function vout = verlet<