°°ħ°°‰	DOUBLE PRECISION PILOW,SLNGTH,PIUP,FACT
°°²°°‰	DOUBLE PRECISION K,SIDES,SUM,SSQ,TERM,TEMP
°°³°°‰	BYTE ITU(80)
°°´°°‰
°°µ°°‰
°°ĥ°°‰	CALL OPEN (6,'RFILE   DAT',0)
°°·°°‰
°°¸°°‰	WRITE(6,198,REC=1)
°°ı°°‰198	FORMAT(10X,' BOUNDS ON PI - DOUBLE PRECISION BINOMIAL THEOREM',
°ħ°°°‰     1  ' VERSION.')
°ħħ°°‰	WRITE(6,199,REC=2)
°ħ²°°‰199	FORMAT('   N     SIDES    SIDE LENGTH   PI - LOWER BOUND',
°ħ³°°‰     1  '   PI - UPPER BOUND')
°ħ´°°‰	SIDES=4.0
°ħµ°°‰	SUM=2.0
°ħĥ°°‰	N=3
°ħ·°°‰1	SIDES=2.0*SIDES
°ħ¸°°‰	SSQ=SUM
°ħı°°‰	SUM=0.0
°²°°°‰	TERM=.25*SSQ
°²ħ°°‰	K=1
°²²°°‰2	TEMP=TERM+SUM
°²³°°‰	IF(TEMP.LE.SUM) GOTO 4
°²´°°‰	SUM=TEMP
°²µ°°‰	FACT=(2.0*K-1.00)/(K+1.0)
°²ĥ°°‰	TERM=FACT*SSQ*TERM/8.0
°²·°°‰	K=K+1.0
°²¸°°‰	GOTO 2
°²ı°°‰4	SLNGTH=DSQRT(SUM)
°³°°°‰	PILOW=0.5*SIDES*SLNGTH
°³ħ°°‰	PIUP=SIDES*SLNGTH/(2.0-SLNGTH)
°³²°°‰	WRITE(6,200,REC=N) N,SIDES,SLNGTH,PILOW,PIUP
°³³°°‰	IF(N.EQ.20) GOTO 8
°³´°°‰	N=N+1
°³µ°°‰	GOTO 1
°³ĥ°°‰200	FORMAT(1X,I3,F9.0,F15.6,2F19.12)
°³·°°‰
°³¸°°‰8	READ(6,300,REC=1) ITU
°³ı°°‰300	FORMAT(80A1)
°´°°°‰	WRITE(1,300) ITU
°´ħ°°‰	WRITE(1,302)
°´²°°‰	READ(6,300,REC=2) ITU
°´³°°‰	WRITE(1,300) ITU
°´´°°‰302	FORMAT(1X)
°´µ°°‰	DO 9 I = 3,20
°´ĥ°°‰	READ(6,300,REC=I) ITU
°´·°°‰9	WRITE(1,300) ITU
°´¸°°‰	WRITE(1,302)
°´ı°°‰	STOP
°µ°°°‰	END