function?possion（nr）
%r為指數分布的參數n為產生隨機數個數?
n=100;r=0.529;
R=exprnd（rn1）;
b=zeros（n+11）;?
b（1）=0;
for?i=2:n+1;
????b（i）=b（i-1）+R（i-1）;
end
t=0:0.01:b（n+1）;?
p=length（t）;
nt=zeros（p1）;
j=1;
for?i=1:p;
????i