Вычисление гиперболического косинуса через разложение в ряд Тейлора с помощью рекурсивной процедуры - C (СИ)

float teylor(float x, int n, float lastVal, float eps) {
    float val = lastVal*x*x / ((2*n-1) * 2*n);
    lastVal += val;
    if (fabs(val) < eps)
        return lastVal;
    else
        return teylor(x, n+1, lastVal, eps);
}