3_3.cpp

清华大学郑莉《C++语言程序设计》课件、C++教材源代码、实验参考程序全集

//3_3.cpp
#include<iostream>
using namespace std;
int main()
{
  double a,b;
  double arctan(double x)  ;
  a=16.0*arctan(1/5.0)  ; 
  b=4.0*arctan(1/239.0)  ; 
    //注意：因为整数相除结果取整，如果参数写1/5，1/239，结果就都是0
  cout<<"PI="<<a-b<<endl;
}
double arctan(double x)
{
  int i;
  double r,e,f,sqr;
  sqr=x*x;
  r=0;
  e=x;
  i=1;
  while(e/i>1e-15)
  {
    f=e/i;
    r=(i%4==1)? r+f : r-f   ;
    e=e*sqr;
    i+=2;
  }
  return  r ;
}