ring Rdp = 0, (c(0..6),b(0..6),a(0..6)), dp;
string s = "ideal k = "+read("eqs");
execute(s); k;
option(redSB);
ideal j = std(k); j;
degree(j);
write(":w gb.dp",j);
write(":a gb.dp",";");


ring Rlp = 0, (c(0..6),b(0..5),a(0..6),b(6)), lp;
ideal i = fglm(Rdp,j); i;
degree(i);
write(":w gb.lp",i);
write(":a gb.lp",";");
list l = facstd(i); l;
write(":w gb.lp.fact.1",l[1]);
write(":a gb.lp.fact.1",";");
write(":w gb.lp.fact.2",l[2]);
write(":a gb.lp.fact.2",";");


ring Rlp = 0, (c(0..6),b(0..6),a(0..6)), lp;
string s = "ideal k = "+read("gb.lp.fact.1");
execute(s); k;
option(redSB);
ideal j = std(k); j;
degree(j);
write(":w gb.lp.1",j);
write(":a gb.lp.1",";");


ring Rlp = 0, (c(0..6),b(0..6),a(0..6)), lp;
string s = "ideal k = "+read("gb.lp.fact.2");
execute(s); k;
option(redSB);
ideal j = std(k); j;
degree(j);
write(":w gb.lp.2",j);
write(":a gb.lp.2",";");