%Ein matlab-Programm zum Weihnachtsfest 1998 f\"ur die H\"orer der Vorlesung 
%"Numerische lineare Algebra mit einer Einf\"uhrung in matlab".

c = pi*(1+sqrt(5))/2;  x=0:70000;
z=cumsum(exp(i*(c*x.*x+.4))); plot(z,'g'); axis image; hold on;
zz = z(1:700:70000); plot(zz,'ro','MarkerSize',4,'LineWidth',2);
[M,I] = max(imag(z));
plot(z(I),'y*','MarkerSize',10,'LineWidth',2);
axis('off');
