原理:arctan x = x - x^3/3 + x^5/5 - ... (x≤1) tan(pi/4)=1所以:arctan(1)=pi/4=1-1/3+1/5- ...;函数:function p=fun(m)%m为计算精度(有效数字)p=0;n=1;while (-1)^(n+1)/(2*n-1)>10^(-m)p=p+(-1)^(n+1)/(2*n-1);n=n+1;endp=4*p;>>fun(100)要是比较的话直接用:>>fun(100)-pi