Plot3D[{Exp[x]*(Cosh[y] + Sqrt[Sinh[y]^2 + Exp[-4*x]]) , Exp[x]*(Cosh[y] - Sqrt[Sinh[y]^2 + Exp[-4*x]])},{x,0,3},{y,-2,2}]