function rate = frate_fn(state)

% this is for use in conjunction with numerical integration schemes
% input vector contains"state" X
% output vector "rate" sends back Xdot

global A F;                       % share A matrix with main program

rate = A*state+F;