Найти приближенное значение определенного интеграла по методу прямоугольников с точностью e=10^-3

#include "stdafx.h"
#include <iostream>
#include <math.h>
#include <Windows.h>
 
using namespace std;
 
double f1(double a, double b, int n,double(*f)(double) )
{
    double f_temp,h;
    int i;
    f_temp=0.0;
    for(i=1;i<=n;i++)
    {
        h=(b-a)/n;
        f_temp+=f(a+h*(i-0.5));
    }
    return f_temp;
}
double integral(double x)
{
    return pow(x+1,2)*pow(2,x);
}
int main()
{
    setlocale(LC_ALL,"Russian");
    double S;
    S=f1(0,1,10,integral);
    cout<<"Значение интеграла: "<<"((x+1)^2)*2^x"<<"\n"<<"По методу прямоугольников равно: "<<S<<"\n";
    return 0;
}

#include "stdafx.h"
#include <iostream>
#include <math.h>
#include <Windows.h>
 
using namespace std;
 
int main()
{
    setlocale(LC_ALL,"Russian");
    double a,b,f,f_temp,h,E,in,x,S;
    int n,i;
    E=0.001;
    f_temp=0;
    x=0;
    S=0;
    cout<<"Введите a: ";cin>>a;
    cout<<"Введите b: ";cin>>b;
    cout<<"Введите кол-во шагов n: ";cin>>n;
    for(i=1;i<=n;i++)
    {
        
        h=(b-a)/n;
        f=(a+h*(i-0.5));
        f_temp+=f;
    }
    in=pow(x+1,2)*pow(2,x);
    
    return 0;
}