f = 70;
for i = 0:511
    x(i+1) = i;
    X(i+1) = 2 * f*pi * i/512;
    Y(i+1) = sin(X(i+1));
end

resultFFT = fft(Y);
plot(abs(resultFFT));