testeForte (n, b) ==
    k := 1
    q:INT:= 1 
    t:= false
    
    if b>n-1 then
        b:=powmod(b,1,n)
    if b=0 then
        return ERRO
        
    while rem(q,2)=0 repeat
        q:=quo((n-1),(2**k))
        k:=k+1
        
    i := 0
    r := powmod(b,q,n)
    if r~=1 then
        while(i<k and r ~= n-1) repeat
            r := powmod(r,2,n)
            i := i + 1
        if i = k then
            t := true

    return t

Add a code snippet to your website: www.paste.org