#include<iostream>
#include<cmath>
using namespace std;
double f(double x) {
return (cos(x) + x);
}
double S(int a, int b) {
return (f(a) + f(b)) / 2;
}
int main() {
int a=1, b=4, m=1500;
double h = 3/ (double)m;
double c = 0;
for (int i = 1; i < m; i++) {
c += f(a + i * h);
}
double result;
result = (S(a,b)+ c) * h;
cout << result;
return 0;
}
输出样例:5.90173