import numpy as np
import matplotlib.pyplot as plt
from  scipy.special import gamma
from scipy.stats import t

N=1000
eps=0.05
x=np.linspace(-5,5,N)
y=np.zeros([N,10])
legend=np.empty(10,dtype=object)

for i in range(1,11):
    y[:,i-1]=t.pdf(x,i)
    plt.plot(x,y)
    legend[i-1]='nu = '+str(i)
plt.xlabel('x')
plt.ylabel('pdf(x)')
plt.legend(legend)
plt.show()