/* 0th order modified bessel function of the first kind.  */
static double bessel(double x)
/* 0th order modified bessel function of the first kind.  */
static double bessel(double x)
{
        lastv=v;
        t *= x/(i*i);
        v += t;
    }
    return v;
}