rodasode_vanderpol.py

assimulo.examples.rodasode_vanderpol.run_example(with_plots=True)

Example for the use of RODAS to solve Van der Pol’s equation

\[\begin{split}\dot y_1 &= y_2 \\ \dot y_2 &= \mu ((1.-y_1^2) y_2-y_1)\end{split}\]

with \(\mu=1 10^6\).

on return:

• exp_mod problem instance
• exp_sim solver instance

---

Final Run Statistics: Van der Pol (explicit) 

 Number of steps                                 : 308
 Number of function evaluations                  : 1968
 Number of Jacobian evaluations                  : 308
 Number of function eval. due to Jacobian eval.  : 0
 Number of error test failures                   : 15
 Number of LU decompositions                     : 332

Solver options:

 Solver                  : Rodas 
 Tolerances (absolute)   : 0.0001
 Tolerances (relative)   : 0.0001

Simulation interval    : 0.0 - 2.0 seconds.
Elapsed simulation time: 0.0452361106873 seconds.

Note

Press [source] (to the top right) to view the example.