Jacobian:
-------------------------------------------------------------------------------
This function computes the Jacobian matrix for a right hand-side.
The central difference approximation good to O(h^2) is used.
For any function f(x,t) it will compute the f0 and a matrices.
f(x,t) ‰ f(xOP,tOP) + a(xOP,tOP) x + Š
funFcn is input as a character string 'xxxxx' and must be of the form
xdot = xxxxx ( x, t, {p1,p2,p3,p4,p5,p6,p7,p8,p9,p10})
which is the same form used in RK2 and RK4. tOP is optional. If not needed
pass [].
-------------------------------------------------------------------------------
Form:
[a, fOP] = Jacobian( funFcn, xOP, tOP, p1, p2, p3, p4, p5, p6, p7, p8, p9, p10 )
-------------------------------------------------------------------------------
------
Inputs
------
funFcn 'function name'
xOP (n,1) State at the operating point
tOP Time
varargin Optional arguments
-------
Outputs
-------
a (n,n) Jacobian matrix of first partials
fOP (n,1) f(xOP,tOP) at the operating point
-------------------------------------------------------------------------------
Children: