I am trying to solve the ODE of the frictional free fall: I think the solution might be wrong. I believe it should be the reciprocal of the returned expression. Therefore, if I try to find the value of the constant for an initial condition of rest (v(0) = 0), I get no real solutions: As a consequence, if I
Tag: ode
Numerical Python: Solving a BVP with a boolean condition?
I am not sure about the best way to ask this question, but I am trying to find the long-term state of an ODE system with an arbitrary extra constraint that needs to be fulfilled. Ex: In the code above, f is the right-hand-side of my ODE model. N is the dimension of my state vector X, t_range is the
How to put initial condition of ODE at a specific time point using odeint in Python?
How to put initial condition of ODE at a specific time point using odeint in Python? So I have y(0) = 5 as initial condition, following code works:: I wanna see the graph in both negative and positive time line. So I change t = np.linspace(0,20) to t = np.linspace(-5,20), but then the initial condition is taken as y(-5) =
How can I stop my Runge-Kutta2 (Heun) method from exploding?
I am currently trying to write some python code to solve an arbitrary system of first order ODEs, using a general explicit Runge-Kutta method defined by the values alpha, gamma (both vectors of dimension m) and beta (lower triangular matrix of dimension m x m) of the Butcher table which are passed in by the user. My code appears to