practica9_6.cpp

Topics Practices: Programming and Numerical Methods Practice 1: Introduction to C Practice 2

/*PRACTICA 9 EJERCICIO 6
Sea f(x) una funci髇 acotada en el intervalo [a,b], y supongamos que f(x) toma 
valores entre c y d en este intervalo. Generalizar el m閠odo del ejercicio 
anterior para calcular la siguiente integral en forma num閞ica:
         
         INTEGRAL ENTRE A Y B DE F(X)
*/

#include<stdio.h>
#include<conio.h>
#include<math.h>
#include<stdlib.h> //biblioteca para el rand y el RAND_MAX..

float f(float x)
{
      return sqrt(x); // Integral entre cero y uno de raiz de x dx...
}

main()
{
      float x,y,integral,a,b,d;
      int n,cant,i;
      
      n=100000; //Cantidad de iteraciones que voy a hacer; cantidad de repeticiones que voy a hacer para calcular la integral..    
      
      cant=0;
      for (i=0; i<n; ++i)
      {
          x=a + (b-a)*rand()/(double)RAND_MAX; //Divido por RAN_MAX; por que retoma con un valor entero entre cero y RAN_MAX(valor entero grande)
          y=d*rand()/(double)RAND_MAX;
          
          if (0<f(x)<d)
             cant++;          
      }
 integral= cant/(double) n; //convierto uno de los terminos a real y me aseguro que asi sea (double)..
 
 printf("El valor estimado de la integral es %f \n\n", integral);
 printf("PRECIONE CUAQUIER TECLA PARA TERMINAR!!!!");
 getch();     
   
}