mappy_ds.vector_field

mappy_ds.vector_field(dim)

Declare the right-hand side of an ODE \(\dot y = f(y, p)\).

Parameters:

dimint

Dimension of the state space.

Examples

>>> import numpy as np
>>> from mappy_ds import vector_field
>>> @vector_field(dim=2)
... def izhikevich(y, p):
...     v, u = y
...     return np.array([0.04 * v**2 + 5 * v + 140 - u + p["I"],
...                      p["a"] * (p["b"] * v - u)])
>>> izhikevich.jacobian([0.0, 0.0], {"a": 0.2, "b": 0.2, "I": 10.0})
array([[ 5.  , -1.  ],
       [ 0.04, -0.2 ]])