function xprime=linear(t,x);
% file name linear.m 
% right hand side of the linearized equation as required by ode45

global A B; 
xprime=A*x+B*u(t); 
