radau5dae_vanderpol.py

assimulo.examples.radau5dae_vanderpol.run_example(with_plots=True)[source]

Example for the use of Radau5DAE 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= 10^6\).

on return:

  • imp_mod problem instance
  • imp_sim solver instance

Final Run Statistics: Van der Pol (implicit) 

 Number of steps                                 : 233
 Number of function evaluations                  : 2144
 Number of Jacobian evaluations                  : 157
 Number of function eval. due to Jacobian eval.  : 314
 Number of error test failures                   : 9
 Number of LU decompositions                     : 235

Solver options:

 Solver                  : Radau5 (implicit)
 Tolerances (absolute)   : 0.0001
 Tolerances (relative)   : 0.0001

Simulation interval    : 0.0 - 2.0 seconds.
Elapsed simulation time: 0.0611169338226 seconds.
_images/radau5dae_vanderpol.png

Note

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