17 lines
400 B
Plaintext
17 lines
400 B
Plaintext
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
|