| LIB "homolog.lib";
int p = printlevel;
printlevel = 1;
ring r = 0,(x,y,z),dp;
ideal i = x2y,y2z,z3x;
module E = Ext_R(1,i); //computes Ext^1(r/i,r)
→ // Computing Ext^1:
→ // Let 0<--coker(M)<--F0<--F1<--F2<--... be a resolution of M,
→ // then F1*-->F2* is given by:
→ x2, -yz,0,
→ 0, z3, -xy,
→ xz2,0, -y2
→ // and F0*-->F1* is given by:
→ y2z,
→ x2y,
→ xz3
→
→ // dimension of Ext^1: -1
→
is_zero(E);
→ 1
qring R = std(x2+yz);
intvec v = 0,2;
printlevel = 2; //shows what is going on
ideal i = x,y,z; //computes Ext^i(r/(x,y,z),r/(x2+yz)), i=0,2
list L = Ext_R(v,i,1); //over the qring R=r/(x2+yz), std and kbase
→ // Computing Ext^0:
→ // Let 0<--coker(M)<--F0<--F1<--F2<--... be a resolution of M,
→ // then F0*-->F1* is given by:
→ z,
→ y,
→ x
→ // and F-1*-->F0* is given by:
→ 0
→
→ // dimension of Ext^0: -1
→
→ // columns of matrix are kbase of Ext^0 in F0*:
→ 0
→
→ // Computing Ext^2:
→ // Let 0<--coker(M)<--F0<--F1<--F2<--... be a resolution of M,
→ // then F2*-->F3* is given by:
→ x,-y,z, 0,
→ z,x, 0, z,
→ 0,0, x, y,
→ 0,0, -z,x
→ // and F1*-->F2* is given by:
→ y,-z,0,
→ x,0, -z,
→ 0,x, -y,
→ 0,z, x
→
→ // dimension of Ext^2: 0
→ // vdim of Ext^2: 1
→
→ // columns of matrix are kbase of Ext^2 in F2*:
→ x,
→ -z,
→ 0,
→ 0
→
printlevel = p;
|