hyper2geom1_stdcplx.cpp

Go to the documentation of this file.

#include "tbci/basics.h"
#include "tbci/specfun_stdcplx.h"
#include <iostream>

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