I was trying to compute an ODE of attitude kinematics eqs in Python using the solve_ivp function, but the problem is that one of the parameters, the a ...

Here is my code which integrates two coupled ODEs In Jupyter, this code works properly, but if I uncomment the second @njit, which is just above de ...

I have code written that discretizes a PDE and solves it using finite difference method I am trying to extend this to a 2D case and am trying to us ...

For training purpose, I try to modeling and analyze a dynamic system using python. I follow the example given on this site. In this example, they are ...

. Answers to this question are eligible for a +50 reputation bounty. Ra ...

I am supposed to solve a ODE as follows: where H_*p, H_m, data2_*sleep are objects that are results of interpolation using scipy.interp1d . For th ...

I wanted to make a bifurcation plot to demonstrate the bistability region for the following system of equations: Toggle switch equations I have alrea ...

I want to pass an array input signal to ordinary differential equations solver 'odeint' instead of defining a function. Much like the linear solver ls ...

I have implemented a simple simulation on a uniform circular motion using the low level API scipy.integrate.RK45() like the following. I tried at f ...

I am trying to solve this equation using Runge Kutta 4th order: applying d2Q/dt2=F(y,x,v) and dQ/dt=u Q=y in my program. I try to run the code but ...

I am trying to use MATLAB to solve the ODE involving step function Here is my code to solve the above ODE: However, the result seems strange and ...

I want to create an ODE file for the equation d^2theta/dt^2 - cos(theta) = 0 I have this: But I don't think this is right as dYdt (1) would be for ...

I was trying to reproduce some results of ode45 solver in Python using solve_ivp. Though all parameters, initial conditions, step size, and 'atol' and ...

For the time, t, from 0 to 30, z0 is a constant value of 6. For the time,t, from 30 to 100, z0 takes on the form of a time variable where z0 = 6exp(-0 ...

For the 2 systems of ODE, I am using RK4 to solve. From 0 <= t <= 30, b is a constant. But at t >= 30, b is a time-dependent variable where b ...

I am trying to replicate a lake food web model published here. The model represents two food chains (littoral vs pelagic), connected by a top predator ...

I'm trying to solve an ODE problem (with Julia), that can stop early when a specific condition is satisfied. Let's say I have a Lorenz system as belo ...

I have been trying to plot a graph that looks like the one given below for a model that I saw on a paper. The R codes for the model are as follows. ...

I am working with an SIR model using R. I need to plot multiple curves of I on the same figure for different values of beta, say for the values of 0.0 ...

I've got this equation from mathematical model to know the thermal behavior of a battery. dTsdt = Ts * a+ Ta * b + dTadt * c + d However, ...