H0=-y*D[f,y]+x0*D[f,x0];
H1=-x2*D[f,x2]+x4*D[f,x4]-x5*D[f,x5]+x7*D[f,x7]+x8*D[f,x8]-x9*D[f,x9];
H2=-x1*D[f,x1]+x2*D[f,x2]-x3*D[f,x3]+x5*D[f,x5]-x6*D[f,x6]+x9*D[f,x9];
H3=-2*f-x0*D[f,x0]+x1*D[f,x1]-x5*D[f,x5]-x7*D[f,x7]-x8*D[f,x8]-x9*D[f,x9];
H4=-x1*D[f,x1]-x2*D[f,x2]+x3*D[f,x3]-x4*D[f,x4]+x5*D[f,x5]+x7*D[f,x7];
H5=-x3*D[f,x3]-x5*D[f,x5]+x6*D[f,x6]-x7*D[f,x7]+x8*D[f,x8]+x9*D[f,x9];
Hom=-6*f-2*y*D[f,y]-x0*D[f,x0]-x1*D[f,x1]-x2*D[f,x2]-x3*D[f,x3]-x4*D[f,x4]-x5*D[f,x5]-x6*D[f,x6]-x7*D[f,x7]-x8*D[f,x8]-x9*D[f,x9];
a1=-x2*D[f,x4]-x5*D[f,x7]-x9*D[f,x8];
a2=-x1*D[f,x2]-x3*D[f,x5]-x6*D[f,x9];
a3=(1/y)*x5*x8*I*f-(1/y)*x7*x9*I*f-x0*D[f,x1];
a4=-x1*D[f,x3]-x2*D[f,x5]-x4*D[f,x7];
a5=-x3*D[f,x6]-x5*D[f,x9]-x7*D[f,x8];
a6=-x1*D[f,x4]-x3*D[f,x7]-x6*D[f,x8];
a7=-(1/y)*x3*x8*I*f+(1/y)*x6*x7*I*f-x0*D[f,x2];
a8=-(1/y)*x2*x8*I*f+(1/y)*x4*x9*I*f-x0*D[f,x3];
a9=-x1*D[f,x6]-x2*D[f,x9]-x4*D[f,x8];
a10=(1/y)*x3*x9*I*f-(1/y)*x5*x6*I*f-x0*D[f,x4];
a11=(1/y)*x1*x8*I*f-(1/y)*x4*x6*I*f-x0*D[f,x5];
a12=(1/y)*x2*x7*I*f-(1/y)*x4*x5*I*f-x0*D[f,x6];
a13=-(1/y)*x1*x9*I*f+(1/y)*x2*x6*I*f-x0*D[f,x7];
a14=-(1/y)*x1*x7*I*f+(1/y)*x3*x4*I*f-x0*D[f,x9];
a15=(1/y)*x1*x5*I*f-(1/y)*x2*x3*I*f-x0*D[f,x8];
b0=y*D[f,x0];
b1=y*D[f,x1];
b2=y*D[f,x2];
b3=y*D[f,x3];
b4=y*D[f,x4];
b5=y*D[f,x5];
b6=y*D[f,x6];
b7=y*D[f,x7];
b8=y*D[f,x8];
b9=y*D[f,x9];
g0=I*x0*f;
g1=I*x1*f;
g2=I*x2*f;
g3=I*x3*f;
g4=I*x4*f;
g5=I*x5*f;
g6=I*x6*f;
g7=I*x7*f;
g8=I*x8*f;
g9=I*x9*f;
om=I*y*f;
am1=x4*D[f,x2]+x7*D[f,x5]+x8*D[f,x9];
am2=x2*D[f,x1]+x5*D[f,x3]+x9*D[f,x6];
am3=-y*D[D[I*f,x8],x5]+y*D[D[I*f,x9],x7]+x1*D[f,x0];
am4=x3*D[f,x1]+x5*D[f,x2]+x7*D[f,x4];
am5=x6*D[f,x3]+x8*D[f,x7]+x9*D[f,x5];
am6=x4*D[f,x1]+x7*D[f,x3]+x8*D[f,x6];
am7=y*D[D[I*f,x8],x3]-y*D[D[I*f,x7],x6]+x2*D[f,x0];
am8=y*D[D[I*f,x8],x2]-y*D[D[I*f,x9],x4]+x3*D[f,x0];
am9=x6*D[f,x1]+x8*D[f,x4]+x9*D[f,x2];
am10=-y*D[D[I*f,x9],x3]+y*D[D[I*f,x6],x5]+x4*D[f,x0];
am11=-y*D[D[I*f,x8],x1]+y*D[D[I*f,x6],x4]+x5*D[f,x0];
am12=-y*D[D[I*f,x7],x2]+y*D[D[I*f,x5],x4]+x6*D[f,x0];
am13=y*D[D[I*f,x9],x1]-y*D[D[I*f,x6],x2]+x7*D[f,x0];
am14=y*D[D[I*f,x7],x1]-y*D[D[I*f,x4],x3]+x9*D[f,x0];
am15=-y*D[D[I*f,x5],x1]+y*D[D[I*f,x3],x2]+x8*D[f,x0];
bm0=-(1/y)*(1/y)*x1*x5*x8*I*f+(1/y)*(1/y)*x1*x7*x9*I*f+(1/y)*(1/y)*x2*x3*x8*I*f-(1/y)*(1/y)*x2*x6*x7*I*f-(1/y)*(1/y)*x3*x4*x9*I*f+(1/y)*(1/y)*x4*x5*x6*I*f-x0*D[f,y];
bm1=2*(1/y)*x1*f+(1/y)*x1*x5*D[f,x5]+(1/y)*x1*x7*D[f,x7]+(1/y)*x1*x8*D[f,x8]+(1/y)*x1*x9*D[f,x9]-(1/y)*x2*x3*D[f,x5]-(1/y)*x2*x6*D[f,x9]-(1/y)*x3*x4*D[f,x7]-(1/y)*x4*x6*D[f,x8]+x0*D[D[I*f,x8],x5]-x0*D[D[I*f,x9],x7]+x1*D[f,y];
bm2=2-*f(1/y)*x1*x5*D[f,x3]-(1/y)*x1*x9*D[f,x6]+2*(1/y)*x2*f+(1/y)*x2*x3*D[f,x3]+(1/y)*x2*x6*D[f,x6]+(1/y)*x2*x7*D[f,x7]+(1/y)*x2*x8*D[f,x8]-(1/y)*x4*x5*D[f,x7]-(1/y)*x4*x9*D[f,x8]-x0*D[D[I*f,x8],x3]+x0*D[D[I*f,x7],x6]+x2*D[f,y];
bm3=-(1/y)*x1*x5*D[f,x2]-(1/y)*x1*x7*D[f,x4]+(1/y)*x2*x3*D[f,x2]+2*(1/y)*x3*f+(1/y)*x3*x4*D[f,x4]+(1/y)*x3*x8*D[f,x8]+(1/y)*x3*x9*D[f,x9]-(1/y)*x5*x6*D[f,x9]-(1/y)*x6*x7*D[f,x8]-x0*D[D[I*f,x8],x2]+x0*D[D[I*f,x9],x4]+x3*D[f,y];
bm4=-(1/y)*x1*x7*D[f,x3]-(1/y)*x1*x8*D[f,x6]-(1/y)*x2*x7*D[f,x5]-(1/y)*x2*x8*D[f,x9]+(1/y)*x3*x4*D[f,x3]+2*(1/y)*x4*f+(1/y)*x4*x5*D[f,x5]+(1/y)*x4*x6*D[f,x6]+(1/y)*x4*x9*D[f,x9]+x0*D[D[I*f,x9],x3]-x0*D[D[I*f,x6],x5]+x4*D[f,y];
bm5=(1/y)*x1*x5*D[f,x1]-(1/y)*x2*x3*D[f,x1]-(1/y)*x2*x7*D[f,x4]-(1/y)*x3*x9*D[f,x6]+(1/y)*x4*x5*D[f,x4]+2*(1/y)*x5*f+(1/y)*x5*x6*D[f,x6]+(1/y)*x5*x8*D[f,x8]-(1/y)*x7*x9*D[f,x8]+x0*D[D[I*f,x8],x1]-x0*D[D[I*f,x6],x4]+x5*D[f,y];
bm6=-(1/y)*x1*x8*D[f,x4]-(1/y)*x1*x9*D[f,x2]+(1/y)*x2*x6*D[f,x2]-(1/y)*x3*x8*D[f,x7]-(1/y)*x3*x9*D[f,x5]+(1/y)*x4*x6*D[f,x4]+(1/y)*x5*x6*D[f,x5]+2*(1/y)*x6*f+(1/y)*x6*x7*D[f,x7]+x0*D[D[I*f,x7],x2]-x0*D[D[I*f,x5],x4]+x6*D[f,y];
bm7=(1/y)*x1*x7*D[f,x1]+(1/y)*x2*x7*D[f,x2]-(1/y)*x3*x4*D[f,x1]-(1/y)*x3*x8*D[f,x6]-(1/y)*x4*x5*D[f,x2]-(1/y)*x5*x8*D[f,x9]+(1/y)*x6*x7*D[f,x6]+2*(1/y)*x7*f+(1/y)*x7*x9*D[f,x9]-x0*D[D[I*f,x9],x1]+x0*D[D[I*f,x6],x2]+x7*D[f,y];
bm8=(1/y)*x1*x8*D[f,x1]+(1/y)*x2*x8*D[f,x2]+(1/y)*x3*x8*D[f,x3]-(1/y)*x4*x6*D[f,x1]-(1/y)*x4*x9*D[f,x2]+(1/y)*x5*x8*D[f,x5]-(1/y)*x6*x7*D[f,x3]-(1/y)*x7*x9*D[f,x5]+2*(1/y)*x8*f+x0*D[D[I*f,x5],x1]-x0*D[D[I*f,x3],x2]+x8*D[f,y];
bm9=(1/y)*x1*x9*D[f,x1]-(1/y)*x2*x6*D[f,x1]-(1/y)*x2*x8*D[f,x4]+(1/y)*x3*x9*D[f,x3]+(1/y)*x4*x9*D[f,x4]-(1/y)*x5*x6*D[f,x3]-(1/y)*x5*x8*D[f,x7]+(1/y)*x7*x9*D[f,x7]+2*(1/y)*x9*f-x0*D[D[I*f,x7],x1]+x0*D[D[I*f,x4],x3]+x9*D[f,y];
gm0=-y*D[D[I*f,x0],y]-y*D[d5*D[f,x8],x1]+y*D[d7*D[f,x9],x1]+y*D[d3*D[f,x8],x2]-y*D[d6*D[f,x7],x2]-y*D[d4*D[f,x9],x3]+y*D[d5*D[f,x6],x4]-x0*D[D[I*f,x0],x0]-6*D[I*f,x0]-x1*D[D[I*f,x1],x0]-x2*D[D[I*f,x2],x0]-x3*D[D[I*f,x3],x0]-x4*D[D[I*f,x4],x0]-x5*D[D[I*f,x5],x0]-x6*D[D[I*f,x6],x0]-x7*D[D[I*f,x7],x0]-x8*D[D[I*f,x8],x0]-x9*D[D[I*f,x9],x0];
gm1=y*D[D[I*f,x1],y]+x0*D[D[I*f,x1],x0]+x1*D[D[I*f,x1],x1]+4*D[I*f,x1]+x2*D[D[I*f,x2],x1]+x3*D[D[I*f,x3],x1]+x4*D[D[I*f,x4],x1]+x5*D[D[I*f,x3],x2]+x5*x8*(1/y)*D[f,x0]+x6*D[D[I*f,x6],x1]+x7*D[D[I*f,x4],x3]-x7*x9*(1/y)*D[f,x0]+x8*D[D[I*f,x6],x4]+x9*D[D[I*f,x6],x2];
gm2=y*D[D[I*f,x2],y]+x0*D[D[I*f,x2],x0]+x1*D[D[I*f,x2],x1]+x2*D[D[I*f,x2],x2]+4*D[I*f,x2]+x3*D[D[I*f,x5],x1]-x3*x8*(1/y)*D[f,x0]+x4*D[D[I*f,x4],x2]+x5*D[D[I*f,x5],x2]+x6*D[D[I*f,x9],x1]+x6*x7*(1/y)*D[f,x0]+x7*D[D[I*f,x5],x4]+x8*D[D[I*f,x9],x4]+x9*D[D[I*f,x9],x2];
gm3=y*D[D[I*f,x3],y]+x0*D[D[I*f,x3],x0]+x1*D[D[I*f,x3],x1]+x2*D[D[I*f,x5],x1]-x2*x8*(1/y)*D[f,x0]+x3*D[D[I*f,x3],x3]+4*D[I*f,x3]+x4*D[D[I*f,x7],x1]+x4*x9*(1/y)*D[f,x0]+x5*D[D[I*f,x5],x3]+x6*D[D[I*f,x6],x3]+x7*D[D[I*f,x7],x3]+x8*D[D[I*f,x7],x6]+x9*D[D[I*f,x6],x5];
gm4=y*D[D[I*f,x4],y]+x0*D[D[I*f,x4],x0]+x1*D[D[I*f,x4],x1]+x2*D[D[I*f,x4],x2]+x3*D[D[I*f,x7],x1]+x3*x9*(1/y)*D[f,x0]+x4*D[D[I*f,x4],x4]+4*D[I*f,x4]+x5*D[D[I*f,x7],x2]-x5*x6*(1/y)*D[f,x0]+x6*D[D[I*f,x8],x1]+x7*D[D[I*f,x7],x4]+x8*D[D[I*f,x8],x4]+x9*D[D[I*f,x8],x2];
gm5=y*D[D[I*f,x5],y]+x0*D[D[I*f,x5],x0]+x1*D[D[I*f,x3],x2]+x1*x8*(1/y)*D[f,x0]+x2*D[D[I*f,x5],x2]+x3*D[D[I*f,x5],x3]+x4*D[D[I*f,x7],x2]-x4*x6*(1/y)*D[f,x0]+x5*D[D[I*f,x5],x5]+4*D[I*f,x5]+x6*D[D[I*f,x9],x3]+x7*D[D[I*f,x7],x5]+x8*D[D[I*f,x9],x7]+x9*D[D[I*f,x9],x5];
gm6=y*D[D[I*f,x6],y]+x0*D[D[I*f,x6],x0]+x1*D[D[I*f,x6],x1]+x2*D[D[I*f,x9],x1]+x2*x7*(1/y)*D[f,x0]+x3*D[D[I*f,x6],x3]+x4*D[D[I*f,x8],x1]-x4*x5*(1/y)*D[f,x0]+x5*D[D[I*f,x9],x3]+x6*D[D[I*f,x6],x6]+4*D[I*f,x6]+x7*D[D[I*f,x8],x3]+x8*D[D[I*f,x8],x6]+x9*D[D[I*f,x9],x6];
gm7=y*D[D[I*f,x7],y]+x0*D[D[I*f,x7],x0]+x1*D[D[I*f,x4],x3]-x1*x9*(1/y)*D[f,x0]+x2*D[D[I*f,x5],x4]+x2*x6*(1/y)*D[f,x0]+x3*D[D[I*f,x7],x3]+x4*D[D[I*f,x7],x4]+x5*D[D[I*f,x7],x5]+x6*D[D[I*f,x8],x3]+x7*D[D[I*f,x7],x7]+4*D[I*f,x7]+x8*D[D[I*f,x8],x7]+x9*D[D[I*f,x8],x5];
gm8=y*D[D[I*f,x8],y]+x0*D[D[I*f,x8],x0]+x1*D[D[I*f,x6],x4]+x1*x5*(1/y)*D[f,x0]-x2*x3*(1/y)*D[f,x0]+x2*D[D[I*f,x9],x4]+x3*D[D[I*f,x7],x6]+x4*D[D[I*f,x8],x4]+x5*D[D[I*f,x9],x7]+x6*D[D[I*f,x8],x6]+x7*D[D[I*f,x8],x7]+x8*D[D[I*f,x8],x8]+4*D[I*f,x8]+x9*D[D[I*f,x9],x8];
gm9=y*D[D[I*f,x9],y]+x0*D[D[I*f,x9],x0]+x1*D[D[I*f,x6],x2]-x1*x7*(1/y)*D[f,x0]+x2*D[D[I*f,x9],x2]+x3*x4*(1/y)*D[f,x0]+x3*D[D[I*f,x6],x5]+x4*D[D[I*f,x8],x2]+x5*D[D[I*f,x9],x5]+x6*D[D[I*f,x9],x6]+x7*D[D[I*f,x8],x5]+x8*D[D[I*f,x9],x8]+x9*D[D[I*f,x9],x9]+4*D[I*f,x9];
mom=-y*D[D[I*f,y],y]-6*(1/y)*I*f-6*D[I*f,y]-x0*D[D[I*f,x0],y]+x0*D[d5*D[f,x8],x1]-x0*D[d7*D[f,x9],x1]-x0*D[d3*D[f,x8],x2]+x0*D[d6*D[f,x7],x2]+x0*D[d4*D[f,x9],x3]-x0*D[d5*D[f,x6],x4]-2*x1*(1/y)*D[I*f,x1]-x1*D[D[I*f,x1],y]-x1*x5*(1/y)*D[D[I*f,x5],x1]+x1*x5*(1/y)*D[D[I*f,x3],x2]+x1*x5*x8*(1/y)*(1/y)*D[f,x0]-x1*x7*(1/y)*D[D[I*f,x7],x1]+x1*x7*(1/y)*D[D[I*f,x4],x3]-x1*x7*x9*(1/y)*(1/y)*D[f,x0]-x1*x8*(1/y)*D[D[I*f,x8],x1]+x1*x8*(1/y)*D[D[I*f,x6],x4]-x1*x9*(1/y)*D[D[I*f,x9],x1]+x1*x9*(1/y)*D[D[I*f,x6],x2]-2*x2*(1/y)*D[I*f,x2]-x2*D[D[I*f,x2],y]+x2*x3*(1/y)*D[D[I*f,x5],x1]-x2*x3*(1/y)*D[D[I*f,x3],x2]-x2*x3*x8*(1/y)*(1/y)*D[f,x0]+x2*x6*(1/y)*D[D[I*f,x9],x1]-x2*x6*(1/y)*D[D[I*f,x6],x2]+x2*x6*x7*(1/y)*(1/y)*D[f,x0]-x2*x7*(1/y)*D[D[I*f,x7],x2]+x2*x7*(1/y)*D[D[I*f,x5],x4]-x2*x8*(1/y)*D[D[I*f,x8],x2]+x2*x8*(1/y)*D[D[I*f,x9],x4]-2*x3*(1/y)*D[I*f,x3]-x3*D[D[I*f,x3],y]+x3*x4*(1/y)*D[D[I*f,x7],x1]-x3*x4*(1/y)*D[D[I*f,x4],x3]+x3*x4*x9*(1/y)*(1/y)*D[f,x0]-x3*x8*(1/y)*D[D[I*f,x8],x3]+x3*x8*(1/y)*D[D[I*f,x7],x6]-x3*x9*(1/y)*D[D[I*f,x9],x3]+x3*x9*(1/y)*D[D[I*f,x6],x5]-2*x4*(1/y)*D[I*f,x4]-x4*D[D[I*f,x4],y]+x4*x5*(1/y)*D[D[I*f,x7],x2]-x4*x5*(1/y)*D[D[I*f,x5],x4]-x4*x5*x6*(1/y)*(1/y)*D[f,x0]+x4*x6*(1/y)*D[D[I*f,x8],x1]-x4*x6*(1/y)*D[D[I*f,x6],x4]+x4*x9*(1/y)*D[D[I*f,x8],x2]-x4*x9*(1/y)*D[D[I*f,x9],x4]-2*x5*(1/y)*D[I*f,x5]-x5*D[D[I*f,x5],y]+x5*x6*(1/y)*D[D[I*f,x9],x3]-x5*x6*(1/y)*D[D[I*f,x6],x5]-x5*x8*(1/y)*D[D[I*f,x8],x5]+x5*x8*(1/y)*D[D[I*f,x9],x7]-2*x6*(1/y)*D[I*f,x6]-x6*D[D[I*f,x6],y]+x6*x7*(1/y)*D[D[I*f,x8],x3]-x6*x7*(1/y)*D[D[I*f,x7],x6]-2*x7*(1/y)*D[I*f,x7]-x7*D[D[I*f,x7],y]+x7*x9*(1/y)*D[D[I*f,x8],x5]-x7*x9*(1/y)*D[D[I*f,x9],x7]-2*x8*(1/y)*D[I*f,x8]-x8*D[D[I*f,x8],y]-2*x9*(1/y)*D[I*f,x9]-x9*D[D[I*f,x9],y];