% spectrum   // Compute fft of y and plot as Fig 2
figure(2)
plot(abs(fft(y)))
