円周率公式

 投稿者:しばっち  投稿日:2014年 1月 4日(土)19時30分19秒
 
LET S=0
LET A=1
FOR N=0 TO 5
   IF N>0 THEN
      LET A=-A*(6*N)*(6*N-1)*(6*N-2)*(6*N-3)*(6*N-4)*(6*N-5)
      LET A=A/(3*N)/(3*N-1)/(3*N-2)/N/N/N/640320/640320/640320
   END IF
   LET S=S+A*(13591409+545140134*N)
NEXT N
PRINT 426880*SQR(10005)/S

DIM U(8)
LET U(1)=4
LET U(2)=0
LET U(3)=0
LET U(4)=-2
LET U(5)=-1
LET U(6)=-1
LET U(7)=0
LET U(8)=0
LET S=0
FOR N=1 TO 50
   LET B=MOD(N,8)
   IF B=0 THEN LET B=8
   LET V=U(B)
   LET S=S+V/N/16^INT(N/8)
NEXT N
PRINT S

LET U(1)=1
LET U(2)=-3
LET U(3)=-2
LET U(4)=-3
LET U(5)=1
LET U(6)=0
LET A=1
LET S=0
FOR N=1 TO 100
   LET A=A/2
   LET B=MOD(N,6)
   IF B=0 THEN LET B=6
   LET V=U(B)
   LET S=S+A*V/N/N
NEXT N
PRINT SQR(ABS(36*S))

LET S=0
LET A=1
FOR N=1 TO 15
   LET S=SQR(2+S)
   LET A=A*2
NEXT N
PRINT A*2*SQR(2-S)

LET S=2*SQR(3)
LET A=3
LET B=1
FOR N=1 TO 50
   LET B=-B
   LET S=S+B*2/SQR(A)/(2*N+1)
   LET A=A*9
NEXT N
PRINT S

LET S=1
FOR N=100 TO 1 STEP-1
   LET S=(3+5*N+(N+1)*(2*N+1)/(3*N+4)/(3*N+5)/3*S)
NEXT N
PRINT 3+S/60

LET S=1
LET A=1
FOR N=1 TO 500
   LET A=A*N/(2*N+1)
   LET S=S+A
NEXT N
PRINT S*2

LET S=0
LET A=2
FOR N=0 TO 100
   IF N>0 THEN LET A=A*N*N*2/(2*N+1)/(2*N)
   LET S=S+A
NEXT N
PRINT S

LET S=0
LET A=1
FOR N=1 TO 10000
   LET S=S+A*1/(2*N)/(2*N+1)/(2*N+2)
   LET A=-A
NEXT N
PRINT 4*S+3

LET S=5*(SQR(5)-2)
LET A=.5
FOR N=0 TO 5
   LET X=5/S-1
   LET Y=(X-1)*(X-1)+7
   LET Z=(.5*X*(Y+SQR(Y*Y-4*X*X*X)))^(1/5)
   LET A=S*S*A-5^N*(.5*(S*S-5)+SQR(S*(S*S-2*S+5)))
   LET S=25/((Z+X/Z+1)*(Z+X/Z+1)*S)
NEXT N
PRINT 1/A

LET A=1/882
LET S=0
FOR N=0 TO 100
   IF N>0 THEN
      LET A=-A/882/882/32
      LET A=A*(2*N-1)*(4*N-1)*(4*N-3)/N/N/N
   END IF
   LET S=S+(1123+21460*N)*A
NEXT N
PRINT 4/S

LET B=(1+SQR(5))/2
LET A=1/B
LET S=1/B
FOR N=1 TO 20
   LET A=A*N*N/B/B/(2*N)/(2*N+1)
   LET S=S+A
NEXT N
PRINT S/2*5*SQR(B+2)

LET S=0
LET A=1/16
LET C=1
FOR N=0 TO 10
   LET C=COMB2(N)
   LET S=S+(42*N+5)*C*C*C*A
   LET A=A/4096
NEXT N
PRINT 1/S

LET S=3.14/2
FOR I=1 TO 10
   LET S=S+COS(S)
NEXT I
PRINT S*2

LET S=0
LET A=1
FOR N=0 TO 100
   LET S=S+COMB2(N)*A/(2*N+1)/(2*N+1)
   LET A=-A/16
NEXT N
PRINT SQR(10*S)

LET S=0
LET K=1
LET A=212175710912*SQR(61)+1657145277365
LET B=13773980892672*SQR(61)+107578229802750
LET C=(5280*(236674+30303*SQR(61)))^3
LET H=SQR(C)
FOR N=0 TO 3
   IF N>0 THEN
      LET K=-K*(6*N)*(6*N-1)*(6*N-2)*(6*N-3)*(6*N-4)*(6*N-5)
      LET K=K/(3*N)/(3*N-1)/(3*N-2)/N/N/N/C
   END IF
   LET S=S+K*(A+B*N)/H
NEXT N
PRINT 1/(12*S)

LET S=0
LET A=1
LET C=SQR(640320^3)
FOR N=0 TO 5
   IF N>0 THEN
      LET A=-A*(6*N)*(6*N-1)*(6*N-2)*(6*N-3)*(6*N-4)*(6*N-5)
      LET A=A/(3*N)/(3*N-1)/(3*N-2)/N/N/N/640320/640320/640320
   END IF
   LET S=S+A*(13591409+545140134*N)/C
NEXT N
PRINT 1/(12*S)

LET S=0
LET B=SQR(3)-SQR(2)
LET A=B*B
FOR N=0 TO 100
   LET S=S+R0(N)*(28*N+7-2*SQR(6))*A
   LET A=-A*B*B*B*B*B*B*B*B
NEXT N
PRINT 2/(S*3)

LET S=0
LET A=1
FOR N=0 TO 100
   LET S=S+R3(N)*(133*N+8)*A
   LET A=A*4*4*4/85/85/85
NEXT N
PRINT 85*SQR(85)/(18*S*SQR(3))

LET S=0
LET A=1
FOR N=0 TO 100
   LET S=S+R3(N)*(11*N+1)*A
   LET A=A*4/125
NEXT N
PRINT 5*SQR(5)/(2*S*SQR(3))

LET S=0
LET A=1
FOR N=0 TO 200
   LET B=-32/(4*N+1)-1/(4*N+3)+256/(10*N+1)-64/(10*N+3)-4/(10*N+5)-4/(10*N+7)+1/(10*N+9)
   LET S=S+B*A
   LET A=-A/1024
NEXT N
PRINT S/64

LET A=0
LET B=2
LET S=0
FOR N=0 TO 100
   LET AA=SQR(2+A)
   LET BB=SQR(2-A)
   LET S=S+B*(1-A/2)*2^N
   LET A=AA
   LET B=BB
NEXT N
PRINT S

LET S=0
LET A=1
FOR N=0 TO 200
   LET B=4/(6*N+1)+1/(6*N+3)+1/(6*N+5)
   LET S=S+B*A
   LET A=-A/8
NEXT N
PRINT S/SQR(2)

LET S=0
LET A=1
FOR N=0 TO 200
   LET B=2/(4*N+1)+2/(4*N+2)+1/(4*N+3)
   LET S=S+B*A
   LET A=-A/4
NEXT N
PRINT S

LET S=0
LET A=1
FOR N=0 TO 200
   LET B=4/(8*N+1)-2/(8*N+4)-1/(8*N+5)-1/(8*N+6)
   LET S=S+B*A
   LET A=A/16
NEXT N
PRINT S

LET A=1/3
LET B=2
FOR N=1 TO 100
   LET BB=4/(1+SQR((4-B)*(2+B)))
   LET AA=B*A+2^(N-1)/3*(1-B)
   LET B=BB
   LET A=AA
NEXT N
PRINT 1/A

LET S=0
LET A=1
FOR N=1 TO 100
   LET A=A*N*N/(2*N-1)/(2*N)
   LET S=S+A/N
NEXT N
PRINT 3*SQR(3)*S

LET S=0
FOR N=100 TO 1 STEP-1
   LET S=N*(2*N-1)/(3*N+1)/3/(3*N+2)*((5*N+3)+S)
NEXT N
PRINT 3+S

LET A=0
FOR N=10000 TO 1 STEP-1
   LET A=N*N/(2*N+1+A)
NEXT N
PRINT 4/(A+1)

LET S=0
LET A=1
FOR N=0 TO 100
   LET S=S+COMB2(N)/(2*N+1)*A
   LET A=A/16
NEXT N
PRINT 3*S

LET S=0
LET A=1
FOR N=1 TO 100
   LET A=A*3
   LET S=S+A/COMB2(N)
NEXT N
PRINT SQR(3)*(S-3)/4

LET S=0
FOR N=1 TO 100
   LET S=S+1/COMB2(N)/N
NEXT N
PRINT 3*SQR(3)*S

LET S=1
LET A=1
FOR N=1 TO 100
   LET A=A*(2*N-1)*(2*N-1)/(4*N)/(4*N+2)
   LET S=S+A
NEXT N
PRINT 3*S

LET S=1
LET A=1
FOR N=1 TO 100
   LET A=A*N*N/(2*N+1)/(2*N+2)
   LET S=S+A
NEXT N
PRINT SQR(9*S)

LET A=1
LET B=(SQR(3)-1)/2
LET S=1
FOR N=0 TO 10
   LET AA=(A+2*B)/3
   LET BB=(B*(A*A+A*B+B*B)/3)^(1/3)
   LET S=S+(3^(N+1)-3^N)*AA*AA
   LET A=AA
   LET B=BB
NEXT N
PRINT 3*A*A/(1-(S-3^N*A*A))

LET A=1
LET B=SQR(.5)
LET S=.25
FOR N=0 TO 10
   LET AA=(A+B)/2
   LET BB=SQR(A*B)
   LET S=S-2^N*(A-AA)*(A-AA)
   LET A=AA
   LET B=BB
NEXT N
PRINT(A+B)*(A+B)/(4*S)

LET A=SQR(.5)
LET B=A
FOR N=0 TO 100
   LET A=SQR(.5+.5*A)
   LET B=B*A
NEXT N
PRINT 2/B

LET S1=3*SQR(3)/4+2
LET A=1
LET S=0
FOR N=0 TO 100
   LET S=S+COMB2(N)/(N+1)/(2*N+5)*A
   LET A=A/16
NEXT N
PRINT S1-.75*S

LET S=0
LET A=1
FOR N=0 TO 100
   LET S=S+R4(N)*(20*N+3)*A
   LET A=-A/4
NEXT N
PRINT 8/S

LET S=0
LET A=1
FOR N=0 TO 100
   LET S=S+R6(N)*(33*N+4)*A
   LET A=A*4/125
NEXT N
PRINT(15*SQR(3))/(2*S)

LET S=0
LET A=1
FOR N=0 TO 100
   LET S=S+R4(N)*(26390*N+1103)*A
   LET A=A/99/99/99/99
NEXT N
PRINT 9801/(2*SQR(2)*S)

LET S=0
LET A=1
FOR N=0 TO 100
   LET S=S+R4(N)*(8*N+1)*A
   LET A=A/9
NEXT N
PRINT 2*SQR(3)/S

LET S=0
LET A=1
FOR N=0 TO 100
   LET S=S+(42*N+5)*A*R0(N)
   LET A=A/64
NEXT N
PRINT 16/S

LET S=0
LET S1=0
FOR N=1 TO 100
   LET C=COMB2(N)
   LET S1=S1+1/(N*N)
   LET S=S+1/(C*(2*N+1))*S1
NEXT N
PRINT(81*SQR(3)*S)^(1/3)

LET S=0
FOR N=1 TO 100
   LET C=COMB2(N)
   LET S=S+1/(C*(2*N+1)*(2*N+2)*(2*N+3))
NEXT N
LET A=1/18
LET B=2/SQR(3)
LET C=-(S+25/6)
PRINT(-B+SQR(B*B-4*A*C))/(2*A)

LET A=2
LET S=0
FOR N=2 TO 100
   LET A=A*2
   LET C=COMB2(N)
   LET S=S+A/C
NEXT N
PRINT 2*S

LET S=0
FOR N=1 TO 100
   LET C=COMB2(N)
   LET S=S+1/C
NEXT N
PRINT(S-1/3)*9*SQR(3)/2

LET A=1
LET S=0
FOR N=1 TO 100
   LET A=A*2
   LET C=COMB2(N)
   LET S=S+A/C/N/N
NEXT N
PRINT SQR(8*S)
END

EXTERNAL FUNCTION COMB2(N)
LET SS=1
IF N=0 THEN
   LET COMB2=1
   EXIT FUNCTION
END IF
FOR II=1 TO N
   LET SS=SS*(2*II-1)*(2*II)/II/II
NEXT II
LET COMB2=SS
END FUNCTION

EXTERNAL FUNCTION R6(N)
LET AA=1
LET BB=1
LET CC=1
FOR II=1 TO N
   LET AA=AA*(2*II-1)/(2*II)
   LET BB=BB*(3*II-2)/(3*II)
   LET CC=CC*(3*II-1)/(3*II)
NEXT II
LET R6=AA*BB*CC
END FUNCTION

EXTERNAL FUNCTION R4(N)
LET AA=1
LET BB=1
LET CC=1
FOR II=1 TO N
   LET AA=AA*(4*II-3)/(4*II)
   LET BB=BB*(2*II-1)/(2*II)
   LET CC=CC*(4*II-1)/(4*II)
NEXT II
LET R4=AA*BB*CC
END FUNCTION

EXTERNAL FUNCTION R3(N)
LET AA=1
LET BB=1
LET CC=1
FOR II=1 TO N
   LET AA=AA*(6*II-5)/(6*II)
   LET BB=BB*(2*II-1)/(2*II)
   LET CC=CC*(6*II-1)/(6*II)
NEXT II
LET R3=AA*BB*CC
END FUNCTION

EXTERNAL FUNCTION R0(N)
LET AA=1
FOR II=1 TO N
   LET AA=AA*(2*II-1)/(2*II)
NEXT II
LET R0=AA*AA*AA
END FUNCTION
 

円周率公式

 投稿者:しばっち  投稿日:2014年 1月 4日(土)19時31分5秒
  LET PP=1
LET X=1/8
FOR I=2 TO 50
   LET PP=PP+X/(I*2-1)
   LET X=X*(I*2-1)/(I*8)
NEXT I
PRINT PP*3

LET PP=1
FOR I=2 TO 5001
   IF I-INT(I/4)*4=1 OR I-INT(I/4)*4=2 THEN
      LET PP=PP+1/(I*2-1)
   ELSE
      LET PP=PP-1/(I*2-1)
   END IF
NEXT I
PRINT PP*2*SQR(2)

LET PP=0
FOR I=1 TO 5000
   LET J=I*8-7
   LET PP=PP+1/J/(J+2)/(J+4)/(J+6)
NEXT I
PRINT PP*96*(2+SQR(2))

LET PA=1
LET PB=1
LET A=1
LET B=1
FOR I=2 TO 20
   LET A=A*(I*2-2)/(I*2-1)/10
   LET B=B*(I*2-2)/(I*2-1)/50
   LET PA=PA+A
   LET PB=PB+B
NEXT I
PRINT PA*12/5+PB*14/25

LET PA=1
LET PB=1
LET A=1
LET B=1
FOR I=2 TO 20
   LET A=A*(I*2-2)/(I*2-1)/50
   LET B=B*(I*2-2)/(I*2-1)*.00144
   LET PA=PA+A
   LET PB=PB+B
NEXT I
PRINT PA*14/5+PB*7584/25000

LET PA=1
LET PB=1
LET A=1
LET B=1
FOR I=2 TO 30
   LET A=A*(I*4-4)/(I*2-1)/10
   LET B=B*(I*4-4)/(I*2-1)/100
   LET PA=PA+A
   LET PB=PB+B
NEXT I
PRINT PA*16/5-PB*14/25

LET PP=0
FOR I=1 TO 5000
   LET PP=PP+(I+1)/(I*2-1)/(I*2-1)/(I*2+1)/(I*2+1)
NEXT I
PRINT SQR((PP+3/8)*16)

LET PP=0
FOR I=1 TO 5000
   LET M=I*4-3
   LET PP=PP+I*I/M/(M+2)/(M+4)/(M+6)
NEXT I
PRINT 256*PP

LET PP=0
FOR I=1 TO 5000
   LET M=I*4-3
   LET PP=PP+1/M/(M+2)/(M+4)/(M+6)
NEXT I
PRINT(PP+1/18)*48

LET X=1
FOR I=1 TO 15
   LET X=2-2*SQR(1-X/4)
NEXT I
LET I=I-1
PRINT SQR(X)*6*(2^(I-1))

LET PP=0
FOR I=0 TO 5000
   LET X=1/(2*I+1)
   IF MOD(I,2)=0 THEN LET PP=PP+X ELSE LET PP=PP-X
NEXT I
PRINT PP*4

LET PP=1
FOR I=1 TO 10000
   LET PP=PP*4*I*I/(2*I-1)/(2*I+1)
NEXT I
PRINT PP*2

LET A=1
LET B=SQR(2)/2
LET C=1/4
LET X=1
FOR I=1 TO 5
   LET Y=A
   LET A=(A+B)/2
   LET B=SQR(B*Y)
   LET C=C-X*(A-Y)*(A-Y)
   LET X=X*2
NEXT I
PRINT (A+B)*(A+B)/4/C

LET PP=0
FOR I=1 TO 5001 STEP 4
   LET PP=PP+1/I/(I+2)
NEXT I
PRINT PP*8

LET X=SQR(2)/2
LET Y=1
FOR I=1 TO 100
   LET Y=Y*X
   LET X=SQR(.5+.5*X)
NEXT I
PRINT 2/Y

LET PP=0
FOR I=1 TO 10000
   LET PP=PP+1/I/I
NEXT I
PRINT SQR(PP*6)

LET PP=0
FOR I=20001 TO 1 STEP-2
   LET PP=I*I/(PP+2)
NEXT I
PRINT 4/(1+PP)

LET PP=0
FOR I=1 TO 10000 STEP 2
   LET PP=PP+2/(2*I-1)/(2*I+1)
NEXT I
PRINT PP*4

LET PP=0
LET A=1
FOR I=1 TO 100
   LET PP=PP+A
   LET A=A*I/(2*I+1)
NEXT I
PRINT PP*2

LET PP=0
FOR I=1 TO 10000
   IF MOD(I,2)=0 THEN
      LET PP=PP-1/I/I
   ELSE
      LET PP=PP+1/I/I
   END IF
NEXT I
PRINT SQR(12*PP)

LET PP=0
FOR I=1 TO 1000
   LET PP=PP+1/I/I/I/I
NEXT I
PRINT SQR(SQR(90*PP))

LET PP=1
FOR I=3 TO 10000 STEP 2
   LET PP=PP*(I*I-1)/I/I
NEXT I
PRINT 4*PP

LET PP=0
LET A=1/2
LET B=1/3
FOR I=1 TO 30
   LET PP=PP+A/(2*I-1)+B/(2*I-1)
   LET A=-A/4
   LET B=-B/9
NEXT I
PRINT 4*PP

LET PP=5/6
LET A=1/2
FOR I=2 TO 10000
   LET A=A*(2*I-3)/(2*I)
   LET PP=PP-A/(2*I+1)
NEXT I
PRINT 4*PP

LET PP=0
FOR I=1 TO 10000
   IF MOD(I,2)=0 THEN
      LET PP=PP-1/(2*I)/(2*I+1)/(2*I+2)
   ELSE
      LET PP=PP+1/(2*I)/(2*I+1)/(2*I+2)
   END IF
NEXT I
PRINT 4*PP+3

LET A=1/6
LET PP=1-A
FOR I=1 TO 50
   LET A=A*I/(2*I+3)
   LET PP=PP-A
NEXT I
PRINT 4*PP

LET PP=1
FOR I=1 TO 10000
   LET PP=PP*SQR(1+1/I)/(1+1/(2*I))
NEXT I
PRINT PP*PP*4

LET A=1
LET PP=1
FOR I=1 TO 1000
   LET A=A*(2*I-1)*(2*I-1)/(4*I)/(2+4*I)
   LET PP=PP+A
NEXT I
PRINT 3*PP

LET A=1
LET B=1/2
LET PP=1/2
FOR I=1 TO 40
   LET A=A*(2*I-1)/(2*I)
   LET B=B/4
   LET PP=PP+A*B/(2*I+1)
NEXT I
PRINT PP*6

LET A=1/16
LET PP=1/2-1/48
FOR I=2 TO 100
   LET A=A*(2*I-3)/(8*I)
   LET PP=PP-A/(2*I+1)
NEXT I
PRINT(PP-SQR(3)/8)*12

LET A=SQR(2)
LET B=0
LET PP=2+SQR(2)
FOR I=1 TO 5
   LET B=(1+B)/(A+B)
   LET A=SQR(A)
   LET B=B*A
   LET A=(A+1/A)/2
   LET PP=PP*B*(A+1)/(1+B)
NEXT I
PRINT PP

LET C=0
LET S=1
FOR K=2 TO 30
   LET C=SQR((1+C)/2)
   LET S=S/C/2
NEXT K
PRINT 2^(K-1)*S*C

LET A=1
LET B=SQR(2)
DO WHILE ABS(A-B)>.000000000000001
   LET C=A
   LET A=(A+B)/2
   LET B=SQR(B*C)
LOOP
LET S=1
LET E=1
FOR I=1 TO 30000
   LET E=E*(2*I-1)/(2*I)
   LET S=S+E/(4*I+1)
NEXT I
PRINT S*A*2

LET A=2*SQR(2)/9801
LET PP=0
FOR I=0 TO 2
   IF I>0 THEN
      LET A=A*(4*I)*(4*I-1)*(4*I-2)*(4*I-3)
      LET A=A/I/I/I/I/396/396/396/396
   END IF
   LET PP=PP+A*(1103+26390*I)
NEXT I
PRINT 1/PP

LET A=1
LET D=.5
FOR I=1 TO 19
   LET A=SQR(2-SQR(4-A*A))
   LET D=D*2
NEXT I
PRINT D*A*6

LET A=6-4*SQR(2)
LET Y=SQR(2)-1
FOR I=0 TO 1
   LET S=SQR(SQR((1-Y^4)))
   LET Y=(1-S)/(1+S)
   LET A=(1+Y)^4*A-(2^(2*I+3))*(1+Y+Y*Y)*Y
NEXT I
PRINT 1/A

LET PP=0
FOR I=1 TO 3000
   LET PP=PP+1/(4*I*I-1)/(4*I*I-1)
NEXT I
PRINT SQR(16*PP+8)

LET PP=0
LET A=1
FOR I=1 TO 30
   LET PP=PP+A/(2*I-1)
   LET A=-A/3
NEXT I
PRINT 2*SQR(3)*PP

LET A=1
LET S=1
FOR I=1 TO 30
   LET A=A*I*I/(2*I+1)/(2*I+2)
   LET S=S+A
NEXT I
PRINT 6*SQR(S/4)

LET PP=3/4*SQR(3)+2-.15
LET A=1.5
LET B=1
LET C=.25
FOR I=1 TO 20
   LET A=A*(2*I-1)/2
   LET B=B/(I+1)
   LET C=C/4
   LET PP=PP-A*B*C/(I+2.5)
NEXT I
PRINT PP

LET PP=0
FOR I=1 TO 10000
   IF MOD(I,2)=1 THEN
      LET PP=PP+1/(I*2-1)/(I*2-1)/(I*2-1)
   ELSE
      LET PP=PP-1/(I*2-1)/(I*2-1)/(I*2-1)
   END IF
NEXT I
PRINT(PP*32)^(1/3)

LET A=SQR(2)
LET B=0
LET PP=2+SQR(2)
FOR I=1 TO 10
   LET AA=(SQR(A)+1/SQR(A))/2
   LET BB=SQR(A)*(B+1)/(A+B)
   LET A=AA
   LET B=BB
   LET PP=PP*B*(1+A)/(1+B)
NEXT I
PRINT PP

LET P=6541681608/(640320^(3/2))
LET A=1
LET PP=0
FOR I=0 TO 3
   IF I>0 THEN
      LET A=A*(6*I)*(6*I-1)*(6*I-2)*(6*I-3)*(6*I-4)*(6*I-5)
      LET A=A/(3*I)/(3*I-1)/(3*I-2)
      LET A=A/I/I/I
   END IF
   LET PP=PP+(13591409/545140134+I)*A*((-1)^I)/(640320^(3*I))
NEXT I
PRINT 1/(P*PP)

LET X=SQR(2)
LET PP=2+SQR(2)
LET Y=2^(1/4)
FOR I=1 TO 3
   LET X=(SQR(X)+1/SQR(X))/2
   LET PP=PP*(X+1)/(Y+1)
   LET Y=SQR(Y)*(SQR(X)+1/SQR(X))/(Y+1)
NEXT I
PRINT PP

LET PP=0
LET B=1
FOR I=0 TO 10
   LET A=(4/(8*I+1)-2/(8*I+4)-1/(8*I+5)-1/(8*I+6))
   LET PP=PP+A*B
   LET B=B/16
NEXT I
PRINT PP

LET A=5/16
LET S=A
FOR N=1 TO 10
   LET A=A*(2*N-1)^3*(42*N+5)/512/N^3/(42*N-37)
   LET S=S+A
NEXT N
PRINT 1/S

LET A=1/3
LET R=(SQR(3)-1)/2
LET S=(1-R^3)^(1/3)
FOR N=0 TO 1
   LET T=1+2*R
   LET U=(9*R*(1+R+R^2))^(1/3)
   LET V=T^2+T*U+U^2
   LET M=27*(1+S+S^2)/V
   LET A=M*A+3^(2*N-1)*(1-M)
   LET S=(1-R)^3/(T+2*U)/V
   LET R=(1-S^3)^(1/3)
NEXT  N
PRINT 1/A

LET A=1/3
LET S=(SQR(3)-1)/2
FOR N=0 TO 5
   LET R=3/(1+2*(1-S^3)^(1/3))
   LET S=(R-1)/2
   LET A=R*R*A-3^N*(R^2-1)
NEXT N
PRINT 1/A

LET S=0
FOR I=1000 TO 1 STEP-1
   LET S=I/(2*I+1)*(1+S)
NEXT I
PRINT (S+1)*2

LET S=0
FOR I=1000 TO 1 STEP-1
   LET S=1/(2*I)+I/(2*I+1)*S
NEXT I
PRINT SQR(16*S)

LET S=0
FOR I=1000 TO 1 STEP-1
   LET S=(2*I-3)^2/(2*I)^2*(1+S)
NEXT I
PRINT 4/(1+S)

LET S=0
FOR I=3000 TO 2 STEP-1
   LET S=(2*I-3)^2/(2+S)
NEXT I
PRINT 4/(1+S)

LET S=0
FOR I=1000 TO 1 STEP-1
   LET S=I^2/(2*I+1)/(2*I+2)*(1+S)
NEXT I
PRINT SQR(9*(S+1))

LET S=0
FOR I=10000 TO 2 STEP-1
   LET S=I*(I-1)/(1+S)
NEXT I
PRINT 2+2/(1+S)

LET S=1
FOR N=5000 TO 1 STEP-1
   LET S=1+S*(2*N-3)^2/(2*N)^2
NEXT N
LET S=S/4
PRINT 1/S

LET A=.5
LET S=0
FOR N=0 TO 100
   IF N>0 THEN LET A=A*(2*N-1)/(2*N)/4
   LET S=S+A/(2*N+1)
NEXT N
PRINT 6*S

LET A=2*SQR(3)
LET B=3
FOR N=1 TO 100
   LET A=2*A*B/(A+B)
   LET B=SQR(A*B)
NEXT N
PRINT A
END
 

円周率公式

 投稿者:しばっち  投稿日:2014年 1月 4日(土)19時31分40秒
  OPTION ARITHMETIC DECIMAL_HIGH
DECLARE EXTERNAL FUNCTION COMB
DECLARE EXTERNAL FUNCTION PP
LET S=0
FOR N=1 TO 500
   LET S=S+3*PP(N)/COMB(N)/2^(N-1)
NEXT N
PRINT(S-20379280)/740025

LET A=63365028312971999585426220+28337702140800842046825600*SQR(5)+384*SQR(5)*SQR(10891728551171178200467436212395209160385656017+4870929086578810225077338534541688721351255040*SQR(5))
LET B=7849910453496627210289749000+3510586678260932028965606400*SQR(5)+2515968*SQR(3110)*SQR(6260208323789001636993322654444020882161+2799650273060444296577206890718825190235*SQR(5))
LET C=-214772995063512240-96049403338648032*SQR(5)-1296*SQR(5)*SQR(10985234579463550323713318473+4912746253692362754607395912*SQR(5))
LET CC=SQR(-C*C*C)
LET X=1
LET S=A
FOR N=1 TO 500
   LET X=X*(6*N)*(6*N-1)*(6*N-2)*(6*N-3)*(6*N-4)*(6*N-5)
   LET X=X/(3*N)/(3*N-1)/(3*N-2)/N/N/N/C/C/C
   LET S=S+X*(A+B*N)
NEXT N
PRINT CC/S

LET A=21242668516504965+15020834958518500*SQR(2)+2*SQR(5)*SQR(45125096427586568251645610141659+31908261685643312902173585434250*SQR(2))
LET B=1839779353703421900+1300920456890691000*SQR(2)+24337404*SQR(10)*SQR(1142912476713024496667+808161162586491705750*SQR(2))
LET J=(71864175655+22725423252*SQR(10)+2808*SQR(5)*SQR(261993316778681+82849561276216*SQR(10)))^3
LET AA=1
LET S=0
FOR N=0 TO 500
   IF N>0 THEN
      LET AA=AA*(6*N-1)*(6*N-3)*(6*N-5)/(6*N)/(6*N)/(6*N)
      LET AA=AA/N/N/N/J
   END IF
   LET S=S+AA*(A+B*N)
NEXT N
PRINT SQR(3*J)/S
END

EXTERNAL FUNCTION COMB(N)
OPTION ARITHMETIC DECIMAL_HIGH
LET SS=1
FOR II=1 TO N
   LET SS=SS*(7*II-6)*(7*II-5)*(7*II-4)*(7*II-3)*(7*II-2)*(7*II-1)*(7*II)/(5*II)/(5*II-1)/(5*II-2)/(5*II-3)/(5*II-4)/(2*II)/(2*II-1)
NEXT II
LET COMB=SS
END FUNCTION

EXTERNAL FUNCTION PP(N)
OPTION ARITHMETIC DECIMAL_HIGH
LET PP=-885673181*N*N*N*N*N+3125347237*N*N*N*N-2942969225*N*N*N+1031962795*N*N-196882274*N+10996648
END FUNCTION
 

戻る