53.c

来自「计算概论阅读材料」· C语言 代码 · 共 23 行

/*用级数展开式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.))) ;
}