practica5_1.cpp

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

/*PRACTICA 5 EJERCICIO 1
Escribir un programa que le pida al usuario dos n鷐eros reales a y b, junto con 
el paso de integraci髇 h > 0, y que calcule la siguiente integral por el m閠odo 
de suma de rect醤gulos, el m閠odo del trapezoide y el m閠odo de Simpson:
        
        I = Integral entre a y b de sqrt(x)*sen(x)
        
Comparar la precisi髇 de los tres m閠odos en funci髇 del paso de integraci髇.*/
        

#include<stdio.h>
#include<conio.h>
#include<math.h>

float f(float x)
{
 return sqrt(x)*sin(x);      
}


float g(float a,float b,float h)
{
float i,suma;
suma=0;
for(i=a;i<=b;i=i+h)
suma=suma + h*f(i);//Suma de rectangulos..
return suma;
}


float t(float a,float b,float h)
{
float i,suma;
suma=0;
for(i=a;i<=b;i=i+h)
suma=suma + (f(i)+ f(i+h))*h/2;//Metodo de trapezoide...
return suma;   
}


float p(float a, float b, float h)
{
float i,suma;
suma=0;
for(i=a;i<=b;i=i+2*h)
suma=suma + (f(i)+ 4*f(i+h) + f(i+2*h))*h/3;//Metodo Simpson..
return suma;   
}


main()
{
      float a,b,h,m,n,z;      
      printf("Calcule la integral entre a y b de sqrt(x)*sen(x) \n");
      printf("Ingrese el valor de el extremo izquierdo(a): ");
      scanf("%f", &a);
      printf("Ingrese el valor de el extremo derecho(b): ");
      scanf("%f", &b);
      printf("Ingrese el paso: ");
      scanf("%f", &h);
      m=g(a,b,h);      
      n=t(a,b,h);
      z=p(a,b,h);
      printf("La integral calculada con Suma de rectangulos  dio %f \n", m);
      printf("La integral calculada con Metodo de trapezoide dio %f \n", n);
      printf("La integral calculada con Metodo Simpson       dio %f \n\n", z);
      printf("PRESIONE CUALQUIER TECLA PARA TERMINAR!!!!!");      
      
      getch();
}