#include <math.h>

/* square region */
#define XSIZE 1000
#define YSIZE XSIZE

#define PI 3.1415927
#define NITER 100