Alpha 1/シェルピンスキー・ガスケット

Alpha 1

https://scratchpad.fandom.com/ja/wiki/Alpha_1

シェルピンスキー・ガスケット

注:rand.datは[0,1]で乱数が収められているファイルです。

#include <stdio.h>

#define CALC_REPEAT 1000000
#define SIDE_LENGTH 1.0
#define X0 0.5
#define Y0 0.1

int main(void)
	{
	FILE *rd,*fs;
	double x=X0,y=Y0,r;
	int count=0;
	
	if(((rd=fopen("rand.dat","r")) != NULL) && ((fs=fopen("plot_gasket.dat","w+"))!=NULL))
		{
		fprintf(fs,"%lf %lf\n",x,y);
		
		while(count<CALC_REPEAT)
			{
			fscanf(rd,"%lf",&r);
			
			switch ((int)(r*3.0))
				{
				case 0:
					x=x/2;
					y=y/2;
					
					break;
				case 1:
					x=(x+SIDE_LENGTH)/2;
					y=y/2;
					
					break;
				case 2:
					x=(x+SIDE_LENGTH/2)/2;
					y=(y+SIDE_LENGTH*1.73205/2)/2;
					
					break;
				}
			
			fprintf(fs,"%lf %lf\n",x,y);
			
			count++;
			}
		
		fclose(rd);
		fclose(fs);
		}
	
	return 0;
	}

コメント 投稿

2008/07/24 13:29

プロフィール

名前:Tommy6
ニックネーム:(　・ω・)
性別:男
生年月日:
出身地:北
血液型:B

2008年07月