#include #include #define NMAX 1000 typedef struct { double re; double im; } CPLX; int main(){ CPLX z,r,c; double x,y,dx,dy; dx = 0.01; dy = dx; int i,j,iter; for(i=0;i<=4/dx;i++) { for(j=0;j<=4/dy;j++) { x=-2+i*dx; y=2-j*dy; //r.re=0; r.im=0; //c.re=x; c.im=y; r.re=x; r.im=y; c.re=0.5; c.im=0.33; iter=0; do { z.re=r.re; z.im=r.im; r.re=z.re*z.re-z.im*z.im+c.re; r.im=2*(z.re*z.im)+c.im; iter++; } while(sqrt(z.re*z.re+z.im*z.im)<4 && iter