Вычислить сумму по формуле с использованием функции

#include <iostream>
 
double sum(double x, double y, double z)
{
    return (x+y)/z; 
}
 
int main()
{
    //Так как величины, на кот. изменяется каждая переменная нам не известны, придется самим посчитать...
    double sum_f=0.0;//суммы 1,2,3 части.
    double sum_s=0.0;
    double sum_t=0.0;
    int n = 15; 
    
    double x_s = 0.2; 
    double x_e = 1.4;
    double x_r = (x_e-x_s)/(n-1);//прирост x
    std::cout << " x_r = " <<x_r;
    
    double y_s = 1; 
    double y_e = 15;
    double y_r = (y_e-y_s)/(n-1); //prirost y
    std::cout << " y_r = " <<y_r;
    
    double c_s = 0.1; 
    double c_e = 2.9;
    double c_r = (c_e - c_s)/(n-1);
    std::cout << " c_r = " <<c_r;
    
    double d_s = 1.0;
    double d_e = 3.8;
    double d_r = (d_e-d_s)/(n-1);
    std::cout << " d_r = " <<d_r;
    
    for(int i = 1; i<n+1; i++, x_s+=x_r, y_s+=y_r,c_s+=c_r, d_s+=d_r)
    {
        sum_f += x_s*y_s;   
        sum_s += c_s*d_s;
        sum_t += x_s*c_s;
    }
    std::cout <<"\nResult = " << sum(sum_f,sum_s,sum_t);
}