loi=[0.4 0.1 0.5]';nbsimul=10;nmax=100;
//trajectoires du processus normalisé (cas surcritique)
m=(0:length(loi)-1)*loi;
trajectoires=zeros(nbsimul,nmax);
for k=1:nbsimul
z=1;
for n=1:nmax-1
if z>0 then
z=(0:length(loi)-1)*grand(1,'mul',z,loi(1:$-1));
trajectoires(k,n+1)=z/m^n;
else break;
end;
end;
end;
plot2d(0:nmax,[zeros(1,nbsimul);trajectoires'])