#1.1-1.2
from scipy.special import gamma, factorial
import matplotlib.pyplot as plt
import numpy as np
from scipy.stats import norm



x = np.linspace(-3.5, 5.5, 2251)
y = gamma(x)
plt.plot(x, y, 'b', alpha=0.6, label='gamma(x)')
plt.xlim(-3.5, 5.5)
plt.ylim(-10, 25)
plt.grid()
plt.xlabel('x')
plt.legend(loc='lower right')
plt.show()

X1 = norm(5,2)
x1 = np.linspace(-0, 10, 1000)
y1= X1.pdf(x1)
plt.plot(x1, y1, "k-")
plt.show()

X2 = norm(5,2)
x2 = np.linspace(-0, 10, 1000)
y2= X2.cdf(x2)
plt.plot(x2, y2, "k-")
plt.show()

X3 = norm(5,2)
x3 = np.linspace(-0, 10, 1000)
y3= X3.sf(x3)
plt.plot(x3, y3, "k-")
plt.show()

X4 = norm(5,2)
x4 = np.linspace(-0, 10, 1000)
y4= X4.ppf(x4)
plt.plot(x4, y4, "k-")
plt.show()

X5 = norm(5,2)
x5 = np.linspace(-0, 10, 1000)
y5= X5.isf(x5)
plt.plot(x5, y5, "k-")
plt.show()