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

ring Rlp = 0, (c(0..9),b(0..9),a(0..9)), lp;
ideal i = fglm(Rdp,j); i;
degree(i);
write(":w gb.lp",i); 
write(":a gb.lp",";");



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

list l = facstd(j); l;
write(":w gb.lp.b0.fact.1",l[1]); write(":a gb.lp.b0.fact.1",";");
write(":w gb.lp.b0.fact.2",l[2]); write(":a gb.lp.b0.fact.2",";");
write(":w gb.lp.b0.fact.3",l[3]); write(":a gb.lp.b0.fact.3",";");
write(":w gb.lp.b0.fact.4",l[4]); write(":a gb.lp.b0.fact.4",";");


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

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

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

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