基于过程的sin函数的计算

基于泰勒展开式，我们可以方便的计算正弦函数的值，那么如果要求不能使用函数呢，可能会带来点小麻烦，下面是纯过程的正弦函数的计算，精度可控

<span style="font-size:24px;">#include <iostream>
#include <stdio.h>
using namespace std;

int main()
{
double x;
double sin_res=0;//正弦结果
double sin_res_before=0;
int sin_presicion=1e-4;//sin循环次数
int PI_precision=1000000;//PI循环次数
int flag;//正负号
int deep;//当前阶数
double mol;//分子
double den;//分母
double PI=0;
//计算PI
{
flag=1;
mol=4.0;
den=1.0;
for(int i=0; i<PI_precision; i++)
{
//printf("num=%lg\n",(flag)*(mol/den));
PI+=(flag)*(mol/den);
flag=-flag;
den+=2.0;
}
}
x=PI*0.25;
//第一问
{
printf("k ");
for(double i=0; i<=1.0; i+=0.1)
{
printf("%-7lg",i);

}
printf("\n");
printf("sin ");
for(double i=0; i<=1.0; i+=0.1)
{
x=i*PI/2.0;
sin_res=0;
sin_res_before=-1;
flag=1;//初始化
deep=1;
mol=x;
den=1.0;
while((sin_res-sin_res_before)>sin_presicion)
{
//printf("mol=%lg,den=%lg\n",mol,den);
sin_res_before=sin_res;
sin_res+=flag*(mol/den);
deep+=2;
mol*=x*x;
den*=deep*(deep-1);
flag=-flag;
}
printf("%.4lf ",sin_res);
}
printf("\n");
}

//第二问
{

}
/*
//计算级数展开式
{
flag=1;//初始化
deep=1;
mol=x;
den=1.0;
for(int p=0; p<sin_presicion; p++)
{
//printf("mol=%lg,den=%lg\n",mol,den);
sin_res+=flag*(mol/den);
deep+=2;
mol*=x*x;
den*=deep*(deep-1);
flag=-flag;
}
}
printf("sin(%lg)=%lg\n",x,sin_res);
*/
return 0;
}
</span>