#include<stdio.h>

void main(void)
{
	FILE *fp;
	float x,y,Newx,Newy;
	int X[640][480];
	int i,j,k;

	for(i=0;i<640;i++){
printf("On iterate %d\n",i);
	for(j=0;j<480;j++){
	X[i][j]=0;
	x=((double) 4*i)/640-2;
	y=((double) -4*j)/480+2;
		for(k=0;k<100;k++){
			if(x<0 && y<=0){
				Newx=2*x+1;
				Newy=2*y+1;
			}
			if(y>0){
				Newx=2*y-1;
				Newy=-2*x+1;
			}
			if(x>=0 && y<=0){
				Newx=2*x-1;
				Newy=2*y+1;
			}
			x=Newx;y=Newy;
			if(x*x+y*y>100000){
				X[i][j]=k;
				break;
			}
		}
	}
	}
	fp=fopen("rug3.dat","w");
	for(i=0;i<640;i++){
		for(j=0;j<480;j++){
		fprintf(fp,"%d ",X[i][j]);
		if(j==479) fprintf(fp,"\n");
		}
	}
	fclose(fp);
}