qui entrent dans l’expression de W et en ne conservant que les termes importants, écrire
log(W) = N*(log(N) — Iog(delta(omega)) — Sigma(rho*log(rho))*delta(omega).
Le premier terme est constant au cours de la transformation, le second varie avec la distribution. En remplaçant par une intégrale la somme qui figure dans ce terme et qui est étendue à tous les éléments d’extension en phase, nous obtenons, en représentant le premier terme par une constante A :
log(W) = A — sum(rho*log(rho)*d(omega)).
Si nous admettons qu’à chaque instant soit réalisée la distribution la plus probable, nous pouvons remplacer p par l’expression (21) et il vient, en tenant compte des conditions (22) :
(25) log(W) = A + U/Thêta — N*log(C).
Différentions cette dernière relation
d(log(W)) = dU/Thêta — U/((Thêta)^2)*d(Thêta) — N*(dC/C).
Différentions également la seconde des conditions (22) ; elle donne
N*(dC/C) + (U/(Thêta^2))*d(Thêta) — (C/Thêta)*sum(exp(-E/Thêta))*dE*d(omega) = 0,
