clc;
close?all;
clear?all;
hold?off;
r=2;n=20;%r為指數分布參數；
temp=0;?
Un=rand（n1）;%共產生n*1個隨機數
for?i=1:1:n
????S（i）=-log（Un（i）/r）/r;?%Y?=?log（X）：對參量X的每個元素計算自然對數（e為底數的對數函數），使X服從指數分布；
end
S=subs（S）;%S?=?subs（S）：替換表達式S中所有的符號變量，其值通過回調函數或MATLAB工作空間獲取。
for?i=1:1:n
????for?j=1:1:i
????????temp=temp+S（j）;????
????en