:: [a] package:numeric-prelude

Example for a linear equation: Setup a differential equation for y with
y   t = (exp (-t)) * (sin t)
y'  t = -(exp (-t)) * (sin t) + (exp (-t)) * (cos t)
y'' t = -2 * (exp (-t)) * (cos t)
Thus the differential equation
y'' = -2 * (y' + y)
holds. The following function generates a power series for exp (-t) * sin t by solving the differential equation.
equalTrunc 500 PSE.expExpl PSE.expODE
equalTrunc 500 PSE.sinExpl PSE.sinODE
equalTrunc 500 PSE.cosExpl PSE.cosODE
equalTrunc 500 PSE.logExpl PSE.logODE
equalTrunc 500 PSE.atanExpl PSE.atanODE
equalTrunc 500 PSE.sinhExpl PSE.sinhODE
equalTrunc 500 PSE.coshExpl PSE.coshODE