ring Rdp = 0, (h0(0..3),h1(0..5)), dp;       // lp; // 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, (h0(0..3),h1(0..5)), lp;       // lp; // dp;

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

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",";");

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

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