infolevel[solve] := 3;
read "CycloSolve.mpl";

read "sys225.txt";
X := CycloSolve(A,b,M,e,5);
# check solutions
convert( map(rem,map(numer,A.X-b),M,e), set );

read "sys196.txt";
A := 182*A: b := 182*b: # clear fractions
X := CycloSolve(A,b,M,e,3);
# check solutions
convert( map(rem,map(numer,A.X-b),M,e), set );