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

[Serjuk]

Здравствуйте, возникла проблемка написания программы вычисления интеграла. У меня есть наклонная прямая, на которой находятся точки. По этой прямой действует синусоидальная функция. q(ksi) = 1000 * sin (Z*ksi), где ksi - от 0 до конца наклонной прямой(координаты конечной точки и координаты начальной точки есть, поэтому вычислить кси не составит труда ),а z - const, ее значение есть.
И надо найти интеграл методом трапеций.
Вот код, что я смог сделать.. По сути просто нашел точки.

* N - количество точек  
DO 10 I=1,N
* Обрабатываемые точки
        XT=MASSIV(NC*(I-1)+1)
        YT=MASSIV(NC*(I-1)+2)
             
  * находим все точки лежащие на наклонной прямой
      EQUAL = 23*XT + 14*YT - 2030          
      IF (XT.GE.0 .AND. XT.LE.35 .AND. ABS(EQUAL).LT.0.01) THEN 

* все точки записываются в массив F
* номер точки
            I1=I1+1
            F(I1)=I

            NR=NR+1
            ENDIF

А как дальше вычислить интеграл, ведь в каждой точке по сути действует функция?