/* 
c console program:
          
----------------------------------------------------
 1. draws parameter plane for Fc(z)=z*z +c
  color shows in which quadrant zn lands. 4 colors are used because there are 4 quadrants : 
-------------------------------         
2. technic of creating ppm file is  based on the code of Claudio Rocchini
http://en.wikipedia.org/wiki/Image:Color_complex_plot.jpg
create 8 bit color ( shades of gray)  graphic file ,  portable pixmap file = PPM 
see http://en.wikipedia.org/wiki/Portable_pixmap
to see the file use external application ( graphic viewer)
---------------------------------

Adam Majewski
fraktal.republika.pl
{{Date|2010|04|02}}
*/
#include <stdio.h>
int main()
{
         /* screen ( integer) coordinate */
       int iX,iY;
       const int iXmax = 4000; 
       const int iYmax = 4000;
       /* 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; 
       static unsigned char color[1]; /* 8 bit color =  1 byte */
       FILE * fp;
       char *filename="d_9.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;
       const int IterationMax=9;
       /* bail-out value , radius of circle ;  */
       const double EscapeRadius=1000;
       double ER2=EscapeRadius*EscapeRadius;
       /*create new file,give it a name and open it in binary mode  */
       fp= fopen(filename,"wb"); /* b -  binary mode */
       /*write ASCII header to the file - 8 bit color*/
       fprintf(fp,"P5\n %s\n %d\n %d\n %d\n",comment,iXmax,iYmax,MaxColorComponentValue);
       /* compute and write image data bytes to the file*/
       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; Iteration++)
                       {
                           Zy=2*Zx*Zy + Cy;
                           Zx=Zx2-Zy2 +Cx;
                           Zx2=Zx*Zx;
                           Zy2=Zy*Zy;
                       };
                       /* compute  pixel color (8 bit = 1 byte) */
                      // if (Iteration==IterationMax) color[0]=0; /*  interior of Mandelbrot set = black */
                       /* exterior of Mandelbrot set = binary decomposition */             
                         // else 
                       if ((Zx>0) && (Zy>0)) color[0]=50;   /* 1-st quadrant */
                       if ((Zx<0) && (Zy>0)) color[0]=100; /* 2-nd quadrant */
                       if ((Zx<0) && (Zy<0)) color[0]=150; /* 3-rd quadrant */
                       if ((Zx>0) && (Zy<0)) color[0]=200; /* 4-th quadrant */
                       /*write color to the file*/
                       fwrite(color,1,1,fp);
               }
       }
       fclose(fp);
       printf("file saved\n");
       getchar();
      return 0;
}