y = 1
for i = 0 to 100
    t = i / 10

    if t = int(t) then
	actual = ((t ^ 2 + 4) ^ 2) / 16
	print "y("; int(t); ") = "; left(string(y), 13), "Error = "; left(string(actual - y), 13)
    end if

    k1 = t * sqr(y)
    k2 = (t + 0.05) * sqr(y + 0.05 * k1)
    k3 = (t + 0.05) * sqr(y + 0.05 * k2)
    k4 = (t + 0.10) * sqr(y + 0.10 * k3)
    y = y + 0.1 * (k1 + 2 * (k2 + k3) + k4) / 6
next i
end