Power Spectrum using fft

I am given the time series data which are sampled at 1000 Hz frequency. The length of the data is also 1000. By using fft, I must obtain the power spectrum of the original signal as a function of frequency. How can I do it?

[Bish Erbas]

Use fft

https://www.mathworks.com/help/matlab/ref/fft.html