/***********************************************************************
 * Soluciones propuestas para resolver el examen de ATDI del día 210907* 
 * Francisco Javier Garcia Castellano 2007                             *
 ***********************************************************************/

/***********************************************************************/ 
//BIBLOTECAS
#include <stdio.h>
#include <math.h>


/***********************************************************************/ 
//CONSTANTES
#define MAX 50 //Maximo tamaño para vectores 


/***********************************************************************/ 
//FUNCIONES

/*---------------------------------------------------------------------*/
/** Procedimiento que calcula la covarianza entre dos series de numeros.
Parametros:
serie1: vector de numeros reales con la primera serie. 
serie2: vector de numeros reales con la segunda serie. 
n: numero de elementos de las series
*/
double covarianza (double serie1[], double serie2[], int n) {

  int i,j; /*Iteradores*/
  double suma1=0,suma2=0; /*para almacenar sumatorias*/
  double media1,media2; /*medias para las dos series*/

  /*Calculamos la sumatoria de ambas series*/
  for (i=0;i<n;i++) {
    suma1+=serie1[i];
    suma2+=serie2[i];
  }

  /*Con la suma de los elementos de cada serie, calculamos las medias*/
  media1=suma1/(double)n;
  media2=suma2/(double)n;
  printf("media1= %lf media2=%lf\n",media1,media2);

  /*Calculamos la sumatoria de la covarianza*/
  for (j=0,suma1=0;j<n;j++) 
    suma1+=(serie1[i]-media1)*(serie2[i]-media2);

  /*Devolvemos la covarianza*/
  return suma1/(double)n;

}/*Fin función covarianza*/

/*---------------------------------------------------------------------*/
/** Función main para el ejercicio de las pesas. Esta función calcula el 
número mínimo de pesas de 50Kg, 20Kg, 10Kg, 5kg y 1Kg que se necesitan en
una balanza para un peso introducido por el usuario.
*/
int main_pesas(){
  int peso; /*peso introducido por el usuario*/

  /*Le pedimos el peso al usuario*/
  do {
    printf ("\nDame el peso (mayor que 0 y positivo):");
    scanf("%d",&peso);
  }while (peso <=0);

  /*Mostramos las pesas de 50Kg  a usar*/
  printf("\nHay que usar %d pesos de 50kg",peso/50);
  peso=peso%50;

  /*Mostramos las pesas de 20Kg a usar, una vez quitado el 
    peso usado por las de 50Kg*/
  printf("\nHay que usar %d pesos de 20kg",peso/20);
  peso=peso%20;

  /*Mostramos las pesas de 10Kg  a usar*/
  printf("\nHay que usar %d pesos de 10kg",peso/10);
  peso=peso%10;

  /*Mostramos las pesas de 5Kg  a usar*/
  printf("\nHay que usar %d pesos de 5kg",peso/5);
  peso=peso%5;

  /*Mostramos las pesas de 1g a usar, son las qeu quedan*/
  printf("\nHay que usar %d pesos de 1kg\n",peso);
 
  /*Todo ha termiando bien*/
  return 0;

}/*Fin main ejercicio pesas*/

/*---------------------------------------------------------------------*/
/**Procedimiento que muestra el polinomio.
Parámetros:
polinomio: vector con los coeficientes del polinomio.
n: grado del polinomio.*/
void muestra(double polinomio[], int n){
  int i; /*Iterador*/

  /*Mostramos cada uno de los elementos del polinomio,
    mostrando el coeficiente y la x con su exponente*/
  for (i=n-1;i>0;i--) 
    printf ("%lf x ^%d + ",polinomio[i],i);
  
  /*Mostramos el término independiente*/
  printf("%lf \n",polinomio[0]);

}/*Fin función muestra polinomio*/

/*---------------------------------------------------------------------*/
/**Función que calcula la derivada de un polinomio y devuelve el grado de la derivada.
Parámetros:
polinomio: vector con los coeficientes del polinomio.
n: grado del polinomio.
derivada: vector donde se guardará los coeficientes de la derivada del polinomio.
*/
int derivada(double polinomio[], int n, double derivada[]){
 int i; /*Iterador*/

   for (i=n-1;i>0;i--) 
     derivada[i-1]=polinomio[i]*i;
 
   return n-1;
  
}/*Fin función que deriva un polinomio*/

/*---------------------------------------------------------------------*/
/**Función que devuelve la evaluación de un polinomio en un punto.
Parámetros:
polinomio: vector con los coeficientes del polinomio.
n: grado del polinomio.
x: punto en el cual se evalúa el polinomio. */
double evalua(double polinomio[], int n, double x){

 int i; /*Iterador*/
 double evaluacion=0; /*Var donde alamacenaremos la suma de la evaluación
		     de los distintos elementos del polinomio*/

 /*Evaluamos cada uno de los términos del polinomio, dándole
   valor a la x*/
 for (i=n-1;i>0;i--) 
   evaluacion+=polinomio[i]*pow(x,i);
 
 /*Añadimos el término independiente*/
 evaluacion+=polinomio[0];
 
 /*Devolvemos el valor calculado*/
 return evaluacion;

}/*Fin función que evaluar un polinomio en un punto*/

/*---------------------------------------------------------------------*/
/*---------------------------------------------------------------------*/
int main(){
  double polinomio[MAX]; /*Vector donde almacenaremos los coeficientes del
			  polinomio*/
  double deriv[MAX]; /*Vector donde almacenaremos la derivada*/
  int n; /*Tamaño del polinomio*/
  int m; /*Tamaño de la derivada*/
  int i; /*Iterador*/
  double x; /*Variable para guardar un valor dado por el usuario*/
  int opcion; /*Opción escogida del menú*/


  /*Pedimos el tamaño del polinomio*/
  do{
    printf ("Dame el el grado del polinomio:");
    scanf("%d",&n);

    /*Le añadimos uno por el término independiente*/
    n++;
  }while  (n<1);

  /*Leemos los coeficientes del polinomio*/
  for (i=0;i<n;i++) {
    printf ("\nDame el coeficiente %d:",i);
    scanf("%lf",&(polinomio[i]) );
  }
  
  /*Mostramos el polinomio introducido*/
  printf("\n\nEl polinomio queda de la siguiente forma:\n");
  muestra (polinomio,n);

  /*Precalculamos la derivada*/
  m=derivada(polinomio,n,deriv);  

  /*MENU*/
  do{
    /*Mostramos el menú*/
    printf("\n----------MENU---------------\n");
    printf("\n1 Mostrar polinomio\n");
    printf("2 Mostrar derivada\n");
    printf("3 Evaluar polinomio en un punto\n");
    printf("4 Evaluar derivada en un punto\n");
    printf("5 Salir\n");
    printf("Opción:");

    /*Leemos la opción del usuario*/
    scanf("%d",&opcion);
    printf("\n");
   
    /*Dependiendo de la opción escogida efectuamos una opción u otra*/
    switch (opcion) {
      /*Primera opción: Mostrar el polinomio*/
      case 1: 
	printf("\nEl polinomio es::\n");
	muestra(polinomio,n);
	break;

      /*Segunda opción: Calcular la derivada y mostrarla*/
      case 2: 
	printf("\n La derivada es:\n");
	muestra(deriv, m);
	break;
	
      /*Tercera opción: Evaluar en un punto*/
      case 3:
	printf("\nDame el punto sobre el cual evaluar el polinomio:");
	scanf("%lf",&x);
  	printf("El resultado de evaluar %lf es =%lf",x,evalua(polinomio,n,x));
	break;

      /*cuarta opción: Evaluar en un punto la derivada*/
      case 4: 
	printf("\nDame el punto sobre el cual evaluar la derivada:");
	scanf("%lf",&x);
  	printf("El resultado de evaluar %lf es =%lf",x,evalua(deriv,m,x));
	break;

      /*Quinta opción: salir*/
      case 5: 
	break;

      /*En otro caso es opción no válida*/
      default:
	printf("Opción Inválida\n");
     
    }/*Fin switch*/

  } while (opcion!=5);

  /*Pruebo el ejercicio de la covarianza*/
  printf("Covarianza=%lf",covarianza(polinomio,deriv,m));
  
  /*Pruebo el ejercicio de las pesas*/
  main_pesas();  

  /*Terminamos*/
  return 0;
}
/*---------------------------------------------------------------------*/