Расчет n значений функции с равномерным шагом - QBasic

Формулировка задачи:

Не знаю как написать код программы и решить задачу: Для заданной подынтегральной функции Y=f(x) на интервале х{[а,b]} рассчитать n значений функции с равномерным шагом изменения аргумента х, значения аргумента и функции предоставить в виде таблицы. Исходные данные: Y=(интеграл)от b до а √(2*x+1)*dx здесь n=40, a=0, b=1.

Решение задачи: «Расчет n значений функции с равномерным шагом»

Листинг программы

DECLARE FUNCTION f! (x!)
 
CLS
a = 0
b = 1
n = 40
h = (b - a) / (n - 1)
 
FOR i = 1 TO n
   x = a + (i - 1) * h
   y = f(x)
   PRINT USING "x = #.##   y = #.### "; x; y
 
   IF i MOD 20 = 0 THEN
      PRINT "PRESS ANY KEY"
      DO: LOOP UNTIL INKEY$ <> ""
      CLS
   END IF
NEXT
END
 
FUNCTION f (x)
   f = SQR(2 * x + 1)
END FUNCTION

Объяснение кода листинга программы

Объявляется функция с именем f(x)
Задаются начальные значения переменных a=0, b=1, n=40
Вычисляется шаг h=(b-a)/(n-1)
В цикле for i от 1 до n происходит вычисление x=a+(i-1)*h и y=f(x)
Выводится на экран значение x и y с интервалом в 20 строк
Если значение i кратно 20, то выводится сообщение PRESS ANY KEY и происходит проверка на нажатие клавиши
При нажатии клавиши, цикл прерывается и программа продолжает работу
Код функции f(x) вычисляет значение f(x)=SQR(2*x+1)