% ODE solver

clear
clf

y0=[1 -2];
tspan=[0 100];
[T,Y]=ode45(@g,tspan,y0);

% plots
figure(1)
plot(Y(1,1),Y(1,2),'*')
hold on
plot(Y(:,1),Y(:,2))

figure(2)
plot(T,Y(:,1),'g')
hold on
plot(T,Y(:,2),'r')


