NumPoints = 1500
MaxIters = 2000
for (double i=0; i < NumPoints; i++) {
    for (double j=0; j < NumPoints; j++) {
      c.real = -2.0 + 3.0*i/NumPoints + 1.0e-7 /* Minimal displacement when i=0 */;
      c.imag = -1 + 2*j/NumPoints) + 1.0e-7 /* Minimal displacement when j=0 */;
      z=c;
      for (int iter=0; iter<NumIters; iter++){
        ztemp=(z.real*z.real)-(z.imag*z.imag)+c.real;
        z.imag=z.real*z.imag*2+c.imag;
        z.real=ztemp;
        if ((z.real*z.real+z.imag*z.imag)>4.0) {
          numoutside++;
          break;
        }
      }
	}
}

double area = 3.0*2*(NumPoints*NumPoints-numoutside)/(NumPoints*NumPoints);
double error=area/NumPoints;