function [y]=Exp(m,n,lambda) // lambda exp(-lambda x) x>=0 // --------------------------- y=(-1/lambda)* log(rand(m,n,'uniform'));
Simulation code:
n=1000; lambda=3; y=Exp(1,n,lambda); histplot([0:0.1:10],y,[1,1],'051',' ',[0,0,5,3]); deff('[y]=f(x)','y=lambda*exp(-lambda*x);'); x=[0:0.1:10]';plot2d(x,f(x),1,"000"); titre= 'macro histplot : Histogram plot'; xtitle(titre,'Classes','N(C)/Nmax');