double tanh1(double x)
{
  return (exp(x)-exp(-x))/(exp(x)+exp(-x));
}