clc
t=[0:0.000001:0.1];
x=sin(20*pi*t);
k=linspace(-1,1,100);
M=hist(x,100);
M=M/length(t);
plot(k,M)