x = linspace(0, 150, 1000); 
y =  2.0 + (x / 50.0) + sin(x) + 2.0 * sin(pi * x / 50.0);
figure(); clf;
plot(x, y, 'b.');

