Вторая производная функции EXP(x) - QBasic

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

Все знаки (10 штук) верные

REM y''(1) = 2.718281828
DEFDBL A-Z
 
CLS
dx = .000005
x = 1
y1 = (EXP(x + 2 * dx) - EXP(x - dx)) / (3 * dx)
y2 = (EXP(x + dx) - EXP(x - 2 * dx)) / (3 * dx)
y = (y1 - y2) / dx
 
PRINT USING "Y'' = #.#########"; y
END

Решение задачи: «Вторая производная функции EXP(x)»

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

REM y''(1) = 2.718281828
DEFDBL A-Z
 
CLS
dx = .000005
x = 1
y1 = (EXP(x + dx) - EXP(x)) / (dx)
y2 = (EXP(x) - EXP(x - dx)) / (dx)
y = (y1 - y2) / (dx)
 
PRINT USING "Y'' = #.#########"; y
END

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

Вторая производная функции EXP(x) находится по формуле y''(x) = EXP(x) EXP(x) - 2 EXP(x) + 1
Задана точность dx = 0.000005 для вычисления производной
Инициализированы переменные x = 1 и dx
Вычислена приближенно первая производная функции EXP(x) по формуле y1 = (EXP(x + dx) - EXP(x)) / (dx)
Вычислена приближенно вторая производная функции EXP(x) по формуле y2 = (EXP(x) - EXP(x - dx)) / (dx)
Вычислена приближенно вторая производная функции EXP(x) по формуле y = (y1 - y2) / (dx)
Выводится результат вычисления y на экран с помощью команды PRINT USING Y'' = #.#########; y
Код выполняется до команды END