%%%%%%%%%%%%%%%%%%%%%
%solve differential equation by numerical integration
%and plot the trajectory in state space
%%%%%%%%%%%%%%%%%%%%%

%numerical integration to solve differential equation
[t,x] = ode23('linearf',t0,tfinal,x0,tol,trace);

%plot the states versus time
clf   %clear the plot
plot(t,x), title('time history of both states')
xlabel('time'), ylabel('states')