53.c

计算概论阅读材料

/*用级数展开式arctan(x)(要求求至某项绝对值<10--6为止)以及*/
/*求p，并将结果与直接计算结果printf("%lf",6*arctan(1.0/sqrt(3)));*/
/* 进行比较*/
#include <math.h>
#include <stdio.h>
void main(void)
{
	int n=1 ;
	float atx=0.0,t ,s ,x ,x2 ;
	x=1.0/sqrt(3.) ;
	x2=x*x ;
	t=x ;
	s=t ;
	while( fabsl(t)>1e-6 ){ //1e-7
		atx+=t ;
		s*=-x2 ;
		n+=2 ;
		t=s/n ;
		}
	printf("pi=%f\n",6*atx ) ;
	printf("pi=%f\n",6*atan(1./sqrt(3.))) ;
}