#/***********************************************************
#    pi1.c -- ±ß¼þΨ
#***********************************************************/

def pi  # Machin¤Î¸ø¼°
    p = 0;  k = 1;  t = 16.0 / 5.0
    begin 
        last = p;  p += t / k;  t /= -5.0*5.0;  k += 2
    end while (p != last)
    k = 1;  t = 4.0 / 239.0
    begin 
        last = p;  p -= t / k;  t /= -239.0*239.0;  k += 2
    end while (p != last)
    return p
end

LDBL_DIG = 18    # float.h

printf("pi = %.*g\n", LDBL_DIG, pi)

exit 0