H0=-y*D[f,y]+x0*D[f,x0];
H1=-1*f-x0*D[f,x0]+x1*D[f,x1]-x2*D[f,x2]-x3*D[f,x3];
H2=-1*f-x0*D[f,x0]-x1*D[f,x1]+x2*D[f,x2]-x3*D[f,x3];
H3=-1*f-x0*D[f,x0]-x1*D[f,x1]-x2*D[f,x2]+x3*D[f,x3];
a1=-(1/y)*x2*x3*I*f-x0*D[f,x1];
a2=-(1/y)*x1*x3*I*f-x0*D[f,x2];
a3=-(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];
g0=I*x0*f;
g1=I*x1*f;
g2=I*x2*f;
g3=I*x3*f;
om=I*y*f;
am1=y*D[D[I*f,x3],x2]+x1*D[f,x0];
am2=y*D[D[I*f,x3],x1]+x2*D[f,x0];
am3=y*D[D[I*f,x2],x1]+x3*D[f,x0];
bm0=(1/y)*(1/y)*x1*x2*x3*I*f-x0*D[f,y];
bm1=(1/y)*x1*f+(1/y)*x1*x2*D[f,x2]+(1/y)*x1*x3*D[f,x3]-x0*D[D[I*f,x3],x2]+x1*D[f,y];
bm2=(1/y)*x1*x2*D[f,x1]+(1/y)*x2*f+(1/y)*x2*x3*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-x0*D[D[I*f,x2],x1]+x3*D[f,y];
gm0=-y*D[D[I*f,x0],y]+y*D[d2*D[f,x3],x1]-x0*D[D[I*f,x0],x0]-3*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];
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];
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]+2*D[I*f,x2];
gm3=y*D[D[I*f,x3],y]+x0*D[D[I*f,x3],x0]-x1*x2*(1/y)*D[f,x0]+x3*D[D[I*f,x3],x3]+2*D[I*f,x3];
mom=-y*D[D[I*f,y],y]-(1/y)*I*f-(1/y)*(1/y)*x1*x2*x3*D[f,x0]-(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)*x2*D[I*f,x2]-(1/y)*x2*x3*D[D[I*f,x3],x2]-(1/y)*x3*D[I*f,x3]-3*D[I*f,y]-x0*D[D[I*f,x0],y]-x0*D[d2*D[f,x3],x1]-x1*D[D[I*f,x1],y]-x2*D[D[I*f,x2],y]-x3*D[D[I*f,x3],y];