ZigangPan/DifferentialEquationSolver

ZigangPan`DifferentialEquationSolver`

nflowV

n

flowV

[f,xc,t0,xc0,ti,tf]

calculates the numerical solution to the system of ordinary differential equation

xc'[t]=f[xc[t]]

;

xc[t0]=xc0

; where

f

is the vector field (a pure function),

xc

is the state vector,

xc0

is the initial condition at time

t0

, and

[ti,tf]

is the time interval on which solution is sought. The return of the function is a function sol[t_] on the interval

[ti,tf]

that satisfies the desired differential equation and the initial condition up to numerical precision. Note that

f

is a function of

xc

but not its indivual elements

xc

=(xc1,xc2,...,xcn)

Details and Options



Examples

(1)

Basic Examples

(1)

In[1]:=

f[x_]:={-x〚1〛,x〚1〛^2-x〚2〛};xc={x1,x2};t0=0;xc0={2,1};ti=t0;tf=10;

In[2]:=

sol=

nflowV

[f,xc,t0,xc0,ti,tf];

In[3]:=

Plot[sol[t],{t,ti,tf},PlotRangeAll]

Out[3]=

SeeAlso

▪

▪

▪

▪

""