import numpy as np
import matplotlib.pyplot as plt
from scipy.special import gamma
from math import factorial

x=np.arange(0,6,0.01)
plt.plot(x,gamma(x),label='$\Gamma (x)$')
#můžeme spečítat a nakreslit také (n-1)!
f=np.empty(6)
k=np.arange(0,6,1)
for i in range(6):
    f[i]=factorial(i)
plt.scatter(k+1,f,c='red',label='(n-1)!')

plt.xlabel('x')
plt.ylabel('$\Gamma (x)$')
plt.legend()
plt.show()