def compute_Cp(T,Cv,V,B0,beta):
    """
    This function computes the isobaric heat capacity from the equation:
    
    :math:`Cp - Cv = T V beta^2 B0`

    where *Cp,Cv* are the isobaric and isocoric heat capacities respectively,
    *T* is the temperature, *V* the unit cell volume, *beta* the volumetric
    thermal expansion and *B0* the isothermal bulk modulus.
    
    """
    Cp = Cv + T * V * beta * beta * B0
    return Cp