ring Rdp = 0, (h3(0..5),h2(0..3),h0(0..3),h1(0..5)), dp;

string str = "ideal k1 = "+read("eqs.B");
execute(str);
k1;


string s2 = "ideal k2 = "+read("gb.lp.fact.2");
execute(s2);
k2;

ideal k = k1,k2;

option(redSB);

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



ring Rlp = 0, (h3(0..5),h2(0..3),h1(0..5),h0(0..3)), lp;

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