Функция egcd ($ а, $ б) {
   $ Х = 0;
   $ У = 1;
   $ И = 1;
   $ V = 0;
   в то время как ($ а! = 0) {
      $ Д = bcdiv ($ Ь, $ а);
      $ Г = bcmod ($ Ь, $ а);
      $ Т = bcsub ($ х, bcmul ($ и, $ д));
      $ П = bcsub ($ у, bcmul ($ v, $ д));
      $ Ь = $ а;
      $ А = $ г;
      $ Х = $ и;
      $ У = $ v;
      $ И = $ м;
      $ V = $ п;
   }
   возвращать массив ($ B $, х, у) $;
}
Функция modinv ($ а, $ м) {
   лист ($ г, $ х, $ у) = egcd ($ а, $ м);
   если ($ г! = 1) {
      вернуть "модульный обратный не существует";
   } Еще {
      вернуться bcmod ($ х, $ м);
   }
}