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

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

%plot trajectory in state space
plot(x(:,1),x(:,2)), title('state space')
hold on