<runge-kutta-4> ( dxn..n/dt delta initial-state t-limit -- seq )

Vocabulary

math.runge-kutta

Inputs

dxn..n/dt	an object
delta	an object
initial-state	an object
t-limit	an object

Outputs

seq

an object

Word description

Simple runge-kutta implementation for generating 4th-order approximated solutions for a set of first order differential equations

Examples

A lorenz attractor is a popular system to model with this:

USING: math.runge-kutta math.runge-kutta.examples ; lorenz.

note that the produced chart is a 2 dimensional representation of a 3 dimensional solution. Similarly, the rabinovich-fabrikant system (stable alpha-gamma limit cycle):

USING: math.runge-kutta math.runge-kutta.examples ; rabinovich-fabrikant.

Definition

USING: kernel sequences ;

IN: math.runge-kutta

: <runge-kutta-4>
    ( dxn..n/dt delta initial-state t-limit -- seq )
    time-limit-predicate
    [ [ (runge-kutta-4) ] [ 3drop f ] if ] compose 2with follow
    ; inline