loi=[0.3 0.2 0.3 0.2]';nbsimul=10000;nmax=100;
//probabilite d’extinction empirique (cas surcritique)
results=zeros(1,nbsimul);
for k=1:nbsimul
z=1;
while z>0&z<nmax
z=(0:length(loi)-1)*grand(1,'mul',z,loi(1:$-1));
end;
if z==0 then results(k)=1;end;
end;
disp('probabilite empirique d''extinction: '+..
string(sum(results)/nbsimul));