我正在使用boost odeint来进化一个微分方程。时间和值是由包装类围绕一个普通数(进行高斯误差传播)给出的。using Time = Number<double>;状态及其导数由带有声明的自定义类给出。但是,您可能会注意到,我没有使用我在第二行中定义的受控步骤,但是我使用了默认的步骤。当我试图在最后一行中使用stepper变量而不是默认的Stepper()时,我
我正在尝试用scipy.integrate.odeint从最终时间的边界条件求解一个常微分方程系统,并一直工作到初始时间(如这里所描述的:Backward integration in time usingscipy odeint)。然而,odeint迭代为负的时间值-超出了我正在寻找的实际解决方案的范围-这会导致错误,因为我的ODE依赖于时间的平方根,并且我的函数返回复数值而不是实数。为什