#include <stdio.h>
#include <stdlib.h>
#include <math.h>

#include "bnc.h"

double dfunc(double x)
{
	return x * x;
}

int main()
{
	double ans, start, end;

	start = 0.0;
	end = 1.0;

	ans = dtrapezoidal_fs(start, end, dfunc, 2);
	printf("%e\n", ans);
	ans = dtrapezoidal_fs(start, end, dfunc, 4);
	printf("%e\n", ans);
	ans = dtrapezoidal_fs(start, end, dfunc, 8);
	printf("%e\n", ans);
	ans = dtrapezoidal_fs(start, end, dfunc, 16);
	printf("%e\n", ans);

	return EXIT_SUCCESS;
}