% circle map in matlan, Eli Tziperman
clear; clf
pi=atan(1.0)*4.0;
omega=0.35;
K=1.75
I=50;

theta(1:I)=0.3001;

for n=1:I
  theta(n+1)=mod(theta(n)+omega+(K/(2*pi))*sin(2*pi*theta(n)),1);
end

plot(theta,'r+')
%hold on