プログラマ専用SNS ミクプラ

コード:

#include<stdio.h>
#include<math.h>

double mysin(double,double);

int main(void) {
	double x,y,PI=atan(1.0)*4.0,k=1.0;

	printf("角度を度数法で入力してください\n");
	scanf("%lf",&y);
	x=(2.0*PI/360.0)*y;

	printf("sin(%lf)=%lf",x,1+mysin(x,k));

	return 0;
}

double mysin(double x,double k) {
	double my_power(double);
	double kaijou(double);

	return (my_power(x)/kaijou(k));
}

double my_power(double x) {
	double a=1.0,b=a,c=0.0;
	int i,j;

	for(j=1;j<=8;++j) {
		for(i=20-j;1<=i;--i)
			a*=x;
		if(j%2==0)
			b=-b;
		c+=b;
	}

	return c;
}

double kaijou(double k) {
	int i,j;
	double a=1.0,b=0.0;

	for(j=1;j<=8;++j) {
		for(i=20-j;1<=i;--i)
			a*=i;
		b+=a;
	}

	return b;
}

コード:

int main()
{
	弧度を取得、

以下を複数回行う(nは試行毎に2ずつ増えていく)
	バッファに、(x^n/n!)を代入する(nは一以上の奇数)

	もし、試行回数が2で割り切れるときは
	答えに、バッファを加算、

	割り切れないときは
	答えに、バッファを減算。
}

コード:

#include <stdio.h>
#include <math.h>

double mysin(double);

int main(void) {
	double result,resultreal;
	double angle;
	/*角度の入力*/
	printf("角度を度数法で入力してください。\n>");
	scanf("%lf",&angle);
	/*弧度法に変換*/
	angle=angle*3.14159265358979323846/180.0;
	/*検証用に本物のsinで計算*/
	resultreal=sin(angle);
	/*自作関数で計算*/
	result=mysin(angle);
	/*結果の表示*/
	printf("本物:%1.15f\n",resultreal);
	printf("自作:%1.15f\n",result);
	return 0;
}

double mysin(double angle) {		/*angleは弧度法*/
	double result;					/*結果を代入する変数*/
	double temp;					/*現在の項の値の変数*/
	double i;						/*ループ変数*/
	double maeresult;				/*前回の結果（終了判定用）*/
	result=temp=angle;				/*一番目の項*/
	for(i=3;;i+=2) {				/*ループ*/
		/*この時点でtempはangle^(i-2)/(i-2)!のはず（符号は除く）*/
		temp=-temp;					/*tempの符号を反転する*/
		temp=temp*angle*angle;		/*tempにangle^2を掛ける*/
		temp=temp/(i-1)/i;			/*tempの分母をi!にする*/
		maeresult=result;			/*前回の結果を保存する*/
		result+=temp;				/*現在の項を加算する*/
		if(result==maeresult)break;	/*結果が変わらなくなったら抜ける*/
	}
	return result;					/*結果を返す*/
}