在我的常微分方程函数中,我需要迭代地求解一个参数的方程,直到在每个时间步收敛。我希望传递最新的参数值作为下一个时间步的初始值,这样当函数对参数进行迭代更新时,花费的时间会更少。但我想不出该怎么做。nn<=100: #update a until convergence在某种语言中,我可以返回要由集成器使用的更新后的参数
我正在使用刚性求解器(ode15s)对常微分方程系统进行时间积分。它是有效的,但我想加快速度。它使用x参数来求解Poisson型二维方程(通过有限体积法)。力F与泊松方程的解成正比。)/aP;
F = sum(Phi(:)); % discrete integration over domain通过迭代方法求解泊松方程</em