function ode1

close all
clear

T0 = 0;
T1 = 20;
N = 500;
t = linspace(T0,T1,N);
x0 = [1,2];
%[t,x] = ode45(@my_model,t,x0);
[t,x] = ode45(@my_model,[T0 T1],x0);

subplot(2,1,1)
plot(t,x(:,1),'r');
xlabel 't'
ylabel 'x'
title 'ODE Respons'

subplot(2,1,2)
plot(x(:,1),x(:,2),'k')
xlabel 'x'
ylabel 'xdot'
title 'Phase Plane'



