11 lines
318 B
Plaintext
11 lines
318 B
Plaintext
100 PROGRAM "Runge.bas"
|
|
110 LET Y=1
|
|
120 FOR T=0 TO 10 STEP .1
|
|
130 IF T=INT(T) THEN PRINT "y(";STR$(T);") =";Y;TAB(21);"Error =";((T^2+4)^2)/16-Y
|
|
140 LET K1=T*SQR(Y)
|
|
150 LET K2=(T+.05)*SQR(Y+.05*K1)
|
|
160 LET K3=(T+.05)*SQR(Y+.05*K2)
|
|
170 LET K4=(T+.1)*SQR(Y+.1*K3)
|
|
180 LET Y=Y+.1*(K1+2*(K2+K3)+K4)/6
|
|
190 NEXT
|