# Author: Fredrik.
# Enhanced by Ika, 2014-02-08

import numpy as np
import math
import matplotlib.pyplot as plt
import scipy.special as sp

plt.figure()
plt.axes([0.17,0.13,0.79,0.8])
plt.hold(True)

Q = []

# Draw the function curve.
X = np.arange(0.001, 5, 0.001)
Y = []
for x in X:
	Y.append(math.gamma(x+1))
a = plt.plot(X, Y, '-', color='red', lw=3)
Q.append(a)
a = plt.plot([1,2,3,4], [1,2,6,24], 'ro')
Q.append(a)

plt.xlim(0,5)
plt.ylim(0,25)
plt.xlabel('n',fontsize=24)
plt.ylabel('n!',fontsize=24)
# Set up labels
plt.text(2,10,'?',fontsize=72)
plt.text(1,20,r'$f(n)=\Gamma(n+1)$',fontsize=36)
# Set up the spines
plt.grid()

plt.savefig("Factorial_Interpolation.svg")