Model function
modelFun(Y, A, K)
Y
time vector
A
the matrix A that is found in: dY/dt = A*Y + F
K
the zero matrix for the seven states
No return value, called for side effects