function [f] = mci_nmm_r2p6_fx (x,u,P,M)
% Flow for two region, six parameter NMM
% FORMAT [f] = mci_nmm_r2p6_fx (x,u,P,M)
%
% x         State
% u         Inputs
% P         Parameters
% M         Model structure
%
% f         Flow, dx/dt
%__________________________________________________________________________
% Copyright (C) 2015 Wellcome Trust Centre for Neuroimaging

% Will Penny and Biswa Sengupta
% $Id: mci_nmm_r2p6_fx.m 6548 2015-09-11 12:39:47Z will $

curr_P=M.can_P; % Canonical parameter set

% Extrinsic connections
curr_P.A{1}(2,1)=P(1); % Forward connection, w_21
curr_P.A{2}(1,2)=P(2); % Backward connection, w_12

% Intrinsic connections
curr_P.H(1)=P(3);
curr_P.H(2)=P(4);
curr_P.H(3)=P(5);
curr_P.H(4)=P(6);

f = mci_nmm_fx_delay(x,u,curr_P,M);