from scipy.special import gamma as Gamma
from pylab import *

def f1(x):
    return Gamma(x)


x = linspace(-10, 10, 512)
y1 = f1(x)
gca().set_autoscale_on(False)

# Matlab-style syntax:
plot(x, y1)

xlabel('x')
ylabel('y')
# legend(r'$\Gamma(x)$')
axis([-10, 10, -100, 1000])
grid(True)

show()