Computes dx/dt for the bilinear neural state equation coupled to the Buxton-Friston hemodynamic model. Mirrors SPM25's spm_fx_fmri.

dcm_fx_fmri(x, u, P, M)

Arguments

x

State matrix (rows = regions, cols = state variables).

u

Driving inputs at the current time.

P

List of model parameters (A, B, C, D, transit, decay, epsilon, ...).

M

Model structure list (optional).

Value

Matrix dx/dt of the same shape as x.