z=1;
phi=(1+sqrt(5))/2;
w=exp(i*(1/(phi+1))*(2*pi));
for j=1:100
z=[z,0,z(end)*w/0.95];
end
plot(z,'g-o')
axis equal
clear all