o = linspace(0, 10* pi, 500);
for i = 1:500;
    r = 0.8*(0.8^o);
    x(i) = r*cos(o);
    y(i) = r*sin(o);
end;
plot(x,y)