In Wolfram Alpha:
plot -EXP(t)/2 LN((EXP(t)-1)/(EXP(t)+1)) from t=0 to t=40
Is this how round off errors appear in solutions when subtraction operand orders of magnitude differ more than 2^24
(32-bit float) ?
Link: Equation Plot
Link: What Every Computer Scientist Should Know About Floating-Point Arithmetic