// ［このプログラムの目的］伊藤氏の方法による4個の卵型曲線を描くこと，2008年3月20日（木）
// a=0.5, b=0.11～0.5 の場合
// file name: C_program_2.c

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

void main(void)
{
	double x,y,a,b,db,bmax,pi;
	double t,dt;
	int i,imax,n,nmax;
	double xx[10001],yy[10][10001];// 配列を格納するメモリ容量の上限に注意

	FILE *fp;

// 定数の設定
        pi=3.1415927;
        a=.5;
	bmax=a;
	db=0.13;// パラメータb の選択間隔

// 他のパラメータ設定
	dt=pi/1000;// t のプロット間隔

	i=0;

// 計算実行
	for(t=-pi;t<=pi;t=t+dt)
	{
		i++;
		n=0;

		for(b=bmax-3*db;b<=bmax+0.01;b=b+db)
		{
			n++;

			x=a*cos(t);
			y=b*cos(t/4)*sin(t);
			xx[i]=x;
			yy[n][i]=y;

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

	imax=i;
	nmax=n;

// 計算データのテキストファイルへの書き込み
	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,%f,%f,%f\n",xx[i],yy[1][i],yy[2][i],yy[3][i],yy[4][i]);
		}

		fflush(fp);
		fclose(fp);
	}

	printf("end\n");
}// the end of the program

～ブラウザの「戻る」で 前に開いたページに戻ります～

～または、下のボタンを押して 本サイトの前ページに戻ります～

戻る