#include <stdio.h>
#include <stdlib.h>
#include <conio.h>
#include <math.h>
#include <ctype.h>
float factorial(float);
float combinacion(float, float);
main()
{
float x, x_1, x_2, number;
float n, p, q, suma = 0;
float binomial = 0, sumatoria = 0;
char opcion, salida;
do{
system("cls");
printf("\tDistribucion binomial \n\n");
printf("[a]-Funcion individual.\n");
printf("[b]-Funcion acumulada.\n\n");
printf("*Selecciona la funcion que desee calcular: ");
scanf("%s", &opcion);
switch(opcion)
{
case 'a': printf("\n\n\a-Funcion individual de distribucion binomial f(x).-\n");
printf("*Ingrese el numero de eventos 'n': ");
scanf("%f", &n);
do{
printf("\n*Ingrese la probabilidad de exito %: ");
scanf("%f", &p);
if(p > 1)
p /= 100;
q = 1 - p;
suma = p + q;
}while(suma != 1);
do{
printf("\n*Ingrese el valor de 'x': ");
scanf("%f", &x);
if(x > n)
printf("ERROR 658: 'x' no puede ser mayor a 'n'\n");
if(x < 0)
printf("ERROR 725: 'x' no puede ser menor a cero.\n");
}while( -((x > n) || (x < 0)) );
printf("\n\n-Resultados del modelo binomial.-\n\n");
binomial = combinacion(n, x) * pow(p, x) * pow(q, n-x);
printf(" B(%2.0f, %2.0f, %2.0f) = %4.6f = %4.2f %", x, n, p, binomial, binomial * 100);
break;
case 'b': printf("\n\n\a-Funcion acumulada de distribucion binomial F(x).-\n");
printf("*Ingrese el numero de eventos ´n´");
scanf("%f", &n);
do{
printf("\n*Ingrese la probabilidad de exito %: ");
scanf("%f", &p);
if(p > 1)
p /= 100;
q = 1 - p;
suma = p + q;
}while(suma != 1);
do{
printf("*Ingresa x_1: ");
scanf("%f", &number);
x_1 = number;
if(x_1 < 0)
printf("\nERROR 458: El intervalor x_1 no puede ser menor a cero.\n");
if(x_1 > n)
printf("\nERROR 218: El intervalo x_1 no puede ser mayor a 'n'.\n\n");
}while( -( (x_1 < 0)||(x_1 > n) ) );
do{
printf("*Ingresa x_2: ");
scanf("%f", &x_2);
if(x_2 < x_1)
printf("\nERROR 147: El intervalo x_2 no puede ser menor a x_1.\n");
if(x_2 > n)
printf("\nERROR 231: El intervalo x_2 no puede ser mayor a 'n'.\n");
}while( -( (x_2 < x_1)||(x_2 > n) ) );
printf("\n\n-Resultados del modelo binomial.-\n\n");
for(x_1 = number; x_1 <= x_2; x_1++){
binomial = combinacion(n, x_1) * pow(p, x_1) * pow(q, n - x_1);
printf(" B(%2.0f, %2.0f, %2.2f) = %4.6f\n", x_1, n, p, binomial);
sumatoria += binomial;
}
printf("\n Probabilidad final.");
printf("\n F(x) = (x; %2.0f, %2.0f) = %4.4f = %4.2f%", n, p, sumatoria, sumatoria * 100);
break;
default: printf("\nERROR 254: caracter de opcion mal introducido, intentelo de nuevo.\n");
break;
}//fin del menu opcion
printf("\n\n¿Desea un nuevo calculo(S/N)?");
salida = getch();
salida = tolower(salida);
binomial = 0;
sumatoria = 0;
}while(salida == 's');
system("cls");
printf("\tModelo binomial \n\n");
printf("\n\nVuelve pronto...");
getch();
return 0;
}//fin de main
float combinacion(float a, float b)
{
float combinacion;
combinacion = ( factorial(a) /( factorial(b) * factorial(a-b) ));
return combinacion;
}
float factorial(float number)
{
float factorial = 1, counter;
for(counter = number; counter >= 1; counter--)
factorial *= counter;
return factorial;
}
