// compare to localgb.m2 (or localgb.m2-hold)
//R = ZZ/101[a..d,MonomialOrder=>Weights=>{-1,-1,-1,-1}]
//I = ideal(a-b*a^2-a^3, b-a*b-c^3, a*b-b*d^4-c^5)

ring R = 101, (a,b,c,d),ds;
ideal I = a-b*a^2-a^3, b-a*b-c^3, a*b-b*d^4-c^5;
matrix m = matrix(I);
m;
matrix m2 = syz(m);
m2;
matrix m2gb = minbase(m2);
m2gb;

std(I);
module M=syz(I);
module M2=syz(M);
M[1];
matrix m=M;
m[1];
resolution C=mres(I,4);
C;
print(C);
minbase(I);
std(I);
size(m2gb[1,3]); // not what I thought it would be...


