/* M A N D E L B R O T */

#include "graphics.h"
#include "math.h"
float xmin,xmax,ymin,ymax;
void MSB(cxmin,cxmax,cymin,cymax,nmax)
float cxmin,cxmax,cymin,cymax;
int nmax;
{
float sx=500.0,sy=400.0;
float cx,cy,x,y,xx,yy,dx,dy,z,L=4.0;
int i,j,n,color;
dx=(cxmax-cxmin)/sx;
dy=(cymax-cymin)/sy;
for(i=0;i<sx;i++)
  {
   cx=cxmin+i*dx;
   for(j=0;j<sy;j++)
    {
     cy=cymin+j*dy;
     x=0;
     y=0;
     for(n=0;n<nmax;n++)
      {
       xx=x*x-y*y+cx;
       yy=2*x*y+cy;
       z=xx+yy;
       if(sqrt(xx*xx+yy*yy)>L)break;
       x=xx;
       y=yy;
       }
       color=n%256;
       putpixel(i+50,j+40,color);
       }
      }
    }
    main()
    {
     int gdriver=DETECT,gmode;
     printf("IMPUT GRAPH SMALLIST X SMALLIST Y BIGGEST X BIGGEST Y!!");
     scanf("%f%f%f%f",&xmin,&xmax,&ymin,&ymax);     
     initgraph(&gdriver,&gmode,"");     
     MSB(xmin,xmax,ymin,ymax,100);
     getch();
    }
/*
何谓Mandelbrot??
Mandelbrot是由几何学家mandelbrot发现的
在复平面上的迭带z<=z*z+z会把复平面上的点分成3种类型：
  1：|z|-〉∞
  2：|z|-> 0
  3:不同于1与2，即运动无规律
“3”的类型就是Mandelbrot集合

-------------------------------------------------------------------------
[url=http://www.fractal.com.cn]分形频道 *** www.fractal.com.cn ***[/url]
-------------------------------------------------------------------------
*/