hyper2geom1.cpp

Go to the documentation of this file.

#include "tbci/cplx.h"
#include "tbci/specfun.h"
#include <iostream>

NAMESPACE_TBCI
TBCI__ cplx<double> hyper2geom1 (const TBCI__ cplx<double> a,
                                 const TBCI__ cplx<double> b,
                                 const TBCI__ cplx<double> c,
                                 const TBCI__ cplx<double> z)
{
        TBCI__ cplx<double> oldres, res, fac;
        res = 1.0; fac = 1.0;
        for (int n = 0; n < 120; n++) {
                oldres = res;
                fac *= (a+n) * (b+n) * z / ((c+n) * (n+1));
                res += fac;
                if (res == oldres)
                        return res;
        }
        STD__ cerr << "hyper2geo1: Accuracy not reached!" << STD__ endl;
        return res;
}
NAMESPACE_END