// ［このプログラムの目的］伊藤氏の方法による１個の卵形原型曲線を描くこと，
// 2010年5月6日（木）
// a=0.5, b=0.37 の場合
// file name: C_program_3.c

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

void main(void)
{
	double x,y,a,b,pi;
	double t,dt;
	int i,imax;
	double xx[10001],yy[10001];//  メモリ容量の上限に注意
	FILE *fp;
// 定数設定
        pi=3.1415927;
        a=.5;
	b=.37;
	printf("a=%f\n\n",a);
// 他のパラメータ設定
	dt=pi/2000;// t のプロット間隔
	i=0;
// 計算実行
	for(t=-2*pi;t<=2*pi;t=t+dt)
	{
		i++;
		x=a*cos(t);
		y=b*cos(t/4)*sin(t);
		xx[i]=x;
		yy[i]=y;
		printf("i=%d,x=%f,y=%f\n",i,x,y);
	}
	imax=i;
// 計算データのテキストファイルへの書き込み
	fp=fopen("egg_shaped_curve.txt","w");
	if(fp==NULL)
	{
		printf("FILE OPEN ERROR\n");
	}
	else
	{
		for(i=1;i<=imax;i++)
		{
			fprintf(fp,"%f,%f\n",xx[i],yy[i]);
		}
	fflush(fp);
	fclose(fp);
	}
	printf("end\n");
}// the end of the program