w = (1+sqrt(5))/2;

t=0:0.05:100;

x = sin(t);
y = sin(w*t);

clf;
plot(x,y);

title (['Lissajous figures: omega = ' num2str(w)]);
xlabel('X');
ylabel('Y');

PRINT;
