Méthode de simulation

> lambda:= 3; mu:=4;

densité

> plot(a(t), t=0..4); pl1:= %:

Pr(t <= T)

> Int('a(t)', t=0..T); tmp:= value(%); plot(%, T=0..3);

On a distrib uniforme
et on veut une distrib exponentielle

> 'ra()'$10;

> z= tmp; solve(%, T); subs(z=1-y, %);
g:= unapply(%,y);

> 'g(ra())'$10;

> ga:= eval(g)@ra;

> subs(-1/lambda=-1/mu, eval(g)): gb:= % @ra;