5-15 计算圆周率 (15分)

#include <stdio.h>
#include <iostream>
using namespace std;

int main()
{
   double n;
   double i=1.0,j=3.0,sum=1,x=1.0,y=1.0,a=1;
   cin>>n;
   while(a>n)
   {
       x*=i;
       y*=j;
       i++;
       j=2*i+1;
       a=x/y;
       sum+=a;
   }
   sum=2*sum;
   printf("%.6f",sum);
}