%generatore per una distribuzione esponenziale
clear;
clf;
a=0;
b=5;
q=2;
n=7000;
fa=exp(-q*a);
fb=exp(-q*b);
cc=fa-fb;
rd=rand(1,n);
xx=-log(fa-cc*rd)/q;
dx=(b-a)/15;
xgr=a:dx:b;
nd=histc(xx,xgr);
nd=nd/(n*dx);
y=exp(-q.*(xgr+dx/2)).*(q/cc);
bar(xgr,nd);
hold on
plot(xgr,y,'rs--','LineWidth',2,'MarkerSize',3);
hold off