#include <stdio.h>
 #include <math.h>
 int main()
 {
  /* screen ( integer) coordinate */
  int iX,iY;
  const int iXmax = 2000, iXmin=0; 
  const int iYmax = 2000, iYmin=0;
  int iWidth=iXmax-iXmin+1,
  iHeight=iYmax-iYmin+1,
  /* number of bytes = number of pixels of image * number of bytes of color */
  iLength=iWidth*iHeight*3,/* 3 bytes of color  */
  index; /* of array */
  /* world ( double) coordinate = parameter plane*/
  double Cx,Cy;
  const double CxMin=-2.5;
  const double CxMax=1.5;
  const double CyMin=-2.0;
  const double CyMax=2.0;
  /* */
  double PixelWidth=(CxMax-CxMin)/iXmax;
  double PixelHeight=(CyMax-CyMin)/iYmax;
  /* color component ( R or G or B) is coded from 0 to 255 */
  /* it is 24 bit color RGB file */
  const int MaxColorComponentValue=255; 
  FILE * fp;
  char *filename="LCMM.ppm";
  char *comment="# ";/* comment should start with # */
  /* Z=Zx+Zy*i  ;   Z0 = 0 */
  double Zx, Zy;
  double Zx2, Zy2; /* Zx2=Zx*Zx;  Zy2=Zy*Zy  */
  /*  */
  int Iteration, PreviousIter;
  const int IterationMax=2000;
  /* bail-out value , radius of circle ;  */
  const double EscapeRadius=1000;
  double ER2=EscapeRadius*EscapeRadius;
  /* dynamic 1D array for 24-bit color values */    
  unsigned char *array;
  /*-------------------------------------------------------------------*/
  array = malloc( iLength * sizeof(unsigned char) );
  if (array == NULL)
    {
    fprintf(stderr,"Could not allocate memory");
    getchar();
    return 1;
    }
    else 
     {   fprintf(stderr,"I'm working. Please wait (:-))\n ");
         /* fill the data array with white points */       
         for(index=0;index<iLength-1;++index) array[index]=255;
     }
  /* ---------------------------------------------------------------*/
  /* first coat of paint */
  for(iY=0;iY<iYmax;iY++)
    {
    Cy=CyMin + iY*PixelHeight;
    if (fabs(Cy)< PixelHeight/2) Cy=0.0; /* Main antenna */
    for(iX=0;iX<iXmax;iX++)
      {         
       Cx=CxMin + iX*PixelWidth;
       /* initial value of orbit = critical point Z= 0 */
       Zx=0.0;
       Zy=0.0;
       Zx2=Zx*Zx;
        Zy2=Zy*Zy;
       /* */
       for (Iteration=0;Iteration<IterationMax && ((Zx2+Zy2)<ER2);Iteration++)
        {
         Zy=2*Zx*Zy + Cy;
         Zx=Zx2-Zy2 +Cx;
         Zx2=Zx*Zx;
         Zy2=Zy*Zy;
        };
       /* plot point of Level Curve */
       if (iX!=0 && Iteration!=PreviousIter)
       { 
         array[((iYmax-iY-1)*iXmax+iX)*3]=0;
         array[((iYmax-iY-1)*iXmax+iX)*3+1]=0;
         array[((iYmax-iY-1)*iXmax+iX)*3+2]=0;
         PreviousIter=Iteration;                       
       }
     }
    }
   /* second coat of paint */
   for(iX=0;iX<iXmax;iX++)
   {         
     Cx=CxMin + iX*PixelWidth;
     for(iY=0;iY<iYmax;iY++)
     {
      Cy=CyMin + iY*PixelHeight;
      if (fabs(Cy)< PixelHeight/2) Cy=0.0; /* Main antenna */
      /* initial value of orbit = critical point Z= 0 */
      Zx=0.0;
      Zy=0.0;
      Zx2=Zx*Zx;
      Zy2=Zy*Zy;
      /* */
      for (Iteration=0;Iteration<IterationMax && ((Zx2+Zy2)<ER2);Iteration++)
       {
        Zy=2*Zx*Zy + Cy;
        Zx=Zx2-Zy2 +Cx;
        Zx2=Zx*Zx;
        Zy2=Zy*Zy;
       };
       /* plot point of Level Curve */
       if (iX!=0 && Iteration!=PreviousIter)
        { 
         array[((iYmax-iY-1)*iXmax+iX)*3]=0;
         array[((iYmax-iY-1)*iXmax+iX)*3+1]=0;
         array[((iYmax-iY-1)*iXmax+iX)*3+2]=0;
         PreviousIter=Iteration;                       
        }
       }
      }
   /* write the whole data array to ppm file in one step */      
   /*create new file,give it a name and open it in binary mode  */
   fp= fopen(filename,"wb"); /* b -  binary mode */
   if (fp == NULL){ fprintf(stderr,"file error"); }
     else
      {
       /*write ASCII header to the file*/
      fprintf(fp,"P6\n %s\n %d\n %d\n %d\n",comment,iXmax,iYmax,MaxColorComponentValue);
      /*write image data bytes to the file*/
      fwrite(array,iLength ,1,fp);
      fclose(fp);
      fprintf(stderr,"file saved\n");
      getchar();
     }
    free(array);
    getchar();
    return 0;
 }