H0=-y*D[f,y]+x0*D[f,x0];
H1=-2*f-x0*D[f,x0]+x1*D[f,x1]-x2*D[f,x2]-x3*D[f,x3]-x4*D[f,x4]-x5*D[f,x5];
H2=-1*f-x0*D[f,x0]-x1*D[f,x1]+x2*D[f,x2]-x3*D[f,x3];
H3=-x2*D[f,x2]+x3*D[f,x3]+x4*D[f,x4]-x5*D[f,x5];
H4=-x2*D[f,x2]+x3*D[f,x3]-x4*D[f,x4]+x5*D[f,x5];
a1=-(1/y)*x2*x3*I*f+(1/y)*x4*x5*I*f-x0*D[f,x1];
a2=-(1/y)*x1*x3*I*f-x0*D[f,x2];
a3=x2*D[f,x4]+x5*D[f,x3];
a4=-x2*D[f,x5]-x4*D[f,x3];
a5=(1/y)*x1*x5*I*f-x0*D[f,x4];
a6=(1/y)*x1*x4*I*f-x0*D[f,x5];
a7=-(1/y)*x1*x2*I*f-x0*D[f,x3];
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];
g0=I*x0*f;
g1=I*x1*f;
g2=I*x2*f;
g3=I*x3*f;
g4=I*x4*f;
g5=I*x5*f;
om=I*y*f;
am1=y*D[D[I*f,x3],x2]-y*D[D[I*f,x5],x4]+x1*D[f,x0];
am2=y*D[D[I*f,x3],x1]+x2*D[f,x0];
am3=-x3*D[f,x5]-x4*D[f,x2];
am4=x3*D[f,x4]+x5*D[f,x2];
am5=-y*D[D[I*f,x5],x1]+x4*D[f,x0];
am6=-y*D[D[I*f,x4],x1]+x5*D[f,x0];
am7=y*D[D[I*f,x2],x1]+x3*D[f,x0];
bm0=(1/y)*(1/y)*x1*x2*x3*I*f-(1/y)*(1/y)*x1*x4*x5*I*f-x0*D[f,y];
bm1=2*(1/y)*x1*f+(1/y)*x1*x2*D[f,x2]+(1/y)*x1*x3*D[f,x3]+(1/y)*x1*x4*D[f,x4]+(1/y)*x1*x5*D[f,x5]-x0*D[D[I*f,x3],x2]+x0*D[D[I*f,x5],x4]+x1*D[f,y];
bm2=(1/y)*x1*x2*D[f,x1]+(1/y)*x2*f+(1/y)*x2*x3*D[f,x3]-(1/y)*x4*x5*D[f,x3]-x0*D[D[I*f,x3],x1]+x2*D[f,y];
bm3=(1/y)*x1*x3*D[f,x1]+(1/y)*x2*x3*D[f,x2]+(1/y)*x3*f-(1/y)*x4*x5*D[f,x2]-x0*D[D[I*f,x2],x1]+x3*D[f,y];
bm4=(1/y)*x1*x4*D[f,x1]-(1/y)*x2*x3*D[f,x5]+(1/y)*x4*f+(1/y)*x4*x5*D[f,x5]+x0*D[D[I*f,x5],x1]+x4*D[f,y];
bm5=(1/y)*x1*x5*D[f,x1]-(1/y)*x2*x3*D[f,x4]+(1/y)*x4*x5*D[f,x4]+(1/y)*x5*f+x0*D[D[I*f,x4],x1]+x5*D[f,y];
gm0=-y*D[D[I*f,x0],y]+y*D[d2*D[f,x3],x1]-y*D[d4*D[f,x5],x1]-x0*D[D[I*f,x0],x0]-4*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];
gm1=y*D[D[I*f,x1],y]+x0*D[D[I*f,x1],x0]+x1*D[D[I*f,x1],x1]+2*D[I*f,x1]-x2*x3*(1/y)*D[f,x0]+x4*x5*(1/y)*D[f,x0];
gm2=y*D[D[I*f,x2],y]+x0*D[D[I*f,x2],x0]-x1*x3*(1/y)*D[f,x0]+x2*D[D[I*f,x2],x2]+3*D[I*f,x2]+x3*D[D[I*f,x5],x4]+x4*D[D[I*f,x4],x2]+x5*D[D[I*f,x5],x2];
gm3=y*D[D[I*f,x3],y]+x0*D[D[I*f,x3],x0]-x1*x2*(1/y)*D[f,x0]+x2*D[D[I*f,x5],x4]+x3*D[D[I*f,x3],x3]+3*D[I*f,x3]+x4*D[D[I*f,x4],x3]+x5*D[D[I*f,x5],x3];
gm4=y*D[D[I*f,x4],y]+x0*D[D[I*f,x4],x0]+x1*x5*(1/y)*D[f,x0]+x2*D[D[I*f,x4],x2]+x3*D[D[I*f,x4],x3]+x4*D[D[I*f,x4],x4]+3*D[I*f,x4]+x5*D[D[I*f,x3],x2];
gm5=y*D[D[I*f,x5],y]+x0*D[D[I*f,x5],x0]+x1*x4*(1/y)*D[f,x0]+x2*D[D[I*f,x5],x2]+x3*D[D[I*f,x5],x3]+x4*D[D[I*f,x3],x2]+x5*D[D[I*f,x5],x5]+3*D[I*f,x5];
mom=-y*D[D[I*f,y],y]-2*(1/y)*I*f-(1/y)*(1/y)*x1*x2*x3*D[f,x0]+(1/y)*(1/y)*x1*x4*x5*D[f,x0]-2*(1/y)*x1*D[I*f,x1]-(1/y)*x1*x2*D[D[I*f,x2],x1]-(1/y)*x1*x3*D[D[I*f,x3],x1]-(1/y)*x1*x4*D[D[I*f,x4],x1]-(1/y)*x1*x5*D[D[I*f,x5],x1]-(1/y)*x2*D[I*f,x2]-(1/y)*x2*x3*D[D[I*f,x3],x2]+(1/y)*x2*x3*D[D[I*f,x5],x4]-(1/y)*x3*D[I*f,x3]-(1/y)*x4*D[I*f,x4]+(1/y)*x4*x5*D[D[I*f,x3],x2]-(1/y)*x4*x5*D[D[I*f,x5],x4]-(1/y)*x5*D[I*f,x5]-4*D[I*f,y]-x0*D[D[I*f,x0],y]-x0*D[d2*D[f,x3],x1]+x0*D[d4*D[f,x5],x1]-x1*D[D[I*f,x1],y]-x2*D[D[I*f,x2],y]-x3*D[D[I*f,x3],y]-x4*D[D[I*f,x4],y]-x5*D[D[I*f,x5],y];