Qbasic "Численное интегрирование"

[PaScalpel]

Прошу проверить программы и написать уравнение F= так как 0 в бейсике
Выдает ошибку в обоих программах "индекс вне диапазона"

Уравнение - x+sqrt(y^2+1.5x^2)
Ввод. числа
x0= 0
y0 = 0.5
B = 1
N= 10

10 PRINT "Решение диф. уравнения"
20 PRINT "усовершенствованным методом Эйлера-Коши"
25 REM
30 INPUT "Ввести нач. значение x0"; x0
40 INPUT "Ввести нач. значение Y0="; Y0
50 INPUT "Ввести конец отрезка интегрирования B="; B
60 INPUT "Ввести число деления отрезка интегрирования N="; N
70 H = (B - x0) / N: H2 = H / 2
80 x = x0: Z = Y0: i = 0: y = Z
90 IF i = N THEN 60
100 GOSUB 160: V = F
110 y = Z + H * F: x = x + H: GOSUB 160
120 y = Z * H2 * (V + F)
130 PRINT "Y(" x ") ="; y
140 Z = y: i = i + l
150 GOTO 90
160 F = F(x, y): F =  x + sqrt(y * y + 1.5 * x * x)
170 RETURN

10 PRINT "Решение задачи Коши"
20 PRINT "методом Рунге-Кутта 4-го порядка точности"
30 PRINT "						"
40 INPUT "Ввести x0="; x0
50 INPUT "Ввести x0="; Y0
60 INPUT "Ввести конец отрезка интегрирования b="; b
70 INPUT "Ввести число деления отрезка N="; N
80 H = (b - x0) / N: H2 = H / 2
90 x = x0: Z = Y0: y = Z: i = 0
100 IF i > N THEN 200
110 GOSUB 210: K1 = H * F
120 x = x + H2: y = Z + K1 / 2
130 GOSUB 210: K2 = H * F
140 y = Z + K2 / 2: GOSUB 210: K3 = H * F
150 x = x + H2: y = Z + K3: GOSUB 210: K4 = H * F
160 y = Z + (K1 + 2 * (K2 + K3) + K4) / 6
170 PRINT "Y(" x ")="; y
180 Z = y: i = i + 1
190 GOTO 100
200 GOTO 70
210 F = x + sqrt(y * y + 1.5 * x * x): RETURN
220 END