X?=?random（‘exp‘1/2510001）;?%構造服從指數分布的Xn序列，個數為1000個
S?=?[?];
S（1）?=?0;?
for?i?=?1:length（X）??%length（?）測出向量長度??
????S（i+1）?=?S（i）+X（i）;
end?%這個函數的作用是生成顧客進入商場的時間序列
N?=?[?];M?=?[];W?=?[];
N（1）=0;?M（1）=0;W（1）=0;?%規定好初始值為0
b?=?rand（1length（S）-1）;
for?i?=?2:length（S）
????N（i）=i-1;??
????if?b（i-1）?<=?0.4
????????M（i）?=?M（i-1）+1;
????????W（i）?=?W（i-1）;
????else?
????????M（i）=M（i-1）;
????????W（i）=W（i-1）+1;