[zsy312]

%先编辑函数，a为已知数组，s为待求密度分布
function s=dens(a)
clc
N=10;
len=ceil(max(a)/N);
s=zeros(1,len);
for i=1:len
s(i)=length(find(a>N*(i-1)&a<=N*i));
end
s=s/length(a);
%%%%%%%%%%%%%%%%%%
主程序调用：
y=dens(a);
bar(100*y);%百分率
%区间我设成10了，可以改成别的。