Program Mandelbrot;
Global
Float cr, ci, ir, ii;
Float cr2, ci2;
Float xmax=-2.5, xmin=1, ymax=-1.2, ymin=1.2;
Float an, bn;
Float a, b;
Float a2, b2;
Float mig;
Float an2, bn2;
Float mig2;
Float Pa, Pb;
Int iter=100, tecla;
Int corx=800, cory=600;
j;
k;
j2;
k2;
n;


Begin
   Full_screen = true;
   set_mode(m800x600);
   set_fps(40,3);
   
ir = (xmax - xmin) / corx;
ii = (ymax - ymin) / cory;   
   
   For(j=0, j2=corx/2; j < corx/2; j++,j2++) 
      Frame;
        For(k=1; k<= cory ; k++)         
                     
                     cr = xmin + (j * ir);
                     ci = ymin + (k * ii); 
        
                     cr2 = xmin + (j2 * ir);
                     ci2 = ymin + (k * ii); 

                     a=Pa;
                     a2=Pa;
                     b=Pb;
                     b2=Pb;
                           
                                For(n=0; n<iter; n++) 
                                        an = (a * a) - (b * b )+ cr;
                                        bn = (2 * a * b)+ci;
                        
                                        a = an;
                                        b = bn;
                        
                                        mig=(a*a)+(b*b);

                                        If(mig>4) 
                                                put_pixel(j, cory-k, n+16);
                                                Break;
                                        End
                                        If(mig<=4&&n==iter-1) 
                                                put_pixel(j, cory-k, 1);
                                                Break;
                                        End
                                End
                
                                For(n=0; n<iter; n++) 
                                        an2 = (a2 * a2) - (b2 * b2 )+ cr2;
                                        bn2 = (2 * a2 * b2)+ci2;
                                                         
                                        a2 = an2;
                                        b2 = bn2;
                        
                                        mig2=(a2*a2)+(b2*b2);

                                        If(mig2>4) 
                                                put_pixel(j2, cory-k, n+16);
                                                Break;
                                        End
                                        If(mig2<=4&&n==iter-1) 
                                                put_pixel(j2, cory-k, 1);
                                                Break;
                                        End
                                End

                If(scan_code!=0)  
                        exit("",0);
                        End
        End

End

Loop
Frame;
                If(scan_code!=0)  
                        exit("",0);
                        End
End   
   
End