lnSec = SECONDS()
lnCnt=10000
lnOddArea=INT(1.2*lnCnt*log(lnCnt)/2)
DIMENSION aDigits(lnOddArea) of boolean 
k=1
FOR i=2 TO lnOddArea
    IF !aDigits[i]
        FOR j=3*i-1 TO lnOddArea STEP 2*i-1
            aDigits[j]=.t.
        ENDFOR 
        k=k+1
        IF k=lnCnt
            ? 2*i-1
            EXIT 
        ENDIF         
    ENDIF 
ENDFOR
? SECONDS() - lnSec

lnSec = SECONDS()  
  lnCnt=10000  
  lnOddArea=INT(1.2*lnCnt*log(lnCnt)/2)  
  DIMENSION aDigits(lnOddArea) of boolean   
  k=2  
  FOR i=1 TO lnOddArea while k<lnCnt
      IF !aDigits[i]  
          FOR j=(i+1)*i*2 TO lnOddArea STEP 2*i+1  
              aDigits[j]=.t.  
          ENDFOR   
          k=k+1  
      ENDIF   
  ENDFOR  
  ? 2*i+1  
  ? SECONDS() - lnSec

lnSec = SECONDS()  
  lnCnt=10000  
  lnOddArea=INT(1.2*lnCnt*log(lnCnt)/2)  
  local aDigits(lnOddArea)   
  k=1  
  FOR i=1 TO lnOddArea
      IF aDigits[i]  
          loop
      ENDIF

    FOR j=(i+1)*i*2 TO lnOddArea STEP 2*i+1  
        aDigits[j]=.t.  
    ENDFOR   
    k=k+1  

      if k=lnCnt
          exit
      endif   
  ENDFOR  
  ? 2*i+1  
  ? SECONDS() - lnSec

CLEAR
PUBLIC npop, nk
zamer('Obrazec')
npop = nk

zamer('leon1')
zamer('leon_ry')
zamer('leon_ry_leon')

CANCEL

PROCEDURE leon_ry
LPARAMETERS lnCnt
  lnOddArea=INT(1.2*lnCnt*log(lnCnt)/2)  
  DIMENSION aDigits(lnOddArea) of boolean   
  k=1  
  FOR i=2 TO lnOddArea  
      IF !aDigits  
          FOR j=3*i-1 TO lnOddArea STEP 2*i-1  
              aDigits[j]=.t.  
          ENDFOR   
          k=k+1  
          IF k=lnCnt  
              EXIT   
          ENDIF           
      ENDIF   
  ENDFOR  
RETURN 2*i-1

PROCEDURE leon_ry_leon
LPARAMETERS lnCnt
    lnOddArea=INT(1.2*lnCnt*log(lnCnt)/2)    
    local aDigits(lnOddArea)     
    k=1    
    FOR i=1 TO lnOddArea  
        IF aDigits[i]    
            loop  
        ENDIF  
      FOR j=(i+1)*i*2 TO lnOddArea STEP 2*i+1    
          aDigits[j]=.t.    
      ENDFOR     
      k=k+1    
        if k=lnCnt  
            exit  
        endif     
    ENDFOR    
RETURN 2*i+1  

PROCEDURE leon1
LPARAMETERS lnCnt
  lnMax=int(1.2*lnCnt*log(lnCnt))  
  local n(lnMax/2)  
  m.ts=1  
  for m.cnt=2 to lnCnt  
      m.tn=m.ts*2+1  
      for tt=int(m.tn*m.tn/2) to lnMax/2 step m.tn  
          n(tt)=.t.  
      next  
      m.ts=m.ts+1  
      do while n(m.ts)  
          m.ts=m.ts+1  
      enddo  
  next  
RETURN m.ts

PROCEDURE Obrazec
LPARAMETERS tnNumber  
  m.tnNumber = INT(m.tnNumber)  
  DO CASE  
  CASE m.tnNumber < 1  
      ERROR "Number must be natural figure"  
      RETURN 0  
  CASE m.tnNumber = 1  
      RETURN 2  
  ENDCASE  
  LOCAL lnCurPos, lnCurValue, llSimple, lnDivisor  
  lnCurPos = 2  
  lnCurValue = 3  
  DO WHILE .T.  
      llSimple = .T.  
      FOR m.lnDivisor = 3 TO INT(SQRT(m.lnCurValue)) STEP 2
          IF m.lnCurValue % m.lnDivisor = 0  
              llSimple = .F.  
              EXIT  
          ENDIF  
      ENDFOR  
      IF m.llSimple  
          IF m.lnCurPos = m.tnNumber  
              EXIT  
          ELSE  
              lnCurPos = m.lnCurPos + 1  
          ENDIF  
      ENDIF  
      lnCurValue = m.lnCurValue + 2  
  ENDDO  
  RETURN m.lnCurValue
  
PROCEDURE zamer
LPARAMETERS cproc
LOCAL nt
nt=SECONDS()
?cproc+' =', LTRIM(STR(&cproc(10000)))
nk=SECONDS()-nt
?IIF(cproc='Obrazec',nk,STR((nk/npop)*100,4)), IIF(cproc='Obrazec','сек.','попугаев')
?

FOR m.lnDivisor = 2 TO INT(SQRT(m.lnCurValue))