Вычислить рекурсивным способом гиперболический косинус

unsigned long long fact (int n)
{
  if (n < 0)
    return 0;
 
  unsigned long long result = 1;
  for (int i = 2; i <= n; i++)
    result *= i;
 
  return result;
}
//---------------------------------------------------------------------------
double func (double x, int n)
{
  if (n < 0)
    return 0;
 
  if (n == 0)
    return 1;
 
  return pow(x, 2 * n) / fact(2 * n) + func(x, n - 1);
}

double S(double x, int n)
{
    if (n==0) return 1.0;
    if (n==1) return 1+0.5*x*x;
    double k=x*x/(2*n*(2*n-1));
    return (1+k)*S(x,n-1)-k*S(x,n-2);
}

double mycosh(double x,double eps)
{
    double s=0,an=1.;
    int n=0;
    while(an>eps)
    {
         s+=an;
         n++;
         an*=x*x/((2.*n-1.)*2.*n); 
    }
    return s;
}

#include <iostream>
#include <cmath>
using namespace std;
double mycosh(double x,double eps=1e-6,int n=0,double an=1.)
{
    if(n==0)
        an=1.;
    if(an<eps)
        return 1.+an;
    an*=x*x/((2.*n+1.)*(2.*n+2.)); 
    return an+mycosh(x,eps,n+1,an);
}
int main()
{
    double x;
    cin>>x;
    cout<<cosh(x)<<' '<<mycosh(x);
    return 0;
}