double tanh2(double x)
{
  if ( x >= 0.0 )
    return (1.0-exp(-2*x))/(1.0+exp(-2*x));
  else
    return (exp(2*x)-1.0)/(exp(2*x)+1.0);
}