Это же целочисленное интегрирование. К тому же самое дубовое - методом прямоугольников. Погрешность, однако. Чем больше n, тем ближе к аналитическому результату. Возьми, например n = 1000, получишь 22.2484 n = 10000 - 22.1442 n = 100000 - 22.1418 и т. д.