17 lines
418 B
Plaintext
17 lines
418 B
Plaintext
LET y = 1
|
|
FOR i = 0 TO 100
|
|
LET t = i / 10
|
|
|
|
IF t = INT(t) THEN
|
|
LET actual = ((t ^ 2 + 4) ^ 2) / 16
|
|
PRINT "y("; STR$(t); ") ="; y ; TAB(20); "Error = "; actual - y
|
|
END IF
|
|
|
|
LET k1 = t * SQR(y)
|
|
LET k2 = (t + 0.05) * SQR(y + 0.05 * k1)
|
|
LET k3 = (t + 0.05) * SQR(y + 0.05 * k2)
|
|
LET k4 = (t + 0.10) * SQR(y + 0.10 * k3)
|
|
LET Y = Y + 0.1 * (k1 + 2 * (k2 + k3) + k4) / 6
|
|
NEXT i
|
|
END
|