Entropie de mélange pour un gaz ou liquide idéal

#! /usr/bin/env python
# -*- coding: utf-8 -*-
"""
représentation graphique de l'entropie de mélange d'un système idéal
"""
 
import numpy as np  # voir http://docs.scipy.org/doc/
import matplotlib.pyplot as plt
 
def s(t):
    return t*np.log(t) + (1-t) * np.log(1-t)
 
x1 = np.arange(0.0, 1., 0.001)
 
plt.plot(x1, s(x1), 'b-')
#plt.plot(x1, x1*np.log(x1) + (1-x1) * np.log(1-x1), 'b-')   #autre façon, n'utilisant pas la fonction
 
plt.show()
 
# des messages du type "RuntimeWarning: divide by zero encountered ..."
# ou "RuntimeWarning: invalid value encountered" sont liés aux valeurs
# extrêmes de x (limites 0. et 1.)