C语言：公式法计算圆周率

#include <stdio.h>
int main()
{
    double a = 1, b = 1, i, j, s = 0, s0 = 0, pi, x;
    scanf("%lf", &x);
    for (i = 1,j = 3; i < 2000; i++,j += 2)
    {
        a *= (i + 1);
        b *= j;
        s0 = a / b;
        s += a / b;
        if (s0 < x)
            break;
    }
    pi = 2 * s;
    printf("%.11lf", pi);
    return 0;
}