double calnum(double (*f)(double), double a, double b, int n)
{
       int i;
       double krok = (b-a)/n; //wysokość każdego z trapezów
       double wynik = 0;
       for (i=0; i<n; i++) //zaczynam od zerowego wyrazu
           wynik += (f(x)*f(x+1)krok)/2;
           x += krok;
           }
           return wynik;
}