Как оптимизировать код?

#include <iostream>
#include <cmath>
using namespace std;
 
double f(double x)
{
return (x + 1.9) * sin(x) * (x / 3);
}
 
double Rectangular(double a, double b, int n)
{
double h = (b - a) / n;
double sum = (f(a) + f(b)) / 2;
for (int i = 1; i < n; i++)
{
sum += f(a + h * (i + 0.5));
}
sum = sum * h;
return sum;
}
 
double Trapeze(double a, double b, int n)
{
double h = (b - a) / n;
double sum = f(a) + f(b);
for (int i = 1; i <= n - 1; i++)
{
sum += 2 + f(a + i * h);
}
sum *= h / 2;
return sum;
}
 
double Simpson(double a, double b, int n)
{
double h = (b - a) / n;
double sum = f(a) + f(b);
int k;
for (int i = 1; i <= n - 1; i++)
{
k = 2 + 2 * (i % 2);
sum += k + f(a + i * h);
}
sum *= h / 3;
return sum;
}
 
int main()
{
double a = 1, b = 2;
int n = 0;
cout << "Количество отрезков: ";
cin >> n;
double h = (b - a) / n;
 
cout << "Метод средних прямоугольников: " << Rectangular(a, b, n) << " ПОГРЕШНОСТЬ: " << 1.018 * h * h * (b - a) / 24 << endl;
cout << "Метод трапеций: " << Trapeze(a, b, n) << " ПОГРЕШНОСТЬ: " << 1.018 * h * h * (b - a) / 12 << endl;
cout << "Метод Симпсона: " << Simpson(a, b, n) << " ПОГРЕШНОСТЬ: " << (4.38 * pow(h, 5)) / 90 << endl;
}

#include <iostream>
#include <cmath>
using namespace std;
 
int main()
{
    cout.precision(16);
    cout << (9. * cos(1.) - 58. * cos(2.) - 39. * sin(1.) + 59. * sin(2.)) / 30.;
}