>b's weblog

double leibnitz(double n) {
    double p=1.0;

    for (int i=3, N=2*n+1; i<=N; i+=2)
        p += ((i&2) ? -1.0 : 1.0) / i;

    return 4*p;
}