（１）　平方数の逆数の和を用いる方法　（ ANSI/ISO/JIS規格Ｃ言語 ）

#include <stdio.h>
#include <math.h>

main()
{
   int i,n;
   double x,y;
   x=0.0;
   n=1;
   for(i=1;i<=10000;i++){
       x+=1.0/(n*n);
       n++;
                        }
   y=sqrt(6.0*x);
   printf("円周率は %1.6f です。\n",y);
   return 0;
}




（２）　マチンの公式を用いる方法　（ JIS Full ＢＡＳＩＣ言語 ）

OPTION ARITHMETIC decimal_HIGH
LET X=1/5
LET F=-X*X
LET T=X
LET S=0
FOR J=1 TO 1000
   LET A=T/(2*J-1)
   LET S=S+A
   LET T=F*T
NEXT J
LET Y=1/239
LET B=-Y*Y
LET D=Y
LET E=0
FOR K=1 TO 1000
   LET C=D/(2*K-1)
   LET E=E+C
   LET D=B*D
NEXT K
PRINT "円周率は次のとおりです。\n"
PRINT 16*S-4*E
END