RK4:

In faustogc/streamChaos: Analysis of Non Linear Data Streams

Usage

RK4(t, syst, y, dt, parms)

Arguments

t
syst
y
dt
parms

Examples

##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function (t, syst, y, dt, parms) 
{
    k1 <- syst(t, y, parms)[[1]] * dt
    k2 <- syst(t + dt/2, y + dt/2 * k1, parms)[[1]] * dt
    k3 <- syst(t + dt/2, y + dt/2 * k2, parms)[[1]] * dt
    k4 <- syst(t + dt, y + dt * k3, parms)[[1]] * dt
    dy <- y + dt/6 * (k1 + 2 * k2 + 2 * k3 + k4)
    list(c(dy))
  }

faustogc/streamChaos documentation built on May 8, 2019, 9:22 a.m.