global a
global b

a=0.02;
b=3;
hold
z=[0 1 1]
for i=1:1000
   z=duffmap(z(2:3));
   plot(z(2),z(3))
end
xlabel('x')
ylabel('y')