ring Rdp = 0, (g(0..7),h(0..7),A), 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, (g(0..7),h(0..7),A), 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.1",l[1]);
write(":a gb.lp.1",";");
write(":w gb.lp.2",l[2]);
write(":a gb.lp.2",";");
write(":w gb.lp.3",l[3]);
write(":a gb.lp.3",";");
write(":w gb.lp.4",l[4]);
write(":a gb.lp.4",";");