/*

  This chunk of code is included in the "arith.c" main() routine ot test the arith-update
  routines.

  */
  /* UPDATE UPDATE    UPDATE UPDATE    UPDATE UPDATE*/

  numExponentiate(&n1, &n2, &mod, &result);
  printf("\nn1 (^) n2 = ");
  numPrintContiguous(&result);

  numExponentiate(&n2, &n1, &mod, &result);
  printf("n2 (^) n1 = ");
  numPrintContiguous(&result);

  numEuclidGCD(&n1, &mod, &result, &xn, &yn, &mod);
  printf("\ngcd(n1,mod)  =   ");
  numPrintContiguous(&result);
  printf("gcd x =?= a^-1 = ");
  numPrintContiguous(&xn);
  printf("gcd y =          ");
  numPrintContiguous(&yn);
  numMultiply(&n1, &xn, &result, &mod);
  printf("a (*) a^-1 =     ");
  numPrintContiguous(&result);

  numEuclidGCD(&n2, &mod, &result, &xn, &yn, &mod);
  printf("\ngcd(n2,mod)  =   ");
  numPrintContiguous(&result);
  printf("gcd x =?= a^-1 = ");
  numPrintContiguous(&xn);
  printf("gcd y =          ");
  numPrintContiguous(&yn);
  numMultiply(&n2, &xn, &result, &mod);
  printf("a (*) a^-1 =     ");
  numPrintContiguous(&result);

  /* END OF UPDATE    END OF UPDATE    END OF UPDATE*/