// ［このプログラムの目的］
// １個のリンゴ形曲線を、
// 媒介変数を用いて計算して描くこと
// 2010年3月29日（月）

// file name: apple.c

#include< stdio.h>

#include< math.h>

void main(void)

{
	double a,pi;
// a は元のカージオイドの定数、pi は円周率

	double b;
// 縦方向に圧縮する圧縮係数

	double p, q;
// 指数関数による横軸と縦軸の伸縮変換率

	double t,dt;
// 変換前のカージオイドの位相角[ラジアン]
// と位相角の変分[ラジアン]

	double R;
// 変換前のカージオイドの曲座標の動径

	double tmin,tmax;
// 変換前のカージオイドの位相角の最小値
// と最大値[ラジアン]

	int i,imax;

	double xx[10001],yy[10001];
// メモリ容量の上限に注意

	FILE *fp;

//  定数設定

	pi=3.14159265;
	a=1.0;
	b=1.2;
	p=0.15;
	q=0.08;

//  他のパラメータ設定

	tmin=-pi/2;
	tmax=3*pi/2;
	dt=(tmax-tmin)/4000;
// 変換前の位相角 t のプロット間隔

//  計算実行

	i=0;
	for(t=tmin;t<=tmax+dt;t=t+dt)
		{
		i++;
		R=a*(1-sin(t));
xx[i]=R*cos(t)*exp(-p*(t-pi/2)*(t-pi/2));
yy[i]=b*R*sin(t)*exp(-q*(t-pi/2)*(t-pi/2));

printf("i=%d,x=%f,y=%f\n",i,xx[i],yy[i]);
		}
	imax=i;

// 計算データ[ｎ角の星形曲線の(x,y)座標]を
// テキストファイル（apple.txt）に書き込む

	fp=fopen("apple.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