from numpy import linspace, append
from math import sinh, exp
from matplotlib.pyplot import *
from mpl_toolkits.axes_grid.axislines import SubplotZero

fig = figure(figsize=(5,7))
ax = SubplotZero(fig,111)
fig.add_subplot(ax)
ax.grid(True)
ax.set_ylim((-13,15))
for direction in ["xzero","yzero"]:
	ax.axis[direction].set_axisline_style("-|>")
	ax.axis[direction].set_visible(True)
for direction in ["left","right","bottom","top"]:
	ax.axis[direction].set_visible(False)

t = linspace(-3,3,50)
H0,H1,H2 = [],[],[]
for i in t:
	H1 = append(H1,exp(i))
	H2 = append(H2,exp(-i))
#        H0 = append(H0,sinh(i))  # either this
        H0 = append(H0,0.5*(exp(i)-exp(-i)))  # or this
ax.plot(t,H0,label=r"$\mathrm{sinh}(x)$")
ax.plot(t,H1,label=r"$e^x$")
ax.plot(t,H2,label=r"$e^{-x}$")

t = linspace(-2.5,2.5,11)
for i in t:
	H0 = sinh(i)
	H1 = exp(i)
	H2 = exp(-i)
	ax.plot([i,i,i],[H0,H1,H2],'yo-.')
ax.text(3,0.5,r"x")
ax.text(-0.5,14.5,r"y")
ax.legend(frameon=False)
ax.minorticks_on()
#fig.show()
fig.savefig("Hyperbolic_and_exponential;_sinh.png",bbox_inches="tight",\
	pad_inches=.15)