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