REM Maclaurin(円周率π Machinの公式)
input N
LET n1=(0.693147+(N+1)*2.302585)/(2*1.609438)+1/2
LET n2=(0.693147+(N+1)*2.302585)/(2*5.476464)+1/2
LET S1=0
LET x1=1/5
LET a1=x1
LET S1=a1
FOR k=2 TO n1
   LET a1=-(2*k-3)/(2*k-1)*x1^2*a1
   LET S1=S1+a1
NEXT k
LET S2=0
LET x2=1/239
LET a2=x2
LET S2=a2
FOR k=2 TO n2
   LET a2=-(2*k-3)/(2*k-1)*x2^2*a2
   LET S2=S2+a2
NEXT k

PRINT (4*S1-S2)*4
print pi
END