Решение пользовательского интеграла методом трапеций - C (СИ)

Формулировка задачи:

Доброго времени суток, товарищи программисты. Задача следующая:решение пользовательского интеграла методом трапеций. Загвоздка, собственно, в самой формуле-не могу понять что в ней не так.

                    Листинг программы
                
 
                        
                          
                    
#include <stdio.h>
#include <math.h>
#include <locale.h>
#include <stdlib.h>
#include "zagolovoc.h"
 
void main()
 
{   
    float a,b,rez,h,k,x[200],n,rez1,rez2;
    int i,t;
    setlocale(LC_ALL,"russian");
    printf("Введите A, B и точность вычисления\n");
    scanf("%f",&a);
    scanf("%f",&b);
    scanf("%f",&n);
    t=(int)n;
    h=(b-a)/n;
    for(i=1;i<t-1;i=i+h)
    { 
        rez=h*((function(a)+function(b))/2+(float)function(x[i]));
    }
 
    printf("Метод Трацеций результат %.4f\n",rez);
}
 
описание функции
 
#include <stdio.h>
 #include <math.h>
 float function(float x)
{
 return tan(x);
 }

Решение задачи: «Решение пользовательского интеграла методом трапеций»

textual
                            Листинг программы