# Author: Ika, 2013-07-26
from pylab import *
import pylab as pl
import numpy as np

# Create a figure of size 8x6 points, 80 dots per inch
pl.figure(figsize=(8,8), dpi=80)

# Create a plot of the cubic function y=-1/(x+4)+4
x = np.linspace(-14.0,-4.0,1024, endpoint=True)
y = -1/(x+4)+4
pl.plot(x,y,color="green", linewidth=3.0, linestyle="-")
x = np.linspace(-4.0,6.0,1024, endpoint=True)
y = -1/(x+4)+4
pl.plot(x,y,color="green", linewidth=3.0, linestyle="-")
# Draw the limit line
x = np.linspace(-14.0,6.0,1024, endpoint=True)
y = 4+0*x
pl.plot(x,y,color="red", linewidth=1.0, linestyle="--")

# Set labels
pl.xlim(-14.0,6.0)
pl.ylim(-6.0,14.0)

# Move the spines
ax = pl.gca()
ax.set_xlabel('x')
ax.set_ylabel('y')
ax.text(-13, 13, r'$f(x)=\frac{-1}{x+4}+4$', fontsize=32)
ax.text(-13, 10, r'$\lim_{x\rightarrow +\infty}f(x)=4$', fontsize=32)

# Set up the spines
ax.spines['right'].set_color('none')
ax.spines['top'].set_color('none')
ax.xaxis.set_ticks_position('bottom')
ax.spines['bottom'].set_position(('data',0))
ax.yaxis.set_ticks_position('left')
ax.spines['left'].set_position(('data',0))

# Save the figure to the output SVG file
plt.savefig("Tamasol_SVG.svg");