let e<sup>x</sup> = u
y = 1/2 (u - 1/u)
2uy = u<sup>2</sup> - 1
u<sup>2</sup> - 2yu - 1 = 0
u = [2y +/- (4y<sup>2</sup> + 4)<sup>1/2</sup>]/2
e<sup>x</sup>=y +/- (y<sup>2</sup> + 1)<sup>1/2</sup>
x=ln[y +/- (y<sup>2</sup> + 1)<sup>1/2</sup>]
=ln[y + sqrt(y<sup>2</sup> + 1)]