function xdot = linodepar(t,x,flag,a,b,u)
%
% b.w. bequette - 4 Sept 13
% pass through the state space a and b matrices
% pass through the input vector
  xdot = a*x + b*u;