Hi Giles. Thank you for your first-rate article. Looking at equations 11 and 12 though, it looks as if a couple of typographical errors may have crept in. I reckon that equation 11 should be 32*sech^2(0.182*t) [where "sech^2" represents 'sech squared'] and that equation 12 should be 967*ln(cosh(0.182*t)) [where "ln" represents 'log to the base e'] Kind regards David Williams