function rate = rate_fn(state)



% this is for use in conjunction with numerical integration schemes

% input vector "state" contains [rx vx]'

% output vector "rate" sends back [rxdot, rvdot]


global A;                       % share A matrix with main program



rate = A*state;