2013 / piprimer / piprim

#include <iostream>
#include <fstream>
#include <iomanip>
#include <cmath>
using namespace std;
ifstream fin;
ofstream fout;

int main ()
{
fout.open("rezpi2.txt");
if (fout.fail())
{cout<<"error open  rezpi1.txt \n"; return(1);}
fout.precision(20);
fout<<"pi="<<4*atan(1.0)<<endl;
long long n=3;
double a=1,a1=1;
double p=3,p1=3,c=sqrt(2.0);
fout<< setw(15)<<"   n    "<<setw(30)<<"pi(плохая формула)"<<
         setw(30)<<"pi(хорошая формула)"<<endl;
fout << setw(15)<<n<<setw(30)<<p<<setw(30)<<p1<<endl;
for(int i=0;i<40;i++)
{n=2*n;
  a=c*sqrt(1-sqrt(1-a*a/4));           // {плохая формула} 
  a1=a1/(c*sqrt(1+sqrt(1-a1*a1/4)));  //  {хорошая формула}
  p=a*n;  p1=a1*n;
  fout << setw(15)<<n<<setw(30)<<p<<setw(30)<<p1<<endl;
  }
return 0;
}