#include <iostream>
#include <kv/defint-singular.hpp>

typedef kv::interval<double> itv;

struct Quadrant {
	template <class T> T operator() (const T& x) {
		return sqrt(1 - x * x);
	}
};

struct Quadrant_s {
	template <class T> T operator() (const T& x) {
		return 1 - x * x;
	}
};

int main() {
	std::cout.precision(17);

	std::cout << kv::defint_power2_autostep_r(Quadrant(), Quadrant_s(), (itv)0., (itv)1., 12, (itv)0.5) << "\n";
}