%yp = [ ];
clear y1max y2max;
y0 = [1.2 .3];                       % initial condition
options = odeset('RelTol',.0000001,'AbsTol',.000000001)     
tspan = [0:1:300];                   % length of time interval
[tt,yy] = ode45('whale',tspan,y0,options);           % call the ODE solver
plot(tt,yy(:,1),'-.',tt,yy(:,2));                  % plot solutions vs. time
%  plot(yy(:,1),yy(:,2),':');                           % plot in y1-y2 plane
axis([0 2 0 2]);