Метод Симпсона С++

[sabaton]

Доброго времени суток! Не работает программка. Нужно вычислить интеграл методом Симпсона с заданным шагом h разбиений интервала a, b

Вот что имеем:
Когда ввожу шаг h=0.1, при a=0, b=1 программа показывает, что интеграл равен 0.. Но это неверно..

 
#include <math.h>
#include <iostream>
#include <conio.h>

 
using namespace std;
 
//Подынтегральная функция
double func(double x)
{
    double r;
    r = 1.5*x/2+cos(2*x);
    return r;
}
double simpson( double (*Fx)(double), double a, double b, int h);
 

int main(int argc, char* argv[])
{
   double I;
   double a,b;
   int h;
 
   cout << "\n Введите парамерты:\n" << endl;
   cout << "\t-> Enter a = ";
   cin >> a;
 
   cout << "\t-> Enter b = ";
   cin >> b;
 
   cout << "\t-> Enter h = ";
   cin >> h;
 
   // вычесление
 
   I = simpson( &func, a, b, h );
 
   // вывод результата
   cout << endl << "      I = " << I;
   getch();
}
 
//---------------------------------------------------------------------------
// сам метод
double simpson( double (*Fx)(double), double a, double b, int h)
{
 
   
   
 
   double I, I2 = 0, I4 = 0;
   I4 = Fx(a+h);
   for(int k = 2; k <=b  ; k += 2 )   // Думаю ошибка в этом b , но не знаю что нужно на самом деле
   {
       I4 += Fx(a+(k+1)*h);
       I2 += Fx(a+k*h);
   }
   I = Fx(a)+Fx(b)+4*I4+2*I2;
   I *= h/3;
 
   return I;
}

Подскажите, буду благодарен)

[Mandrivnyk]

Когда ввожу шаг h=0.1

int h

Поменяй тип у h -)

[sabaton]

блин...вот это тормоз))) сцука..
Спасибо) а то целый день уже сижу и не знаю, что делать))