min et err sont des constantes
un échelon, ben c'est un échelon, toutes les valeurs sont nulles, sauf entre 0 et 1 ou c'est constant
En fait le problème vient d'ailleurs.
Je fais ces deux boucles :
do i=1,n
x(i)=-log(uniform())
tmp(i)=min(i)+x(i)*err(i)
S=S+tmp(i)
enddo
do i=1,n
RB(i)=tmp(i)/S
enddo
où RB et tmp sont des vecteurs de dimension n.
Le but de ces boucles est d'obtenir les RB, tels que leur somme fasse 1.
Uniform() est un tirage aléatoire uniforme entre 0 et 1.
Le problème est que la répartition des RB(i) est censée être uniforme également, alors que quand je trace le graphe j'ai plutôt une répartition normale, voire log-normale.