#include <oxstd.h>

#include <lib/quantile.ox>

ProbN(const x, ...)
{
	return probn(x);
}
DensN(const x, ...)
{
	return densn(x);
}

main()
{
    decl t = range(1,10);

    print("\nF-distribution", "%c", {"t", "probf(t)", "quanf(probf(t))"},
		"%20.10g", t' ~ probf(t,10,10)'
		~ Quantile(probf(t,10,10),10,10, densf, probf)' );

	t -= 5;
    print("\nNormal", "%c", {"t", "probn(t)", "quann(probn(t))"},
		"%20.10g", t' ~ probn(t)'
		~ Quantile(probn(t),0,0, DensN, ProbN)' );
}

