Такс.... значение функции на заданном промежутке не превышает 1.0457.
Если кто-то из математиков поможет найти абсолютный максимум (он же единственный локальный) было бы здорово.
Но будем исходить из этого.
Функция строго положительна на заданном промежутке, отсюда можем сказать, что ошибка для любого из прямоугольников не будет превышать величины h*1.0457.
Отсюда получаем h*1.0457<=t*(1/1.0457).
h = t/(1.0457^2).
Таким образом мы нашли максимально допустимый шаг.
Находим целое количество шагов и округляем его вверх.
m = -Int(-1/h) '-Int(-x) - это округление вверх в VB.
Для такого количества шагов рассчитываем новую величину шага
h=1/m
Вот с этим шагом и рассчитываем. Дальше разберетесь?

Осторожно, у меня тут в многих местах единички только потому, что у Вас промежуток равен 1.