#include <oxstd.h>
#include <quadpack.h>

namespace gauss
{

gausscall2ox_func(const a, const b)
{
	return a + unit(b);
}
qng(const func, const a, const b)
{
	decl result, abserr;
	println("In Ox: qng(): func=", func, " a=", a, " b=", b);
	::QNG(func, a, b, &result, &abserr);
	return {result, abserr};
}

}