麦克劳林公式求sin

 

#include<iostream>
#include<iomanip>
#include<cmath>
using namespace std;
#define pi 3.1415926535
int main()
{
	long double x=0;
	cin>>x;
	while(x>2*pi)
	{
		x=x-2*pi;
	}
	while(x<-2*pi)
	{
		x=x+2*pi;
	}//缩小x范围让求值更精确；
	int sign=-1,j=0;//sign-符号，j-两项间阶乘变化量
	long double sinx=x;//sinx为输出值
	long double v=x;//v保存最后一项
	while(fabs(v)>=0.5e-13)//最后一项不满足时持续循环
	{
		j=j+2;//j每次循环都改变
		v=sign*v*(x*x/(j*(j+1)));//v进行循环改变
		sinx=sinx+v;//每次sinx加上新增的最后一项
	}
	cout<<fixed<<setprecision(13)<<sinx;
}

注意：1.define Pi值

2.利用三角变换缩小x值减小误差

3.fabus()取绝对值来进行循环判断

4.利用最后一项变换，来每次都加最后一项v

5.给出最后一项的变化公式