function pgcd(nb1,nb2)

{

    var nb3;

    while(nb2>0)

    {
        nb3=nb1;
        nb1=nb2;

        nb2=(nb3-nb1*Math.floor(nb3/nb1));

    }
    return nb1;
}

/* Récursivité utilisée   */
int pgcd(int m, int n)
{
  if (n==0) return m;
  else return pgcd(n, m%n);
}
Sans récursivité car dangereux pour la stack:
DWORD pgcd(DWORD a, DWORD b)
{
  DWORD t;
  if(!a || !b) return 0;
  if(a < b) {t = a; a = b; b = t;}
  for(;;) {
    if((t = a % b) == 0) return b;
    a = b; b = t;
  }
}

(* Calcule le pgcd des entiers a et b. Exemple d'utilisation : pgcd 4 8 donne 4*)
let rec pgcd a = function
    | b when a mod b = 0 -> b
    | b -> pgcd b (a mod b)

(* pgcd de deux nombres. Exemple d'utilisation : pgcd 2 4 renvoie 2 *)
let rec pgcd m = function
      | 0 -> m
      | n -> pgcd n (m mod n)

    Public Function PGCD(ByVal a As Integer, ByVal b As Integer) As Integer
        If a = b Then
            Return a
        End If
        Dim c As Integer
        If b > a Then
            c = a
            a = b
            b = c
        End If
        Dim rest As Integer
        Math.DivRem(a, b, rest)
        If rest = 0 Then a = b
        While rest <> 0
            Math.DivRem(a, b, rest)
            a = b
            b = rest
        End While
        Return a
    End Function

Function Pgcd2(ByVal A As Long, ByVal B As Long) As Long
'Calcul le PGCD de deux nombres
'Variante optimisant l'Algo d'Euclide
'Déclarations
Dim R As Long
'validité paramètre
If A = 0 And B = 0 Then Exit Function 'Pas de solution
If A < B Then R = B: B = A: A = R 'inverse les valeurs
'Algo
Do While Abs(B) >= 1
    R = A - Int(A / B + 0.5) * B 'Ici le +0.5 permet de réduire mathématiquement le nb de boucle nécessaire
    A = B
    B = R
Loop
Pgcd2 = Abs(A)
End Function

function PGCD(a, b: Integer): Integer;
begin
  if b = 0 then
    Result := a
  else
    Result := PGCD(b, a mod b);
end;