%
%   IVPNEX1A.M  Function evaluation for ODE step in IVPNEX1.M 
%     
%   given differential eqn.   z' = Az + g   where g = B*u
%
%   uses function file ivpnex1b.m to evaluate forcing function at time t
%
%   File prepared by J. R. White, UMass-Lowell (Aug. 2003).
%
      function zp = odefile(t,z,flag,AA,BB)
      u = ivpnex1b(t);   g = BB*u;
      zp = AA*z + g;
%
%   end of function
