{

try

{ // Try a zero order v.

float dodgy_root = boost::math::cyl_bessel_j_zero(0.F, 0);

std::cout << "boost::math::cyl_bessel_j_zero(0.F, 0) " << dodgy_root << std::endl;

// Thrown exception Error in function boost::math::cyl_bessel_j_zero<double>(double, int):

// Requested the 0'th zero of J0, but the rank must be > 0 !

}

catch (std::exception& ex)

{

std::cout << "Thrown exception " << ex.what() << std::endl;

}

try

{ // order v = inf

std::cout << "boost::math::cyl_bessel_j_zero(inf, 1) " << std::endl;

double inf = std::numeric_limits<double>::infinity();

double inf_root = boost::math::cyl_bessel_j_zero(inf, 1);

std::cout << "boost::math::cyl_bessel_j_zero(inf, 1) " << inf_root << std::endl;

// Throw exception Error in function boost::math::cyl_bessel_j_zero<long double>(long double, unsigned):

// Order argument is 1.#INF, but must be finite >= 0 !

}

catch (std::exception& ex)

{

std::cout << "Thrown exception " << ex.what() << std::endl;

}

try

{ // order v = NaN, rank m = 1

std::cout << "boost::math::cyl_bessel_j_zero(nan, 1) " << std::endl;

double nan = std::numeric_limits<double>::quiet_NaN();

double nan_root = boost::math::cyl_bessel_j_zero(nan, 1);

std::cout << "boost::math::cyl_bessel_j_zero(nan, 1) " << nan_root << std::endl;

// Throw exception Error in function boost::math::cyl_bessel_j_zero<long double>(long double, unsigned):

// Order argument is 1.#QNAN, but must be finite >= 0 !

}

catch (std::exception& ex)

{

std::cout << "Thrown exception " << ex.what() << std::endl;

}

try

{ // Try a zero rank m.

std::cout << "boost::math::cyl_neumann_zero(0.0, 0) " << std::endl;

double dodgy_root = boost::math::cyl_bessel_j_zero(0.0, 0);

// warning C4146: unary minus operator applied to unsigned type, result still unsigned.

std::cout << "boost::math::cyl_neumann_zero(0.0, -1) " << dodgy_root << std::endl;

// boost::math::cyl_neumann_zero(0.0, -1) 6.74652e+009

// This *should* fail because m is unreasonably large.

}

catch (std::exception& ex)

{

std::cout << "Thrown exception " << ex.what() << std::endl;

}

try

{ // m = inf

std::cout << "boost::math::cyl_bessel_j_zero(0.0, inf) " << std::endl;

double inf = std::numeric_limits<double>::infinity();

double inf_root = boost::math::cyl_bessel_j_zero(0.0, inf);

// warning C4244: 'argument' : conversion from 'double' to 'int', possible loss of data.

std::cout << "boost::math::cyl_bessel_j_zero(0.0, inf) " << inf_root << std::endl;

// Throw exception Error in function boost::math::cyl_bessel_j_zero<long double>(long double, int):

// Requested the 0'th zero, but must be > 0 !

}

catch (std::exception& ex)

{

std::cout << "Thrown exception " << ex.what() << std::endl;

}

try

{ // m = NaN

double nan = std::numeric_limits<double>::quiet_NaN();

double nan_root = boost::math::airy_ai_zero<double>(nan);

// warning C4244: 'argument' : conversion from 'double' to 'int', possible loss of data.

std::cout << "boost::math::airy_ai_zero<double>(nan) " << nan_root << std::endl;

// Thrown exception Error in function boost::math::airy_ai_zero<double>(double,double):

// The requested rank of the zero is 0, but must be 1 or more !

}

catch (std::exception& ex)

{

std::cout << "Thrown exception " << ex.what() << std::endl;

}