% Lissajous figures by CYM on Apr. 3, 2007
% plots Lissajous for 100 different phases
t=[0:0.01:10];
% b is phase difference
b=pi/2;

xa = input('Type frequency of x-axis sinusoid (1<= fx <= 20)');
ya = input('Type frequency of y-axis sinusoid (1<= fy <= 20)');

for i = 1:100
    x = cos(xa.*t);
    y = cos( ya.*t + b);
    plot(x,y);
    b=b+0.1;
    
    xx = input('Press any key');

end;