SODE Interface

type

errorType=(constAbsolute,percent);

variables=array[1..16] of Double;

realFunc=Function(x:variables): Double;

funcArray=array[1..16] of realFunc;

T_ODE=Class(TObject)

numEqu:Integer;

dydx:funcArray;

RKVdx:Double;

tolerance:double;

Constructor Create(nEqu:Integer;func:FuncArray);

Procedure EulerStep(dx:Double;var x:variables);

Procedure RK2Step(dx:Double;var x:variables);

Procedure RK4Step(dx:Double;var x:variables);

Procedure RK4VStep(var x:variables);

Procedure RK4VReport;

end;