! Juria set 
OPTION ARITHMETIC NATIVE 
LET t0=TIME 
LET a=-0.2                    
LET b=0.75    
!  INPUT a,b      
LET left = -1.4                    ! left                 
LET right = 1.4                    ! right               
LET h = (right - left)             ! height
SET WINDOW left,right,-h/2,h/2
DRAW axes    
ASK PIXEL SIZE(left,-h/2;right,h/2) px,py
LET px=px-1
LET py=py-1     
SET POINT STYLE 1
FOR u= left TO right step (right-left)/px                               
   FOR v = -h/2 to h/2 step h/py                            
      LET x1 = u                                
      LET y1 = v                                
      FOR n = 1 TO 1000                                         
         LET x = x1                           
         LET y = y1                           
         LET x1 = x*x - y*y + a       
         LET y1 = 2 * x * y + b 
         IF x1^2+y1^2>4 THEN EXIT FOR        ! determine to divergent
      NEXT n                                  
      IF n>1000 THEN PLOT POINTS: u,v
   NEXT v                                      
NEXT u     
PRINT TIME -t0                                   
END                                             

