pgcd2(a,b):= {
  local r,u,v,r1,u1,v1,q;
    (r,u,v,r1,u1,v1):=(a,1,0,b,0,1);

  while (r1!=0) {
    q:=iquo(r,r1);
    (r,u,v,r1,u1,v1):=(r1,u1,v1,r-q*r1,u-q*u1,v-q*v1);
    }
  return(r,u,v);
  }