Проверьте правильность решения на C++ (вычислить интеграл)

[mov]

Задание

Решение

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

double f (double t,double u)
{
if(u>0) return t+u;
else if (u<=0) return t*u;
else return -1;
}

main ()
{
double j,h,e=0,i=1,a=3,b=8,n=10;
h=(b-a)/n;

while (i<10)
{
e+=f(pow(a+i*h,3)/(pow(a+i*h,4)+1),sin(pow(a+i*h,2))) ;
i++;
}

j=h*((f(pow(a,3)/(pow(a,4)+1),sin(pow(a,2))) +
f(pow(b,3)/(pow(b,4)+1),sin(pow(b,2))))/2+e);

cout<<"j="<<j;
}

С математикой плохо, наверное есть ошибки, подскажите кто разбирается.