Método dos trapézios compostos

Método dos trapézios compostos

#include<iostream>
using namespace std;
float F(float x){
return 1/x;
}
int main(){
    float a; float b;
    int n,i;
    float I; float H;
    float soma1=0,soma2=0;
    cout<<"Digite o limite de integracao a esquerda: ";
    cin>>a;
    cout<<"Digite o limite de integracao a direita: ";
    cin>>b;
    cout<<"Digite a quantidade de subintervalos: ";
    cin>>n;
    H=(b-a)/(float)n;
    float p=a;
    for(i=1;i<n;i++){
        soma1+=F(p+H);
        p=p+H;
    }
I=(H/2)*(F(a)+(2*soma1)+F(b));
cout<<"A integral de "<<a<<" ate "<<b<<" de (1/x dx) correspode a:"<<I<<endl;
}

Comentários