double besselpoly(double a, double lambda, double nu);