16 lines
428 B
Plaintext
16 lines
428 B
Plaintext
y = 1.0
|
|
FOR i% = 0 TO 100
|
|
t = i% / 10
|
|
|
|
IF t = INT(t) THEN
|
|
actual = ((t^2 + 4)^2) / 16
|
|
PRINT "y("; t ") = "; y TAB(20) "Error = "; actual - y
|
|
ENDIF
|
|
|
|
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 += 0.1 * (k1 + 2 * (k2 + k3) + k4) / 6
|
|
NEXT i%
|