n=10000;
X=repmat(0,n,1);
i=0;
while i<=n
u1=unifrnd(0,1);
u2=exprnd(2);
x=u2;
if u1<u2*exp(-u2/2)/(2*exp(-1))
i=i+1;
X(i)=u2;
end
end
hist(X)
Envién sus sugerencias
En la parte inferior de la página se ha incorporado un formulario en el cual ustedes pueden enviar sus sugerencias, sea de contenido, sea de forma, etc., todo esto para mejorar la página.
No hay comentarios:
Publicar un comentario