#reg 0 n
#reg c integral
#reg 1 xi
#reg 2 h
#reg 3 k
#reg 4 temp1
#reg 5 temp2


DO
        =n				; n = input value
        2 * 1/x @B @A - *  =h		; h=(b-a)/2n
        $0.6 SQR *  =k			; k=h*srq(3/5)
        0 =integral			; integral = 0
        
        @a @h 
        
        For n
		+ =xi			; center xi=a+h
        	Gosub f(x) 8 * =temp1	; 8f(xi)
        	@xi @k -	
        	Gosub f(x) 5 * =temp2	; 5f(xi-k)
        	@xi @k +
        	Gosub f(x) 5 * 		; 5f(xi+k)
		@temp1 + @temp2 +
		sum integral		; integral = integral+5f(xi-k)+8f(xi)+5f(xi+k)
        	@xi @h 2 *		; xi=xi+2h, next center
        Next
        
        @integral 9 / @h * =integral	; integral = integral*h/9
        
        R/S				; display and wait for new n
LOOP					; go on with new value

f(x):
RTN

	
	