14 lines
387 B
Plaintext
14 lines
387 B
Plaintext
decimals(8)
|
|
y = 1.0
|
|
for i = 0 to 100
|
|
t = i / 10
|
|
if t = floor(t)
|
|
actual = (pow((pow(t,2) + 4),2)) / 16
|
|
see "y(" + t + ") = " + y + " error = " + (actual - y) + nl ok
|
|
k1 = t * sqrt(y)
|
|
k2 = (t + 0.05) * sqrt(y + 0.05 * k1)
|
|
k3 = (t + 0.05) * sqrt(y + 0.05 * k2)
|
|
k4 = (t + 0.10) * sqrt(y + 0.10 * k3)
|
|
y += 0.1 * (k1 + 2 * (k2 + k3) + k4) / 6
|
|
next
|