#include <stdio.h>

#include "betaDF.h"

void main(void)
{
	char str[10];
	double trials, successes, p, x1, x2, y, level = 0.95;

   printf("\nChecking Binomial:\n");
   while(1) {
      printf("\nEnter number of trials (0 to exit): ");
		gets(str);      
		sscanf(str, "%lf", &trials);
		if (trials == 0) break;
      printf("Enter number of successes: ");
		gets(str);      
		sscanf(str, "%lf", &successes);
      printf("Enter probability of success: ");
		gets(str);      
		sscanf(str, "%lf", &p);
		x1 = binomialDF(trials, successes, p);
		x2 = binomialDF(trials, successes-1, p);
		printf("\tBinomial(<=%.0lf|%.0lf,%.4f)=%lg\n", successes, trials, p , x1);
		printf("\tBinomial(==%.0lf|%.0lf,%.4f)=%lg\n", successes, trials, p , x1-x2);
		y = rev_binom(trials, successes, x1);
		printf("\tEstimated p=%lg\n", y);
		y = binom_leftCI(trials, successes, level);
		printf("\tCI(%.2lg): %lg", level, y);
		y = binom_rightCI(trials, successes, level);
		printf(" <= %.2lg <= %lg\n", successes/trials, y);
   }
}