Powered By Blogger

martes, 11 de diciembre de 2012

METODO DIRECTO DE LOS TRAPECIOS


//Torres Martinez Gustavo Alan
//Grupo:4CM1
//Programa que realiza la solucion de una integral por el metodo
//directo de los TRAPECIOS.

#include<stdio.h>
#include<conio.h>
#include<dos.h>
#include<math.h>
#include<stdlib.h>
#include<iostream.h>

void main(void)
{
int i;
double a,b,n,in,h,vv=0.3413447461,e,ea,eb,ev,erv,f,sf=0.0,fx,fxa,fxb;
clrscr();
cout<<"Dame el valor de n:";
cin>>n;
cout<<endl<<"Dame el valor del l°mite inferior:";
cin>>a;
cout<<endl<<"Dame el valor del l°mite superior:";
cin>>b;
h=(b-a)/n;
printf("\nValor de h:%1.2f\n",h);
getch();
clrscr();
//e=a+(1*h);
//f=exp(e);
//printf("%2.8f",f);
for(i=1.0;i<n;i++)
{
e=pow(a+i*h,2)/2.0;
f=exp(-e);
sf=sf+f;
e=(a+i*h);
printf("\nValor de la funcion evaluada en (%0.2f):%0.8f\n",e,f);
if(i%24==0&&i!=0){getch(); clrscr();}
}
in=2*sf;        //printf("%1.8f",in);
ea=pow(a,2)/2.0;
fxa=exp(-ea);
eb=pow(b,2)/2.0;
fxb=exp(-eb);
fx=(h/(2.0*sqrt(6.28318531))*(fxa+fxb+in));
printf("\n\nValor de f(x):%1.8f\n",fx);
printf("\nValor verdadero de la funci¢n:%0.10f\n",vv);
ev=vv-fx;
erv=(ev/vv)*100;
printf("\nEv=%0.8f\tErv=%2.3f%",fabs(ev),fabs(erv));
cout<<endl<<endl;
textcolor(10+128);
cprintf("             Presiona cualquier tecla....");
getch();
}

No hay comentarios:

Publicar un comentario