函数库 持续更新

阶乘函数库 

double fact(int n)
{
    int i = 1;
    double sum;
    sum = 1;
    for (i = 1; i <= n; i++)
    {
        sum = sum * i;
    }
    return sum;
}

判断素数函数库

int prime(int num)
{
	int i,ret=1 ;
	if (num == 1 || (num != 2 && num % 2 == 0))
	{
		ret = 0;
	}
	for (i = 2; i <= sqrt(num); i++)
	{
		if (num % i == 0)//不是素数
		{
			return 0;
			break;
		}
	}
	return ret;
  }

最大公约数和最小公倍数

int gcd(int a, int b)//最大公约数
{
	//return a % b ? gcd(b, a % b) : b;
	if (a % b)
	{
		return gcd(b, a % b);
	}
	else
	{
		return b;
	}

}
int lcm(int a, int b)//最小公倍数
{
	int t;
	t = a * b / gcd(a, b);
	return t;
}

Sin函数

//通项求法 
double Sin(double x)
{
	double sum, t;
	int n;
	t = x;
	sum = t;
	for (n = 2; n < 1000; n++)
	{
		t = -1 * t * pow(x,2) / ((2 * n - 1) * (2 * n - 2));
		sum = sum + t;
	}
	return sum;
}