#include "interval.h"

int main()
{
	interval a, b, c, d, r, t1, t2;

	a = interval_init(1);
	b = interval_init(1e15);
	c = interval_init(1e14);

	/* r = (-b + sqrt(b * b - 4 * a * c)) / (2 * a); */

	t1 = interval_mul(b, b);
	t2 = interval_mul(interval_mul(interval_init(4), a), c);
	d = interval_sub(t1, t2);

	d = interval_sqrt(d);

	t1 = interval_sub(d, b);
	t2 = interval_mul(interval_init(2), a);

	r = interval_div(t1, t2);

	interval_print(r);

	/* r = 2 * c / (-b - sqrt(b * b - 4 * a * c)); */

	t1 = interval_mul(b, b);
	t2 = interval_mul(interval_mul(interval_init(4), a), c);
	d = interval_sub(t1, t2);
	d = interval_sqrt(d);

	t1 = interval_mul(interval_init(2), c);
	t2 = interval_minus(interval_add(b, d));

	r = interval_div(t1, t2);

	interval_print(r);
}