clear all
format long


% specify the time domain with starting time t_o
% and final time t_f

   t_o=
   t_f=
% specify the initial positon of the trajectory

    x_o=
    y_o=

% call the integration routine
        
[tt,y] = ode45(@dble,[t_o t_f],[x_o y_o]);

% plot the trajectory just calculated in the x-y domain

line(y(:,1),y(:,2))
% modify the axis to the domain 0<x<2 and 0<y<1
axis([0 2 0 1]);