Quelle  cosh.c

/* cosh.c
 *
 * Hyperbolic cosine
 *
 *
 *
 * SYNOPSIS:
 *
 * double x, y, cosh();
 *
 * y = cosh( x );
 *
 *
 *
 * DESCRIPTION:
 *
 * Returns hyperbolic cosine of argument in the range MINLOG to
 * MAXLOG.
 *
 * cosh(x)  =  ( exp(x) + exp(-x) )/2.
 *
 *
 *
 * ACCURACY:
 *
 *                      Relative error:
 * arithmetic   domain     # trials      peak         rms
 *    DEC       +- 88       50000       4.0e-17     7.7e-18
 *    IEEE     +-MAXLOG     30000       2.6e-16     5.7e-17
 *
 *
 * ERROR MESSAGES:
 *
 *   message         condition      value returned
 * cosh overflow    |x| > MAXLOG       MAXNUM
 *
 *
 */

/* cosh.c */

/*
Cephes Math Library Release 2.8:  June, 2000
Copyright 1985, 1995, 2000 by Stephen L. Moshier
*/

#include "mconf.h"
#ifdef ANSIPROT
extern double exp(double);
extern int isnan(double);
extern int isfinite(double);
#else
double exp();
int isnan(), isfinite();
#endif
extern double MAXLOG, INFINITY, LOGE2;

double cosh(x) double x;
{
  double y;

#ifdef NANS
  if (isnan(x))
    return (x);
#endif
  if (x < 0)
    x = -x;
  if (x > (MAXLOG + LOGE2)) {
    mtherr("cosh", OVERFLOW);
    return (INFINITY);
  }
  if (x >= (MAXLOG - LOGE2)) {
    y = exp(0.5 * x);
    y = (0.5 * y) * y;
    return (y);
  }
  y = exp(x);
  y = 0.5 * (y + 1.0 / y);
  return (y);
}