phaseR is an R package for the qualitative analysis of one and two dimensional autonomous ODE systems, using phase plane methods. Programs are available to identify and classify equilibrium points, plot the direction field, and plot trajectories for multiple initial conditions. In the one dimensional case, a program is also available to plot the phase portrait. Whilst in the two dimensional case, additionally programs are available to plot nullclines and stable/unstable manifolds of saddle points. Many example systems are provided for the user.
The package contains nine main functions for performing phase plane analysis:
drawManifolds: Draws the stable and unstable manifolds
of a saddle point in a two dimensional autonomous ODE system.
findEquilibrium: Identifies a nearby equilibrium point of an
autonomous ODE system based on a specified starting point.
flowField: Plots the flow or velocity field of a one or two
dimensional autonomous ODE system.
the nullclines of a two dimensional autonomous ODE system.
numericalSolution: Numerically solves a two dimensional
autonomous ODE system in order to plot the two dependent variables against
the independent variable.
a simple means of performing a phase plane analysis by typing only numbers
in to the command line.
phasePortrait: Plots the phase
portrait of a one dimensional autonomous ODE system, for use in classifying
stability: Performs stability, or
perturbation, analysis in order to classify equilibria.
trajectory: Numerically solves a one or two dimensional ODE
system to plot trajectories in the phase plane.
In addition, the package contains over 25 derivative functions for example systems. Links to these can be found in the packages index.
An accompanying guide with further information, examples, and exercises, can be found in the doc directory of the packages installation.
This package makes use of the ode function in the package
Michael J. Grayling
Maintainer: Michael J. Grayling <[email protected]>
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.