先 #include<math.h>
反3角函数有 acos(double),asin(double),atan(double),atan(double,double),返回值 double
cos(double x)
在C/C++中使用pi (π) 值
马青公式:
π=16arctan1/5-4arctan1/239
走进无限美妙的数学世界 http://amuseum.cdstm.cn/AMuseum/math/4/46/4_46_1004.htm
马青公式计算圆周率程序 - CSDN博客 https://blog.csdn.net/MHL_1208980380/article/details/54284631
用马青公式计算圆周率最精简的C++代码(万进制):
#include<stdio.h> #include<iostream> using namespace std; int main(void) { //本程序每四位数输出结果,如果请求计算的位数不是4的整数倍,最后输出可能会少1~3位 long a[2]={956,80},b[2]={57121,25},i=0,j,k,p,q,r,s=2,t,u,v,N,M=10000; printf("%9cMachin%6cpi=16arctan(1/5)-4arctan(1/239) Please input a number. ",32,32); cin>>N,N=N/4+3; long *pi=new long[N],*e=new long[N]; while(i<N)pi[i++]=0; while(--s+1) { for(*e=a[k=s],i=N;--i;)e[i]=0; for(q=1;j=i-1,i<N;e[i]?0:++i,q+=2,k=!k) for(r=v=0;++j<N;pi[j]+=k?u:-u)u=(t=v*M+(e[j]=(p=r*M+e[j])/b[s]))/q,r=p%b[s],v=t%q; } while(--i)(pi[i]=(t=pi[i]+s)%M)<0?pi[i]+=M,s=t/M-1:s=t/M; for(cout<<"3.";++i<N-2;)printf("%04ld",pi[i]); delete []pi,delete []e,cin.ignore(),cin.ignore(); return 0; }