read(`/amaterasu/sd2e/labm/coffey/escher/test/specific.txt`);
g := 6;
read(`/amaterasu/sd2e/labm/coffey/escher/test/specific_6_L2.m`);
read(`/amaterasu/sd2e/labm/coffey/escher/test/initial_data_6.m`);
tosolve := TDerivativeAtZero(L,g);
save(tosolve,`tosolve6.m`);
index := indets(tosolve,name):
for i in index do
  assume(i,real);
od;
index;
nops(");
ans := solve(tosolve);
save(ans,`ans6.m`);

restart;
read(`/amaterasu/sd2e/labm/coffey/escher/test/specific.txt`);
read(`/amaterasu/sd2e/labm/coffey/escher/test/ans6.m`);
test_set := select(x->rhs(x)<>lhs(x) and rhs(x)<>0,ans);
ZeroTerms := select(x->rhs(x) = 0,ans);
bad := select(y->member(true,map(x->has(x,I),[coeffs(expand(rhs(y)))])),test_set);
t := proc(i) if binary[i] then i fi end;
L := NULL:
for j from 1 to nops(bad) do
if  type(expand(rhs(bad[j])),`+`) then 
  d := [op(expand(rhs(bad[j])))];
else 
  d := [rhs(bad[j])];
fi;
binary := map(x->has(x,I),map(x->coeffs(x),d));
Binary := [seq(t(i),i=1..nops(binary))]:
L := L,convert({seq(op(Binary[i],d),i=1..nops(Binary))},`+`);
od:
{L};
ToZero := map(x->x=0,");
ans2 := solve(ToZero union ans);
bad2 := select(y->member(true,map(x->has(x,I),[coeffs(expand(rhs(y)))])),select(x->rhs(x)<>lhs(x) and rhs(x)<>0,ans2));
dim := nops(select(x->rhs(x)=lhs(x) and rhs(x)<>0,ans2));
save(dim,`dim6.txt`);

done;