为什么我从matplotlib获得两次传奇？

Question

这是我的数字：

这是我的代码：

#!/usr/bin/env python 
import numpy as np
import matplotlib.pyplot as plt
from scipy.integrate import quad 

mu = 0 
def fermions_integrand(x):
    return np.sqrt(x)*((np.log(1+np.exp(-x+mu))) + (1/(np.exp(x-mu)+1))*(x-mu))
def s_fermions(mu):
    return quad(fermions_integrand, 0, np.inf)

def bosons_integrand(x):
    return np.sqrt(x)*((-np.log(1-np.exp(-x+mu))) + (1/(np.exp(x-mu)-1))*(x-mu))
def s_bosons(mu):
    return quad(bosons_integrand, 0, np.inf)


plt.figure(figsize=(8, 5))
plt.title('Fermions Vs Bosons')
plt.xlabel(r'$\mu$')
plt.ylabel('s/g')
plt.plot(np.linspace(-10,10), [s_fermions(mu) for mu in np.linspace(-10,10)], 'b-', lw=3, color='red', label='Fermions')
plt.plot(np.linspace(-10,10), [s_bosons(mu) for mu in np.linspace(-10,10)], 'b-', lw=3, color='blue', label='Bosons')
plt.legend()
plt.grid(True)
plt.show() 

我认为代码很简单，但为什么我会得到这两个传说呢？如何解决这个问题？

进一步的信息：实际上，我已经使用过这个绘图代码一次就另一个问题没有问题（在Ubuntu上）但是对于这个我只是转移到OSX我不确定它是否与这个问题有关不是吗？

Answer 1

这是因为scipy.integrate.quad返回估计的积分值和绝对误差的估计值，作为元组(y,abserr)

然后是两个理解列表

[s_fermions(mu) for mu in np.linspace(-10,10)]
[s_bosons(mu) for mu in np.linspace(-10,10)]

是值列表的列表，您可以绘制积分和绝对误差（实际上在整个域中非常接近于零）。由于您一次绘制两行，因此您有两个图例条目。

只需在[0]和s_bosons的返回行中添加尾随s_fermions即可绘制错误。