根據數值獲得概率密度pdf和累積密度分布cdf（MATLAB語言）

y=randn(1,3000); % 生成1-by-3000的標准正態分布隨機數
ymin=min(y);
ymax=max(y);
x=linspace(ymin,ymax,20); %將最大最小區間分成20個等分點(19等分),然后分別計算各個區間的個數
yy=hist(y,x); %計算各個區間的個數
yy=yy/length(y); %計算各個區間樣本數所占總樣本數的比例（概率）
bar(x,yy) %畫出概率密度分布圖

s=0;
for i=2:length(x)
s=[s,trapz(x([1:i]),yy([1:i]))]; % 使用trapz模擬積分
end
figure;
plot(x,s,'*-')