基于过程的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>