Вычисление интеграла на С++

[ExZent]

На плоскости своими координатами xk, yk заданы n точек r0, . . . , rn−1, причем
xk строго монотонно возрастают. Заданные точки определяют на отрезке x ∈
[x0, xn−1] кусочно-линейную функцию f(x). Вычислить определенный интеграл функции f(x) в пределе [x0, xn−1]

Помогите пожалуйста решить задачу, на языке С++.
Ну или покажите пример решения на бумаге, что бы можно было понять что к чему