solve -t = 2/a*(sqrt(h) -b/a*ln(1+a*sqrt(h)/b)) for h