z=0.001
p=0
a=exp(i*4*pi/(3+sqrt(5))+2*pi*i/100*p)
for j=1:500
z=[z,0.995*exp(-j*0.000000001)*z(1,end)*a];
end
plot(z,'o')
clear all