// ［このプログラムの目的］１個の卵型曲線を描くこと，2009年09月17日（木）

// a=1.5, b=1.2, c=8 の場合

// file name: egg_shaped_curve_SUUDUKU_2.c

#include< stdio.h>

#include< math.h>

void main(void)

{

	double x,y,a,b,c;

	int i,imax,j,m;

	double xmax,dx;

	double xx[10001],yy[10001];// 配列を格納するメモリ容量の上限に注意

	FILE *fp;

//  定数設定

        a=1.5;

	b=1.2;

        c=8;

	printf("a=%f\n\n",a);


//  他のパラメータ設定

	xmax=a;// x の最大値

	dx=xmax/2000;// x のプロット間隔

//  計算実行

	i=0;

	for(x=-xmax;x<=xmax+dx/10;x=x+dx)

	{

		i++;

		xx[i]=x;

	        if(x>xmax)

        	{

         		y=0;

        	}

        	else

        	{
         		y=(1+(x+a)/c)*sqrt(1-(x/a)*(x/a))/b;

        	}

		yy[i]=y;

		printf("i=%d,x=%f,y=%f\n",i,x,y);

	}

	imax=i;

	j=0;

	for(i=imax+1;i<=2*imax;i++)

	{
		j++;

		m=imax-j;

		xx[i]=xx[m];

		yy[i]=-yy[m];

	}

// 計算データのテキストファイルへの書き込み

	fp=fopen("egg-shaped curve.txt","w");

	if(fp==NULL)

	{

		printf("FILE OPEN ERROR\n");

	}

	else

	{

		for(i=1;i<= 2*imax;i++)

		{

			fprintf(fp,"%f,%f\n",xx[i],yy[i]);

		}

		fflush(fp);

		fclose(fp);

	}

	printf("end\n");

}// the end of the program