ode: Solve arbitrary first order ODE systems

In shill1729/odeSolveR: Implementation of simple Euler schemes and fourth-order Runge-Kutta methods for arbitrary first-order ODE systems

Description

A solver for arbitrary first-order ODE systems of finite dimension, using either RK4 (default) or an Euler scheme.

Usage

ode(dynamics, IC, t0 = 0, tn = 1, n = 1000, engine = "RK4")

Arguments

dynamics	a list of functions with as many variables as IC, and in the same order with the names
IC	a list of variables with the same names as the non-time input as f, in the same order with initial values
t0	initial time
tn	ending time
n	number of sub-intervals in time-grid
engine	which scheme to use to solve the ODE; 'RK4' or 'Euler', see details.

Details

The list of functions must have syntax e.g. f(t,x,y,z) for each element matching the named elements of IC as list(x = y0, y = y0, z = z0) for example. Further, the engine must be either "RK4" for a fourth-order Runge-Kutta method, or "Euler" for an Euler scheme.

Value

data.frame

shill1729/odeSolveR documentation built on March 31, 2021, 10:52 a.m.