Plot a Histogram of the Amplitudes of various frequencies of a Signal

Hello, I need to plot a bar Histogram of a signal after performing a fft. It is an ENF data, so most amplitudes should be around 60Hz region and around within +/-0.5 Hz. I want an output like the figure I attached. The signal 'x' has length L=10000

code:

    x = x(1:1000,1);
    T = 1/fs;
    L=length(x);
    t = (0:L-1)*T;

    Y = fft(x);
    P2 = abs(Y/L);
    P1 = P2(1:L/2+1);
    P1(2:end-1) = 2*P1(2:end-1);

    f = fs*(0:(L/2))/L;
    plot(f,P1)
    title('Single-Sided Amplitude Spectrum of X(t)')
    xlabel('f (Hz)')

I tried using the code

figure(2) 
hist(P1,f);

But it doesn't work as I expect. Please help out. thanks.

[rodrigo figueiredo]

X = randi ([59.5,60.5], 100,1);

[Walter Roberson]

X = randi ([59.5,60.5], 100,1);

Error using randi
First input must be a positive scalar integer value IMAX, or two integer values [IMIN IMAX] with IMIN less than or equal to IMAX.

Nope, not that.

If you want a uniform distribution 59.5 to 60.5 then

X = 59.5 + rand(100,1);

[Yazan]

Can you provide the data?

[Adam Danz]

See histfit.