And ode15s can reach the end of the time span I want to calculate. I also write a Python version and use solve_ivp to solve this equation. mkm_electro4.py.txt However, the current density jr decrease ...

Describe your issue. I’m currently trying to solve a stiff ODE problem in Python using the LSODA method from solve_ivp. The governing equation is u' = Lu. Here, u is defined as complex 2M-vector, and ...

The analytical solution of the Initial Value Problem (IVP) and the Boundary Value Problem (BVP) is an essential issue in numerous engineering applications. However, it cannot usually be realized by ...