ring Rdp = 0, (h1(0..5),h0(0..3)), dp;       // lp; // dp;

string s = "ideal k = "+read("eqs");

execute(s); k;

option(redSB);

ideal j = std(k); j;
write(":w gb.dp",j);

ring Rlp = 0, (h1(0..5),h0(0..3)), lp;       // lp; // dp;

ideal i = fglm(Rdp,j); i;
write(":w gb.lp",i);
write(":a gb.lp",";");		// append ";"

list l = facstd(i); l;

write(":w gb.lp.fact",l[1]);
write(":a gb.lp.fact",";");		// append ";"