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

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

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

На СИ, к сожалению, не пишу, а вот "что к чему" - попробую объяснить.
Искомый интеграл представляет собой сумму площадей вертикально стоящих прямоугольных трапеций с основаниями Y(i), Y(i+1) и высотой (X(i+1)-X(i)). Вот и напишите цикл вычисления и суммирования площадей этих трапеций.