1 DIM X(20),F(10),K(4,10)
3 GOTO 100
5 REM --SUBROUTINE 
99 RETURN 
210 PRINT "T              ",  
220 READ N,T 
230 FOR I=1 TO N  
240 READ X(I)
250 PRINT "X";I;,
260 NEXT I  
270 READ D,B 
280 A=T
285 FOR J=A TO B STEP D 
286 PRINT 
288 PRINT T, 
290 FOR I=1 TO N  
292 PRINT X(I), 
294 NEXT I  
300 FOR I=1 TO N  
310 X(N+I)=X(I) 
320 NEXT I  
330 GOSUB 5
340 FOR I=1 TO N  
350 K(1,I)=F(I)
360 NEXT I  
370 D1=D/2 
380 T=T+D1  
390 FOR L=1 TO 3 
400 T=T+INT(L/3)*D1 
410 FOR I=1 TO N  
420 X(I)=X(N+I)+D1*K(L,I)*INT(L/3+1)  
430 NEXT I  
440 GOSUB 5
450 FOR I=1 TO N  
460 K(L+1,I)=F(I) 
470 NEXT I  
480 NEXT L  
610 FOR I=1 TO N  
620 X(I)=X(N+I)+D/6*(K(1,I)+2*K(2,I)+2*K(3,I)+K(4,I))
630 NEXT I  
640 NEXT J  
999 END 

