import matplotlib.pyplot as plt import numpy as np from scipy.special import factorial x = np.array(np.arange(0, 10, 0.2)) y_factorial = factorial(x) y_n = x**x plt.plot(x, y_factorial, color='r') plt.plot(x, y_n, color='g') plt.show()