clear;
close;

a=0.1;
%a=0.14;
b=0.5;
options=[];

[t y]=ode23('cyclefunc',[0 50],[0.6 1.4],options,a,b);
%plot(y(:,1),y(:,2));