Программа

' Решение дифура y'=f(x,y)
'  при известных начальных условиях
'  x0, y0=f(x0) методом Рунге-Кутта
DEF FNF! (X!) 'Вычисление производной
    FNF! = 3 * X! * X!
END DEF
DEF FNDifRunKut! (X, Y, Break%, Xresult)
                                    'Вычисление значения исходной функции
    S! = (Xresult - X) / Break%     'Определение шага разбиения
    FOR Meter% = 0 TO Break% - 1 'Цикл вычисления первообразной в точке Xresul
        k1 = S! * (FNF!(X + Meter% * S!))
        k2 = S! * (FNF!(X + Meter% * S! + S! / 2))
        k3 = S! * (FNF!(X + Meter% * S! + S! / 2))
        k4 = S! * (FNF!(X + Meter% * S! + S!))
        Y = Y + (k1 + 2 * k2 + 2 * k3 + k4) / 6
     NEXT Meter%
    FNDifRunKut! = Y                    'Определение значения функции
END DEF

PRINT (FNDifRunKut!(0, 0, 10, 10))  'Вычисление первообразной
END