Trailing-Edge
-
PDP-10 Archives
-
decuslib20-01
-
decus/20-0020/romint.clc
There are 2 other files named romint.clc in the archive. Click here to see a list.
10000' NAME--ROMINT
10010'
10020' DESCRIPTION--SUBROUTINE FOR PERFORMING ROMBERG'S EQUATION.
10030'
10040' SOURCE--UNKNOWN
10050'
10060' INSTRUCTIONS--THE SUBROUTINE ASSUMES THAT L IS EQUAL
10070' TO THE LEFT POINT END POINT, THAT R IS EQUAL TO THE
10080' RIGHT END POINT, THAT O IS THE DESIRED ORDER OF THE INTEGRATION,
10090' AND THET FNC IS THE FUNCTION TO BE INTEGRATED. THE ANSWER IS
10100' V. THE ROUTINE USES THE LIST Z, AND THE VARIABLES
10110' Z5 THROUGH Z9. THE SUBROUTINE PROPER STARTS IN LINE 10170
10120' AND CONTINUES THROUGH LINE 10350.
10130'
10140'
10150' * * * * * * SUBROUTINE * * * * * * * * * * *
10160'
10170 LET Z7 = (R - L)/2
10180 LET Z8 = 1
10190 LET Z(1) = (FNC(L) + FNC(R))/2
10200 FOR Z5 = 1 TO O
10210 LET Z6 = 0
10220 FOR Z9 = 1 TO 2*Z8 - 1 STEP 2
10230 LET Z6 = Z6 + FNC(L + Z9*Z7)
10240 NEXT Z9
10250 LET Z(Z5 + 1) = ( Z6/Z8 + Z(Z5) )/2
10260 LET Z6 = 1
10270 FOR Z9 = Z5 TO 1 STEP -1
10280 LET Z6 = 4*Z6
10290 LET Z(Z9) = Z(Z9 + 1) + (Z(Z9 + 1) - Z(Z9))/(Z6 - 1)
10300 NEXT Z9
10310 LET Z8 = Z8*2
10320 LET Z7 = Z7/2
10330 NEXT Z5
10340 LET V = Z(1)*(R - L)
10350 RETURN