kat9:=[x1 + 2*x2 + 2*x3 + 2*x4 + 2*x5 + 2*x6 + 2*x7 + 2*x8 + 2*x9 - 1, x1^2 + 2*x2^2 + 2*x3^2 + 2*x4^2 + 2*x5^2 + 2*x6^2 + 2*x7^2 + 2*x8^2 + 2*x9^2 - x1, 2*x1*x2 + 2*x2*x3 + 2*x3*x4 + 2*x4*x5 + 2*x5*x6 + 2*x6*x7 + 2*x7*x8 + 2*x8*x9 - x2, x2^2 + 2*x1*x3 + 2*x2*x4 + 2*x3*x5 + 2*x4*x6 + 2*x5*x7 + 2*x6*x8 + 2*x7*x9 - x3, 2*x2*x3 + 2*x1*x4 + 2*x2*x5 + 2*x3*x6 + 2*x4*x7 + 2*x5*x8 + 2*x6*x9 - x4, x3^2 + 2*x2*x4 + 2*x1*x5 + 2*x2*x6 + 2*x3*x7 + 2*x4*x8 + 2*x5*x9 - x5, 2*x3*x4 + 2*x2*x5 + 2*x1*x6 + 2*x2*x7 + 2*x3*x8 + 2*x4*x9 - x6, x4^2 + 2*x3*x5 + 2*x2*x6 + 2*x1*x7 + 2*x2*x8 + 2*x3*x9 - x7, 2*x4*x5 + 2*x3*x6 + 2*x2*x7 + 2*x1*x8 + 2*x2*x9 - x8];
//time(G:=gbasis(kat9,[x1,x2,x3,x4,x5,x6,x7,x8,x9],rur));
//G;
//S:=fsolve(kat9,[x1,x2,x3,x4,x5,x6,x7,x8,x9]);
//max(seq(normal(l2norm(subst(kat9,[x1,x2,x3,x4,x5,x6,x7,x8,x9],S[j]))),j,0,size(S)-1));
S:=cfsolve(kat9,[x1,x2,x3,x4,x5,x6,x7,x8,x9]):;
size(S);
max(seq(normal(l2norm(subst(kat9,[x1,x2,x3,x4,x5,x6,x7,x8,x9],S[j]))),j,0,size(S)-1));