H0=-y*D[f,y]+x0*D[f,x0];
H1=-x2*D[f,x2]+x3*D[f,x3]-x4*D[f,x4]+x5*D[f,x5]-x6*D[f,x6]+x8*D[f,x8]-x9*D[f,x9]+x12*D[f,x12];
H2=-3*f-x0*D[f,x0]+x1*D[f,x1]-x7*D[f,x7]-x10*D[f,x10]-x11*D[f,x11]-x13*D[f,x13]-x14*D[f,x14]-x15*D[f,x15];
H3=-x1*D[f,x1]+x2*D[f,x2]-x5*D[f,x5]+x7*D[f,x7]-x8*D[f,x8]+x11*D[f,x11]-x12*D[f,x12]+x15*D[f,x15];
H4=-x2*D[f,x2]-x3*D[f,x3]+x4*D[f,x4]+x5*D[f,x5]-x11*D[f,x11]+x13*D[f,x13]+x14*D[f,x14]-x15*D[f,x15];
H5=-x4*D[f,x4]-x5*D[f,x5]+x6*D[f,x6]-x7*D[f,x7]+x8*D[f,x8]+x10*D[f,x10]+x11*D[f,x11]-x13*D[f,x13];
H6=-x6*D[f,x6]-x8*D[f,x8]+x9*D[f,x9]-x11*D[f,x11]+x12*D[f,x12]+x13*D[f,x13]-x14*D[f,x14]+x15*D[f,x15];
Hom=-9*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]-x10*D[f,x10]-x11*D[f,x11]-x12*D[f,x12]-x13*D[f,x13]-x14*D[f,x14]-x15*D[f,x15];
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];
b10=y*D[f,x10];
b11=y*D[f,x11];
b12=y*D[f,x12];
b13=y*D[f,x13];
b14=y*D[f,x14];
b15=y*D[f,x15];
g0=x0*I*f;
g1=x1*I*f;
g2=x2*I*f;
g3=x3*I*f;
g4=x4*I*f;
g5=x5*I*f;
g6=x6*I*f;
g7=x7*I*f;
g8=x8*I*f;
g9=x9*I*f;
g10=x10*I*f;
g11=x11*I*f;
g12=x12*I*f;
g13=x13*I*f;
g14=x14*I*f;
g15=x15*I*f;
om=y*I*f;
a2=(1/y)*x7*x10*I*f-(1/y)*x11*x13*I*f+(1/y)*x14*x15*I*f-x0*D[f,x1];
a8=-(1/y)*x5*x10*I*f+(1/y)*x8*x13*I*f-(1/y)*x12*x14*I*f-x0*D[f,x2];
a12=(1/y)*x4*x10*I*f-(1/y)*x6*x13*I*f+(1/y)*x9*x14*I*f-x0*D[f,x3];
a14=(1/y)*x3*x10*I*f-(1/y)*x8*x15*I*f+(1/y)*x11*x12*I*f-x0*D[f,x4];
a17=-(1/y)*x2*x10*I*f+(1/y)*x6*x15*I*f-(1/y)*x9*x11*I*f-x0*D[f,x5];
a19=-(1/y)*x3*x13*I*f+(1/y)*x5*x15*I*f-(1/y)*x7*x12*I*f-x0*D[f,x6];
a21=(1/y)*x1*x10*I*f-(1/y)*x6*x12*I*f+(1/y)*x8*x9*I*f-x0*D[f,x7];
a22=(1/y)*x2*x13*I*f-(1/y)*x4*x15*I*f+(1/y)*x7*x9*I*f-x0*D[f,x8];
a24=(1/y)*x3*x14*I*f-(1/y)*x5*x11*I*f+(1/y)*x7*x8*I*f-x0*D[f,x9];
a30=(1/y)*x1*x7*I*f-(1/y)*x2*x5*I*f+(1/y)*x3*x4*I*f-x0*D[f,x10];
a25=-(1/y)*x1*x13*I*f+(1/y)*x4*x12*I*f-(1/y)*x5*x9*I*f-x0*D[f,x11];
a26=-(1/y)*x2*x14*I*f+(1/y)*x4*x11*I*f-(1/y)*x6*x7*I*f-x0*D[f,x12];
a29=-(1/y)*x1*x11*I*f+(1/y)*x2*x8*I*f-(1/y)*x3*x6*I*f-x0*D[f,x13];
a27=(1/y)*x1*x15*I*f-(1/y)*x2*x12*I*f+(1/y)*x3*x9*I*f-x0*D[f,x14];
a28=(1/y)*x1*x14*I*f-(1/y)*x4*x8*I*f+(1/y)*x5*x6*I*f-x0*D[f,x15];
a1=x2*D[f,x3]+x4*D[f,x5]+x6*D[f,x8]+x9*D[f,x12];
a3=x1*D[f,x2]+x5*D[f,x7]+x8*D[f,x11]+x12*D[f,x15];
a4=x2*D[f,x4]+x3*D[f,x5]+x11*D[f,x14]+x15*D[f,x13];
a5=x4*D[f,x6]+x5*D[f,x8]+x7*D[f,x11]+x13*D[f,x10];
a6=x6*D[f,x9]+x8*D[f,x12]+x11*D[f,x15]+x14*D[f,x13];
a7=x1*D[f,x3]-x4*D[f,x7]-x6*D[f,x11]-x9*D[f,x15];
a9=x1*D[f,x4]-x3*D[f,x7]+x8*D[f,x14]+x12*D[f,x13];
a10=x2*D[f,x6]+x3*D[f,x8]-x7*D[f,x14]+x15*D[f,x10];
a11=-x4*D[f,x9]-x5*D[f,x12]-x7*D[f,x15]+x14*D[f,x10];
a13=x1*D[f,x5]+x2*D[f,x7]-x6*D[f,x14]-x9*D[f,x13];
a15=x1*D[f,x6]-x3*D[f,x11]-x5*D[f,x14]+x12*D[f,x10];
a16=-x2*D[f,x9]-x3*D[f,x12]+x7*D[f,x13]+x11*D[f,x10];
a18=-x1*D[f,x8]-x2*D[f,x11]-x4*D[f,x14]+x9*D[f,x10];
bm0=-(1/y)*(1/y)*x1*x7*x10*I*f+(1/y)*(1/y)*x1*x11*x13*I*f-(1/y)*(1/y)*x1*x14*x15*I*f+(1/y)*(1/y)*x2*x5*x10*I*f-(1/y)*(1/y)*x2*x8*x13*I*f+(1/y)*(1/y)*x2*x12*x14*I*f-(1/y)*(1/y)*x3*x4*x10*I*f+(1/y)*(1/y)*x3*x6*x13*I*f-(1/y)*(1/y)*x3*x9*x14*I*f+(1/y)*(1/y)*x4*x8*x15*I*f-(1/y)*(1/y)*x4*x11*x12*I*f-(1/y)*(1/y)*x5*x6*x15*I*f+(1/y)*(1/y)*x5*x9*x11*I*f+(1/y)*(1/y)*x6*x7*x12*I*f-(1/y)*(1/y)*x7*x8*x9*I*f-x0*D[f,y];
am2=-y*D[D[I*f,x10],x7]+y*D[D[I*f,x13],x11]-y*D[D[I*f,x15],x14]+x1*D[f,x0];
am8=y*D[D[I*f,x10],x5]-y*D[D[I*f,x13],x8]+y*D[D[I*f,x14],x12]+x2*D[f,x0];
am12=-y*D[D[I*f,x10],x4]+y*D[D[I*f,x13],x6]-y*D[D[I*f,x14],x9]+x3*D[f,x0];
am14=-y*D[D[I*f,x10],x3]+y*D[D[I*f,x15],x8]-y*D[D[I*f,x12],x11]+x4*D[f,x0];
am17=y*D[D[I*f,x10],x2]-y*D[D[I*f,x15],x6]+y*D[D[I*f,x11],x9]+x5*D[f,x0];
am19=y*D[D[I*f,x13],x3]-y*D[D[I*f,x15],x5]+y*D[D[I*f,x12],x7]+x6*D[f,x0];
am21=-y*D[D[I*f,x10],x1]+y*D[D[I*f,x12],x6]-y*D[D[I*f,x9],x8]+x7*D[f,x0];
am22=-y*D[D[I*f,x13],x2]+y*D[D[I*f,x15],x4]-y*D[D[I*f,x9],x7]+x8*D[f,x0];
am24=-y*D[D[I*f,x14],x3]+y*D[D[I*f,x11],x5]-y*D[D[I*f,x8],x7]+x9*D[f,x0];
am30=-y*D[D[I*f,x7],x1]+y*D[D[I*f,x5],x2]-y*D[D[I*f,x4],x3]+x10*D[f,x0];
am25=y*D[D[I*f,x13],x1]-y*D[D[I*f,x12],x4]+y*D[D[I*f,x9],x5]+x11*D[f,x0];
am26=y*D[D[I*f,x14],x2]-y*D[D[I*f,x11],x4]+y*D[D[I*f,x7],x6]+x12*D[f,x0];
am29=y*D[D[I*f,x11],x1]-y*D[D[I*f,x8],x2]+y*D[D[I*f,x6],x3]+x13*D[f,x0];
am27=-y*D[D[I*f,x15],x1]+y*D[D[I*f,x12],x2]-y*D[D[I*f,x9],x3]+x14*D[f,x0];
am28=-y*D[D[I*f,x14],x1]+y*D[D[I*f,x8],x4]-y*D[D[I*f,x6],x5]+x15*D[f,x0];
am1=-x3*D[f,x2]-x5*D[f,x4]-x8*D[f,x6]-x12*D[f,x9];
am3=-x2*D[f,x1]-x7*D[f,x5]-x11*D[f,x8]-x15*D[f,x12];
am4=-x4*D[f,x2]-x5*D[f,x3]-x13*D[f,x15]-x14*D[f,x11];
am5=-x6*D[f,x4]-x8*D[f,x5]-x10*D[f,x13]-x11*D[f,x7];
am6=-x9*D[f,x6]-x12*D[f,x8]-x13*D[f,x14]-x15*D[f,x11];
am7=-x3*D[f,x1]+x7*D[f,x4]+x11*D[f,x6]+x15*D[f,x9];
am9=-x4*D[f,x1]+x7*D[f,x3]-x13*D[f,x12]-x14*D[f,x8];
am10=-x6*D[f,x2]-x8*D[f,x3]-x10*D[f,x15]+x14*D[f,x7];
am11=x9*D[f,x4]-x10*D[f,x14]+x12*D[f,x5]+x15*D[f,x7];
am13=-x5*D[f,x1]-x7*D[f,x2]+x13*D[f,x9]+x14*D[f,x6];
am15=-x6*D[f,x1]-x10*D[f,x12]+x11*D[f,x3]+x14*D[f,x5];
am16=x9*D[f,x2]-x10*D[f,x11]+x12*D[f,x3]-x13*D[f,x7];
am18=x8*D[f,x1]-x10*D[f,x9]+x11*D[f,x2]+x14*D[f,x4];
gm0=-y*D[D[I*f,x0],y]-y*D[d7*D[f,x10],x1]+y*D[d11*D[f,x13],x1]-y*D[d14*D[f,x15],x1]+y*D[d5*D[f,x10],x2]-y*D[d8*D[f,x13],x2]+y*D[d12*D[f,x14],x2]-y*D[d4*D[f,x10],x3]+y*D[d6*D[f,x13],x3]-y*D[d9*D[f,x14],x3]+y*D[d8*D[f,x15],x4]-y*D[d11*D[f,x12],x4]-y*D[d6*D[f,x15],x5]+y*D[d9*D[f,x11],x5]+y*D[d7*D[f,x12],x6]-y*D[d8*D[f,x9],x7]-x0*D[D[I*f,x0],x0]-9*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]-x10*D[D[I*f,x10],x0]-x11*D[D[I*f,x11],x0]-x12*D[D[I*f,x12],x0]-x13*D[D[I*f,x13],x0]-x14*D[D[I*f,x14],x0]-x15*D[D[I*f,x15],x0];
bm1=3*(1/y)*x1*f+(1/y)*x1*x7*D[f,x7]+(1/y)*x1*x10*D[f,x10]+(1/y)*x1*x11*D[f,x11]+(1/y)*x1*x13*D[f,x13]+(1/y)*x1*x14*D[f,x14]+(1/y)*x1*x15*D[f,x15]-(1/y)*x2*x5*D[f,x7]-(1/y)*x2*x8*D[f,x11]-(1/y)*x2*x12*D[f,x15]+(1/y)*x3*x4*D[f,x7]+(1/y)*x3*x6*D[f,x11]+(1/y)*x3*x9*D[f,x15]-(1/y)*x4*x8*D[f,x14]-(1/y)*x4*x12*D[f,x13]+(1/y)*x5*x6*D[f,x14]+(1/y)*x5*x9*D[f,x13]-(1/y)*x6*x12*D[f,x10]+(1/y)*x8*x9*D[f,x10]+x0*D[D[I*f,x10],x7]-x0*D[D[I*f,x13],x11]+x0*D[D[I*f,x15],x14]+x1*D[f,y];
bm2=-(1/y)*x1*x7*D[f,x5]-(1/y)*x1*x11*D[f,x8]-(1/y)*x1*x15*D[f,x12]+3*(1/y)*x2*f+(1/y)*x2*x5*D[f,x5]+(1/y)*x2*x8*D[f,x8]+(1/y)*x2*x10*D[f,x10]+(1/y)*x2*x12*D[f,x12]+(1/y)*x2*x13*D[f,x13]+(1/y)*x2*x14*D[f,x14]-(1/y)*x3*x4*D[f,x5]-(1/y)*x3*x6*D[f,x8]-(1/y)*x3*x9*D[f,x12]-(1/y)*x4*x11*D[f,x14]-(1/y)*x4*x15*D[f,x13]+(1/y)*x6*x7*D[f,x14]-(1/y)*x6*x15*D[f,x10]+(1/y)*x7*x9*D[f,x13]+(1/y)*x9*x11*D[f,x10]-x0*D[D[I*f,x10],x5]+x0*D[D[I*f,x13],x8]-x0*D[D[I*f,x14],x12]+x2*D[f,y];
bm3=(1/y)*x1*x7*D[f,x4]+(1/y)*x1*x11*D[f,x6]+(1/y)*x1*x15*D[f,x9]-(1/y)*x2*x5*D[f,x4]-(1/y)*x2*x8*D[f,x6]-(1/y)*x2*x12*D[f,x9]+3*(1/y)*x3*f+(1/y)*x3*x4*D[f,x4]+(1/y)*x3*x6*D[f,x6]+(1/y)*x3*x9*D[f,x9]+(1/y)*x3*x10*D[f,x10]+(1/y)*x3*x13*D[f,x13]+(1/y)*x3*x14*D[f,x14]-(1/y)*x5*x11*D[f,x14]-(1/y)*x5*x15*D[f,x13]+(1/y)*x7*x8*D[f,x14]+(1/y)*x7*x12*D[f,x13]-(1/y)*x8*x15*D[f,x10]+(1/y)*x11*x12*D[f,x10]+x0*D[D[I*f,x10],x4]-x0*D[D[I*f,x13],x6]+x0*D[D[I*f,x14],x9]+x3*D[f,y];
bm4=(1/y)*x1*x7*D[f,x3]-(1/y)*x1*x13*D[f,x12]-(1/y)*x1*x14*D[f,x8]-(1/y)*x2*x5*D[f,x3]-(1/y)*x2*x13*D[f,x15]-(1/y)*x2*x14*D[f,x11]+(1/y)*x3*x4*D[f,x3]+3*(1/y)*x4*f+(1/y)*x4*x8*D[f,x8]+(1/y)*x4*x10*D[f,x10]+(1/y)*x4*x11*D[f,x11]+(1/y)*x4*x12*D[f,x12]+(1/y)*x4*x15*D[f,x15]-(1/y)*x5*x6*D[f,x8]-(1/y)*x5*x9*D[f,x12]-(1/y)*x6*x7*D[f,x11]-(1/y)*x6*x13*D[f,x10]-(1/y)*x7*x9*D[f,x15]+(1/y)*x9*x14*D[f,x10]+x0*D[D[I*f,x10],x3]-x0*D[D[I*f,x15],x8]+x0*D[D[I*f,x12],x11]+x4*D[f,y];
bm5=-(1/y)*x1*x7*D[f,x2]+(1/y)*x1*x13*D[f,x9]+(1/y)*x1*x14*D[f,x6]+(1/y)*x2*x5*D[f,x2]-(1/y)*x3*x4*D[f,x2]-(1/y)*x3*x13*D[f,x15]-(1/y)*x3*x14*D[f,x11]-(1/y)*x4*x8*D[f,x6]-(1/y)*x4*x12*D[f,x9]+3*(1/y)*x5*f+(1/y)*x5*x6*D[f,x6]+(1/y)*x5*x9*D[f,x9]+(1/y)*x5*x10*D[f,x10]+(1/y)*x5*x11*D[f,x11]+(1/y)*x5*x15*D[f,x15]-(1/y)*x7*x8*D[f,x11]-(1/y)*x7*x12*D[f,x15]-(1/y)*x8*x13*D[f,x10]+(1/y)*x12*x14*D[f,x10]-x0*D[D[I*f,x10],x2]+x0*D[D[I*f,x15],x6]-x0*D[D[I*f,x11],x9]+x5*D[f,y];
bm6=-(1/y)*x1*x10*D[f,x12]+(1/y)*x1*x11*D[f,x3]+(1/y)*x1*x14*D[f,x5]-(1/y)*x2*x8*D[f,x3]-(1/y)*x2*x10*D[f,x15]+(1/y)*x2*x14*D[f,x7]+(1/y)*x3*x6*D[f,x3]-(1/y)*x4*x8*D[f,x5]-(1/y)*x4*x10*D[f,x13]-(1/y)*x4*x11*D[f,x7]+(1/y)*x5*x6*D[f,x5]+3*(1/y)*x6*f+(1/y)*x6*x7*D[f,x7]+(1/y)*x6*x12*D[f,x12]+(1/y)*x6*x13*D[f,x13]+(1/y)*x6*x15*D[f,x15]-(1/y)*x8*x9*D[f,x12]-(1/y)*x9*x11*D[f,x15]-(1/y)*x9*x14*D[f,x13]-x0*D[D[I*f,x13],x3]+x0*D[D[I*f,x15],x5]-x0*D[D[I*f,x12],x7]+x6*D[f,y];
bm7=(1/y)*x1*x7*D[f,x1]-(1/y)*x2*x5*D[f,x1]+(1/y)*x2*x13*D[f,x9]+(1/y)*x2*x14*D[f,x6]+(1/y)*x3*x4*D[f,x1]+(1/y)*x3*x13*D[f,x12]+(1/y)*x3*x14*D[f,x8]-(1/y)*x4*x11*D[f,x6]-(1/y)*x4*x15*D[f,x9]-(1/y)*x5*x11*D[f,x8]-(1/y)*x5*x15*D[f,x12]+(1/y)*x6*x7*D[f,x6]+3*(1/y)*x7*f+(1/y)*x7*x8*D[f,x8]+(1/y)*x7*x9*D[f,x9]+(1/y)*x7*x10*D[f,x10]+(1/y)*x7*x12*D[f,x12]-(1/y)*x11*x13*D[f,x10]+(1/y)*x14*x15*D[f,x10]+x0*D[D[I*f,x10],x1]-x0*D[D[I*f,x12],x6]+x0*D[D[I*f,x9],x8]+x7*D[f,y];
bm8=(1/y)*x1*x10*D[f,x9]-(1/y)*x1*x11*D[f,x2]-(1/y)*x1*x14*D[f,x4]+(1/y)*x2*x8*D[f,x2]-(1/y)*x3*x6*D[f,x2]-(1/y)*x3*x10*D[f,x15]+(1/y)*x3*x14*D[f,x7]+(1/y)*x4*x8*D[f,x4]-(1/y)*x5*x6*D[f,x4]-(1/y)*x5*x10*D[f,x13]-(1/y)*x5*x11*D[f,x7]-(1/y)*x6*x12*D[f,x9]+(1/y)*x7*x8*D[f,x7]+3*(1/y)*x8*f+(1/y)*x8*x9*D[f,x9]+(1/y)*x8*x13*D[f,x13]+(1/y)*x8*x15*D[f,x15]-(1/y)*x11*x12*D[f,x15]-(1/y)*x12*x14*D[f,x13]+x0*D[D[I*f,x13],x2]-x0*D[D[I*f,x15],x4]+x0*D[D[I*f,x9],x7]+x8*D[f,y];
bm9=(1/y)*x1*x10*D[f,x8]+(1/y)*x1*x13*D[f,x5]+(1/y)*x1*x15*D[f,x3]+(1/y)*x2*x10*D[f,x11]-(1/y)*x2*x12*D[f,x3]+(1/y)*x2*x13*D[f,x7]+(1/y)*x3*x9*D[f,x3]+(1/y)*x4*x10*D[f,x14]-(1/y)*x4*x12*D[f,x5]-(1/y)*x4*x15*D[f,x7]+(1/y)*x5*x9*D[f,x5]-(1/y)*x6*x12*D[f,x8]-(1/y)*x6*x13*D[f,x14]-(1/y)*x6*x15*D[f,x11]+(1/y)*x7*x9*D[f,x7]+(1/y)*x8*x9*D[f,x8]+3*(1/y)*x9*f+(1/y)*x9*x11*D[f,x11]+(1/y)*x9*x14*D[f,x14]+x0*D[D[I*f,x14],x3]-x0*D[D[I*f,x11],x5]+x0*D[D[I*f,x8],x7]+x9*D[f,y];
bm10=(1/y)*x1*x10*D[f,x1]+(1/y)*x2*x10*D[f,x2]+(1/y)*x3*x10*D[f,x3]+(1/y)*x4*x10*D[f,x4]+(1/y)*x5*x10*D[f,x5]-(1/y)*x6*x12*D[f,x1]-(1/y)*x6*x13*D[f,x4]-(1/y)*x6*x15*D[f,x2]+(1/y)*x7*x10*D[f,x7]+(1/y)*x8*x9*D[f,x1]-(1/y)*x8*x13*D[f,x5]-(1/y)*x8*x15*D[f,x3]+(1/y)*x9*x11*D[f,x2]+(1/y)*x9*x14*D[f,x4]+3*(1/y)*x10*f+(1/y)*x11*x12*D[f,x3]-(1/y)*x11*x13*D[f,x7]+(1/y)*x12*x14*D[f,x5]+(1/y)*x14*x15*D[f,x7]+x0*D[D[I*f,x7],x1]-x0*D[D[I*f,x5],x2]+x0*D[D[I*f,x4],x3]+x10*D[f,y];
bm11=(1/y)*x1*x11*D[f,x1]-(1/y)*x2*x8*D[f,x1]+(1/y)*x2*x10*D[f,x9]-(1/y)*x2*x14*D[f,x4]+(1/y)*x3*x6*D[f,x1]+(1/y)*x3*x10*D[f,x12]-(1/y)*x3*x14*D[f,x5]+(1/y)*x4*x11*D[f,x4]+(1/y)*x5*x11*D[f,x5]-(1/y)*x6*x7*D[f,x4]-(1/y)*x6*x15*D[f,x9]-(1/y)*x7*x8*D[f,x5]-(1/y)*x7*x10*D[f,x13]-(1/y)*x8*x15*D[f,x12]+(1/y)*x9*x11*D[f,x9]+3*(1/y)*x11*f+(1/y)*x11*x12*D[f,x12]+(1/y)*x11*x13*D[f,x13]-(1/y)*x14*x15*D[f,x13]-x0*D[D[I*f,x13],x1]+x0*D[D[I*f,x12],x4]-x0*D[D[I*f,x9],x5]+x11*D[f,y];
bm12=-(1/y)*x1*x10*D[f,x6]-(1/y)*x1*x13*D[f,x4]-(1/y)*x1*x15*D[f,x2]+(1/y)*x2*x12*D[f,x2]-(1/y)*x3*x9*D[f,x2]+(1/y)*x3*x10*D[f,x11]+(1/y)*x3*x13*D[f,x7]+(1/y)*x4*x12*D[f,x4]-(1/y)*x5*x9*D[f,x4]+(1/y)*x5*x10*D[f,x14]-(1/y)*x5*x15*D[f,x7]+(1/y)*x6*x12*D[f,x6]+(1/y)*x7*x12*D[f,x7]-(1/y)*x8*x9*D[f,x6]-(1/y)*x8*x13*D[f,x14]-(1/y)*x8*x15*D[f,x11]+(1/y)*x11*x12*D[f,x11]+3*(1/y)*x12*f+(1/y)*x12*x14*D[f,x14]-x0*D[D[I*f,x14],x2]+x0*D[D[I*f,x11],x4]-x0*D[D[I*f,x7],x6]+x12*D[f,y];
bm13=(1/y)*x1*x13*D[f,x1]+(1/y)*x2*x13*D[f,x2]+(1/y)*x3*x13*D[f,x3]-(1/y)*x4*x10*D[f,x6]-(1/y)*x4*x12*D[f,x1]-(1/y)*x4*x15*D[f,x2]+(1/y)*x5*x9*D[f,x1]-(1/y)*x5*x10*D[f,x8]-(1/y)*x5*x15*D[f,x3]+(1/y)*x6*x13*D[f,x6]+(1/y)*x7*x9*D[f,x2]-(1/y)*x7*x10*D[f,x11]+(1/y)*x7*x12*D[f,x3]+(1/y)*x8*x13*D[f,x8]-(1/y)*x9*x14*D[f,x6]+(1/y)*x11*x13*D[f,x11]-(1/y)*x12*x14*D[f,x8]+3*(1/y)*x13*f-(1/y)*x14*x15*D[f,x11]-x0*D[D[I*f,x11],x1]+x0*D[D[I*f,x8],x2]-x0*D[D[I*f,x6],x3]+x13*D[f,y];
bm14=(1/y)*x1*x14*D[f,x1]+(1/y)*x2*x14*D[f,x2]+(1/y)*x3*x14*D[f,x3]-(1/y)*x4*x8*D[f,x1]+(1/y)*x4*x10*D[f,x9]-(1/y)*x4*x11*D[f,x2]+(1/y)*x5*x6*D[f,x1]+(1/y)*x5*x10*D[f,x12]-(1/y)*x5*x11*D[f,x3]+(1/y)*x6*x7*D[f,x2]-(1/y)*x6*x13*D[f,x9]+(1/y)*x7*x8*D[f,x3]+(1/y)*x7*x10*D[f,x15]-(1/y)*x8*x13*D[f,x12]+(1/y)*x9*x14*D[f,x9]-(1/y)*x11*x13*D[f,x15]+(1/y)*x12*x14*D[f,x12]+3*(1/y)*x14*f+(1/y)*x14*x15*D[f,x15]+x0*D[D[I*f,x15],x1]-x0*D[D[I*f,x12],x2]+x0*D[D[I*f,x9],x3]+x14*D[f,y];
bm15=(1/y)*x1*x15*D[f,x1]-(1/y)*x2*x10*D[f,x6]-(1/y)*x2*x12*D[f,x1]-(1/y)*x2*x13*D[f,x4]+(1/y)*x3*x9*D[f,x1]-(1/y)*x3*x10*D[f,x8]-(1/y)*x3*x13*D[f,x5]+(1/y)*x4*x15*D[f,x4]+(1/y)*x5*x15*D[f,x5]+(1/y)*x6*x15*D[f,x6]-(1/y)*x7*x9*D[f,x4]+(1/y)*x7*x10*D[f,x14]-(1/y)*x7*x12*D[f,x5]+(1/y)*x8*x15*D[f,x8]-(1/y)*x9*x11*D[f,x6]-(1/y)*x11*x12*D[f,x8]-(1/y)*x11*x13*D[f,x14]+(1/y)*x14*x15*D[f,x14]+3*(1/y)*x15*f+x0*D[D[I*f,x14],x1]-x0*D[D[I*f,x8],x4]+x0*D[D[I*f,x6],x5]+x15*D[f,y];
gm1=y*D[D[I*f,x1],y]+x0*D[D[I*f,x1],x0]+x1*D[D[I*f,x1],x1]+6*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,x5],x1]+x6*D[D[I*f,x6],x1]+x7*D[D[I*f,x5],x2]-x7*D[D[I*f,x4],x3]+x7*x10*(1/y)*D[f,x0]+x8*D[D[I*f,x8],x1]+x9*D[D[I*f,x9],x1]+x10*D[D[I*f,x12],x6]-x10*D[D[I*f,x9],x8]+x11*D[D[I*f,x8],x2]-x11*D[D[I*f,x6],x3]-x11*x13*(1/y)*D[f,x0]+x12*D[D[I*f,x12],x1]+x13*D[D[I*f,x12],x4]-x13*D[D[I*f,x9],x5]+x14*D[D[I*f,x8],x4]-x14*D[D[I*f,x6],x5]+x14*x15*(1/y)*D[f,x0]+x15*D[D[I*f,x12],x2]-x15*D[D[I*f,x9],x3];
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]+6*D[I*f,x2]+x3*D[D[I*f,x3],x2]+x4*D[D[I*f,x4],x2]+x5*D[D[I*f,x7],x1]+x5*D[D[I*f,x4],x3]-x5*x10*(1/y)*D[f,x0]+x6*D[D[I*f,x6],x2]+x7*D[D[I*f,x7],x2]+x8*D[D[I*f,x11],x1]+x8*D[D[I*f,x6],x3]+x8*x13*(1/y)*D[f,x0]+x9*D[D[I*f,x9],x2]+x10*D[D[I*f,x15],x6]-x10*D[D[I*f,x11],x9]+x11*D[D[I*f,x11],x2]+x12*D[D[I*f,x15],x1]+x12*D[D[I*f,x9],x3]-x12*x14*(1/y)*D[f,x0]+x13*D[D[I*f,x15],x4]-x13*D[D[I*f,x9],x7]+x14*D[D[I*f,x11],x4]-x14*D[D[I*f,x7],x6]+x15*D[D[I*f,x15],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,x3],x2]+x3*D[D[I*f,x3],x3]+6*D[I*f,x3]-x4*D[D[I*f,x7],x1]+x4*D[D[I*f,x5],x2]+x4*x10*(1/y)*D[f,x0]+x5*D[D[I*f,x5],x3]-x6*D[D[I*f,x11],x1]+x6*D[D[I*f,x8],x2]-x6*x13*(1/y)*D[f,x0]+x7*D[D[I*f,x7],x3]+x8*D[D[I*f,x8],x3]-x9*D[D[I*f,x15],x1]+x9*D[D[I*f,x12],x2]+x9*x14*(1/y)*D[f,x0]+x10*D[D[I*f,x15],x8]-x10*D[D[I*f,x12],x11]+x11*D[D[I*f,x11],x3]+x12*D[D[I*f,x12],x3]+x13*D[D[I*f,x15],x5]-x13*D[D[I*f,x12],x7]+x14*D[D[I*f,x11],x5]-x14*D[D[I*f,x8],x7]+x15*D[D[I*f,x15],x3];
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*D[D[I*f,x5],x2]+x3*x10*(1/y)*D[f,x0]+x4*D[D[I*f,x4],x4]+6*D[I*f,x4]+x5*D[D[I*f,x5],x4]+x6*D[D[I*f,x6],x4]+x7*D[D[I*f,x7],x4]+x8*D[D[I*f,x14],x1]+x8*D[D[I*f,x6],x5]-x8*x15*(1/y)*D[f,x0]+x9*D[D[I*f,x9],x4]+x10*D[D[I*f,x13],x6]-x10*D[D[I*f,x14],x9]+x11*D[D[I*f,x14],x2]+x11*D[D[I*f,x7],x6]+x11*x12*(1/y)*D[f,x0]+x12*D[D[I*f,x13],x1]+x12*D[D[I*f,x9],x5]+x13*D[D[I*f,x13],x4]+x14*D[D[I*f,x14],x4]+x15*D[D[I*f,x13],x2]+x15*D[D[I*f,x9],x7];
gm5=y*D[D[I*f,x5],y]+x0*D[D[I*f,x5],x0]+x1*D[D[I*f,x5],x1]+x2*D[D[I*f,x7],x1]+x2*D[D[I*f,x4],x3]-x2*x10*(1/y)*D[f,x0]+x3*D[D[I*f,x5],x3]+x4*D[D[I*f,x5],x4]+x5*D[D[I*f,x5],x5]+6*D[I*f,x5]-x6*D[D[I*f,x14],x1]+x6*D[D[I*f,x8],x4]+x6*x15*(1/y)*D[f,x0]+x7*D[D[I*f,x7],x5]+x8*D[D[I*f,x8],x5]-x9*D[D[I*f,x13],x1]+x9*D[D[I*f,x12],x4]-x9*x11*(1/y)*D[f,x0]+x10*D[D[I*f,x13],x8]-x10*D[D[I*f,x14],x12]+x11*D[D[I*f,x14],x3]+x11*D[D[I*f,x8],x7]+x12*D[D[I*f,x12],x5]+x13*D[D[I*f,x13],x5]+x14*D[D[I*f,x14],x5]+x15*D[D[I*f,x13],x3]+x15*D[D[I*f,x12],x7];
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,x6],x2]-x3*D[D[I*f,x11],x1]+x3*D[D[I*f,x8],x2]-x3*x13*(1/y)*D[f,x0]+x4*D[D[I*f,x6],x4]-x5*D[D[I*f,x14],x1]+x5*D[D[I*f,x8],x4]+x5*x15*(1/y)*D[f,x0]+x6*D[D[I*f,x6],x6]+6*D[I*f,x6]-x7*D[D[I*f,x14],x2]+x7*D[D[I*f,x11],x4]-x7*x12*(1/y)*D[f,x0]+x8*D[D[I*f,x8],x6]+x9*D[D[I*f,x9],x6]+x10*D[D[I*f,x10],x6]+x11*D[D[I*f,x11],x6]+x12*D[D[I*f,x10],x1]+x12*D[D[I*f,x9],x8]+x13*D[D[I*f,x10],x4]+x13*D[D[I*f,x14],x9]+x14*D[D[I*f,x14],x6]+x15*D[D[I*f,x10],x2]+x15*D[D[I*f,x11],x9];
gm7=y*D[D[I*f,x7],y]+x0*D[D[I*f,x7],x0]+x1*D[D[I*f,x5],x2]-x1*D[D[I*f,x4],x3]+x1*x10*(1/y)*D[f,x0]+x2*D[D[I*f,x7],x2]+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,x14],x2]+x6*D[D[I*f,x11],x4]-x6*x12*(1/y)*D[f,x0]+x7*D[D[I*f,x7],x7]+6*D[I*f,x7]-x8*D[D[I*f,x14],x3]+x8*D[D[I*f,x11],x5]+x8*x9*(1/y)*D[f,x0]-x9*D[D[I*f,x13],x2]+x9*D[D[I*f,x15],x4]+x10*D[D[I*f,x13],x11]-x10*D[D[I*f,x15],x14]+x11*D[D[I*f,x11],x7]-x12*D[D[I*f,x13],x3]+x12*D[D[I*f,x15],x5]+x13*D[D[I*f,x13],x7]+x14*D[D[I*f,x14],x7]+x15*D[D[I*f,x15],x7];
gm8=y*D[D[I*f,x8],y]+x0*D[D[I*f,x8],x0]+x1*D[D[I*f,x8],x1]+x2*D[D[I*f,x11],x1]+x2*D[D[I*f,x6],x3]+x2*x13*(1/y)*D[f,x0]+x3*D[D[I*f,x8],x3]+x4*D[D[I*f,x14],x1]+x4*D[D[I*f,x6],x5]-x4*x15*(1/y)*D[f,x0]+x5*D[D[I*f,x8],x5]+x6*D[D[I*f,x8],x6]-x7*D[D[I*f,x14],x3]+x7*D[D[I*f,x11],x5]+x7*x9*(1/y)*D[f,x0]+x8*D[D[I*f,x8],x8]+6*D[I*f,x8]-x9*D[D[I*f,x10],x1]+x9*D[D[I*f,x12],x6]+x10*D[D[I*f,x10],x8]+x11*D[D[I*f,x11],x8]+x12*D[D[I*f,x12],x8]+x13*D[D[I*f,x10],x5]+x13*D[D[I*f,x14],x12]+x14*D[D[I*f,x14],x8]+x15*D[D[I*f,x10],x3]+x15*D[D[I*f,x12],x11];
gm9=y*D[D[I*f,x9],y]+x0*D[D[I*f,x9],x0]+x1*D[D[I*f,x9],x1]+x2*D[D[I*f,x9],x2]-x3*D[D[I*f,x15],x1]+x3*D[D[I*f,x12],x2]+x3*x14*(1/y)*D[f,x0]+x4*D[D[I*f,x9],x4]-x5*D[D[I*f,x13],x1]+x5*D[D[I*f,x12],x4]-x5*x11*(1/y)*D[f,x0]+x6*D[D[I*f,x9],x6]-x7*D[D[I*f,x13],x2]+x7*D[D[I*f,x15],x4]+x7*x8*(1/y)*D[f,x0]-x8*D[D[I*f,x10],x1]+x8*D[D[I*f,x12],x6]+x9*D[D[I*f,x9],x9]+6*D[I*f,x9]+x10*D[D[I*f,x10],x9]-x11*D[D[I*f,x10],x2]+x11*D[D[I*f,x15],x6]+x12*D[D[I*f,x12],x9]+x13*D[D[I*f,x13],x9]-x14*D[D[I*f,x10],x4]+x14*D[D[I*f,x13],x6]+x15*D[D[I*f,x15],x9];
gm10=y*D[D[I*f,x10],y]+x0*D[D[I*f,x10],x0]+x1*D[D[I*f,x12],x6]+x1*x7*(1/y)*D[f,x0]-x1*D[D[I*f,x9],x8]-x2*x5*(1/y)*D[f,x0]+x2*D[D[I*f,x15],x6]-x2*D[D[I*f,x11],x9]+x3*x4*(1/y)*D[f,x0]+x3*D[D[I*f,x15],x8]-x3*D[D[I*f,x12],x11]+x4*D[D[I*f,x13],x6]-x4*D[D[I*f,x14],x9]+x5*D[D[I*f,x13],x8]-x5*D[D[I*f,x14],x12]+x6*D[D[I*f,x10],x6]+x7*D[D[I*f,x13],x11]-x7*D[D[I*f,x15],x14]+x8*D[D[I*f,x10],x8]+x9*D[D[I*f,x10],x9]+x10*D[D[I*f,x10],x10]+6*D[I*f,x10]+x11*D[D[I*f,x11],x10]+x12*D[D[I*f,x12],x10]+x13*D[D[I*f,x13],x10]+x14*D[D[I*f,x14],x10]+x15*D[D[I*f,x15],x10];
gm11=y*D[D[I*f,x11],y]+x0*D[D[I*f,x11],x0]+x1*D[D[I*f,x8],x2]-x1*D[D[I*f,x6],x3]-x1*x13*(1/y)*D[f,x0]+x2*D[D[I*f,x11],x2]+x3*D[D[I*f,x11],x3]+x4*D[D[I*f,x14],x2]+x4*D[D[I*f,x7],x6]+x4*x12*(1/y)*D[f,x0]+x5*D[D[I*f,x14],x3]+x5*D[D[I*f,x8],x7]-x5*x9*(1/y)*D[f,x0]+x6*D[D[I*f,x11],x6]+x7*D[D[I*f,x11],x7]+x8*D[D[I*f,x11],x8]-x9*D[D[I*f,x10],x2]+x9*D[D[I*f,x15],x6]+x10*D[D[I*f,x11],x10]+x11*D[D[I*f,x11],x11]+6*D[I*f,x11]-x12*D[D[I*f,x10],x3]+x12*D[D[I*f,x15],x8]+x13*D[D[I*f,x10],x7]+x13*D[D[I*f,x15],x14]+x14*D[D[I*f,x14],x11]+x15*D[D[I*f,x15],x11];
gm12=y*D[D[I*f,x12],y]+x0*D[D[I*f,x12],x0]+x1*D[D[I*f,x12],x1]+x2*D[D[I*f,x15],x1]+x2*D[D[I*f,x9],x3]-x2*x14*(1/y)*D[f,x0]+x3*D[D[I*f,x12],x3]+x4*D[D[I*f,x13],x1]+x4*D[D[I*f,x9],x5]+x4*x11*(1/y)*D[f,x0]+x5*D[D[I*f,x12],x5]+x6*D[D[I*f,x10],x1]-x6*x7*(1/y)*D[f,x0]+x6*D[D[I*f,x9],x8]-x7*D[D[I*f,x13],x3]+x7*D[D[I*f,x15],x5]+x8*D[D[I*f,x12],x8]+x9*D[D[I*f,x12],x9]+x10*D[D[I*f,x12],x10]-x11*D[D[I*f,x10],x3]+x11*D[D[I*f,x15],x8]+x12*D[D[I*f,x12],x12]+6*D[I*f,x12]+x13*D[D[I*f,x13],x12]-x14*D[D[I*f,x10],x5]+x14*D[D[I*f,x13],x8]+x15*D[D[I*f,x15],x12];
gm13=y*D[D[I*f,x13],y]+x0*D[D[I*f,x13],x0]+x1*D[D[I*f,x12],x4]-x1*D[D[I*f,x9],x5]-x1*x11*(1/y)*D[f,x0]+x2*D[D[I*f,x15],x4]-x2*D[D[I*f,x9],x7]+x2*x8*(1/y)*D[f,x0]+x3*D[D[I*f,x15],x5]-x3*x6*(1/y)*D[f,x0]-x3*D[D[I*f,x12],x7]+x4*D[D[I*f,x13],x4]+x5*D[D[I*f,x13],x5]+x6*D[D[I*f,x10],x4]+x6*D[D[I*f,x14],x9]+x7*D[D[I*f,x13],x7]+x8*D[D[I*f,x10],x5]+x8*D[D[I*f,x14],x12]+x9*D[D[I*f,x13],x9]+x10*D[D[I*f,x13],x10]+x11*D[D[I*f,x10],x7]+x11*D[D[I*f,x15],x14]+x12*D[D[I*f,x13],x12]+x13*D[D[I*f,x13],x13]+6*D[I*f,x13]+x14*D[D[I*f,x14],x13]+x15*D[D[I*f,x15],x13];
gm14=y*D[D[I*f,x14],y]+x0*D[D[I*f,x14],x0]+x1*D[D[I*f,x8],x4]-x1*D[D[I*f,x6],x5]+x1*x15*(1/y)*D[f,x0]+x2*D[D[I*f,x11],x4]-x2*D[D[I*f,x7],x6]-x2*x12*(1/y)*D[f,x0]+x3*D[D[I*f,x11],x5]-x3*D[D[I*f,x8],x7]+x3*x9*(1/y)*D[f,x0]+x4*D[D[I*f,x14],x4]+x5*D[D[I*f,x14],x5]+x6*D[D[I*f,x14],x6]+x7*D[D[I*f,x14],x7]+x8*D[D[I*f,x14],x8]-x9*D[D[I*f,x10],x4]+x9*D[D[I*f,x13],x6]+x10*D[D[I*f,x14],x10]+x11*D[D[I*f,x14],x11]-x12*D[D[I*f,x10],x5]+x12*D[D[I*f,x13],x8]+x13*D[D[I*f,x14],x13]+x14*D[D[I*f,x14],x14]+6*D[I*f,x14]-x15*D[D[I*f,x10],x7]+x15*D[D[I*f,x13],x11];
gm15=y*D[D[I*f,x15],y]+x0*D[D[I*f,x15],x0]+x1*D[D[I*f,x12],x2]-x1*D[D[I*f,x9],x3]+x1*x14*(1/y)*D[f,x0]+x2*D[D[I*f,x15],x2]+x3*D[D[I*f,x15],x3]+x4*D[D[I*f,x13],x2]+x4*D[D[I*f,x9],x7]-x4*x8*(1/y)*D[f,x0]+x5*D[D[I*f,x13],x3]+x5*x6*(1/y)*D[f,x0]+x5*D[D[I*f,x12],x7]+x6*D[D[I*f,x10],x2]+x6*D[D[I*f,x11],x9]+x7*D[D[I*f,x15],x7]+x8*D[D[I*f,x10],x3]+x8*D[D[I*f,x12],x11]+x9*D[D[I*f,x15],x9]+x10*D[D[I*f,x15],x10]+x11*D[D[I*f,x15],x11]+x12*D[D[I*f,x15],x12]+x13*D[D[I*f,x15],x13]-x14*D[D[I*f,x10],x7]+x14*D[D[I*f,x13],x11]+x15*D[D[I*f,x15],x15]+6*D[I*f,x15];
mom=-y*D[D[I*f,y],y]-15*(1/y)*I*f+(1/y)*(1/y)*x1*x7*x10*D[f,x0]-(1/y)*(1/y)*x1*x11*x13*D[f,x0]+(1/y)*(1/y)*x1*x14*x15*D[f,x0]-(1/y)*(1/y)*x2*x5*x10*D[f,x0]+(1/y)*(1/y)*x2*x8*x13*D[f,x0]-(1/y)*(1/y)*x2*x12*x14*D[f,x0]+(1/y)*(1/y)*x3*x4*x10*D[f,x0]-(1/y)*(1/y)*x3*x6*x13*D[f,x0]+(1/y)*(1/y)*x3*x9*x14*D[f,x0]-(1/y)*(1/y)*x4*x8*x15*D[f,x0]+(1/y)*(1/y)*x4*x11*x12*D[f,x0]+(1/y)*(1/y)*x5*x6*x15*D[f,x0]-(1/y)*(1/y)*x5*x9*x11*D[f,x0]-(1/y)*(1/y)*x6*x7*x12*D[f,x0]+(1/y)*(1/y)*x7*x8*x9*D[f,x0]-3*(1/y)*x1*D[I*f,x1]-(1/y)*x1*x7*D[D[I*f,x7],x1]+(1/y)*x1*x7*D[D[I*f,x5],x2]-(1/y)*x1*x7*D[D[I*f,x4],x3]-(1/y)*x1*x10*D[D[I*f,x10],x1]+(1/y)*x1*x10*D[D[I*f,x12],x6]-(1/y)*x1*x10*D[D[I*f,x9],x8]-(1/y)*x1*x11*D[D[I*f,x11],x1]+(1/y)*x1*x11*D[D[I*f,x8],x2]-(1/y)*x1*x11*D[D[I*f,x6],x3]-(1/y)*x1*x13*D[D[I*f,x13],x1]+(1/y)*x1*x13*D[D[I*f,x12],x4]-(1/y)*x1*x13*D[D[I*f,x9],x5]-(1/y)*x1*x14*D[D[I*f,x14],x1]+(1/y)*x1*x14*D[D[I*f,x8],x4]-(1/y)*x1*x14*D[D[I*f,x6],x5]-(1/y)*x1*x15*D[D[I*f,x15],x1]+(1/y)*x1*x15*D[D[I*f,x12],x2]-(1/y)*x1*x15*D[D[I*f,x9],x3]-3*(1/y)*x2*D[I*f,x2]+(1/y)*x2*x5*D[D[I*f,x7],x1]-(1/y)*x2*x5*D[D[I*f,x5],x2]+(1/y)*x2*x5*D[D[I*f,x4],x3]+(1/y)*x2*x8*D[D[I*f,x11],x1]-(1/y)*x2*x8*D[D[I*f,x8],x2]+(1/y)*x2*x8*D[D[I*f,x6],x3]-(1/y)*x2*x10*D[D[I*f,x10],x2]+(1/y)*x2*x10*D[D[I*f,x15],x6]-(1/y)*x2*x10*D[D[I*f,x11],x9]+(1/y)*x2*x12*D[D[I*f,x15],x1]-(1/y)*x2*x12*D[D[I*f,x12],x2]+(1/y)*x2*x12*D[D[I*f,x9],x3]-(1/y)*x2*x13*D[D[I*f,x13],x2]+(1/y)*x2*x13*D[D[I*f,x15],x4]-(1/y)*x2*x13*D[D[I*f,x9],x7]-(1/y)*x2*x14*D[D[I*f,x14],x2]+(1/y)*x2*x14*D[D[I*f,x11],x4]-(1/y)*x2*x14*D[D[I*f,x7],x6]-3*(1/y)*x3*D[I*f,x3]-(1/y)*x3*x4*D[D[I*f,x7],x1]+(1/y)*x3*x4*D[D[I*f,x5],x2]-(1/y)*x3*x4*D[D[I*f,x4],x3]-(1/y)*x3*x6*D[D[I*f,x11],x1]+(1/y)*x3*x6*D[D[I*f,x8],x2]-(1/y)*x3*x6*D[D[I*f,x6],x3]-(1/y)*x3*x9*D[D[I*f,x15],x1]+(1/y)*x3*x9*D[D[I*f,x12],x2]-(1/y)*x3*x9*D[D[I*f,x9],x3]-(1/y)*x3*x10*D[D[I*f,x10],x3]+(1/y)*x3*x10*D[D[I*f,x15],x8]-(1/y)*x3*x10*D[D[I*f,x12],x11]-(1/y)*x3*x13*D[D[I*f,x13],x3]+(1/y)*x3*x13*D[D[I*f,x15],x5]-(1/y)*x3*x13*D[D[I*f,x12],x7]-(1/y)*x3*x14*D[D[I*f,x14],x3]+(1/y)*x3*x14*D[D[I*f,x11],x5]-(1/y)*x3*x14*D[D[I*f,x8],x7]-3*(1/y)*x4*D[I*f,x4]+(1/y)*x4*x8*D[D[I*f,x14],x1]-(1/y)*x4*x8*D[D[I*f,x8],x4]+(1/y)*x4*x8*D[D[I*f,x6],x5]-(1/y)*x4*x10*D[D[I*f,x10],x4]+(1/y)*x4*x10*D[D[I*f,x13],x6]-(1/y)*x4*x10*D[D[I*f,x14],x9]+(1/y)*x4*x11*D[D[I*f,x14],x2]-(1/y)*x4*x11*D[D[I*f,x11],x4]+(1/y)*x4*x11*D[D[I*f,x7],x6]+(1/y)*x4*x12*D[D[I*f,x13],x1]-(1/y)*x4*x12*D[D[I*f,x12],x4]+(1/y)*x4*x12*D[D[I*f,x9],x5]+(1/y)*x4*x15*D[D[I*f,x13],x2]-(1/y)*x4*x15*D[D[I*f,x15],x4]+(1/y)*x4*x15*D[D[I*f,x9],x7]-3*(1/y)*x5*D[I*f,x5]-(1/y)*x5*x6*D[D[I*f,x14],x1]+(1/y)*x5*x6*D[D[I*f,x8],x4]-(1/y)*x5*x6*D[D[I*f,x6],x5]-(1/y)*x5*x9*D[D[I*f,x13],x1]+(1/y)*x5*x9*D[D[I*f,x12],x4]-(1/y)*x5*x9*D[D[I*f,x9],x5]-(1/y)*x5*x10*D[D[I*f,x10],x5]+(1/y)*x5*x10*D[D[I*f,x13],x8]-(1/y)*x5*x10*D[D[I*f,x14],x12]+(1/y)*x5*x11*D[D[I*f,x14],x3]-(1/y)*x5*x11*D[D[I*f,x11],x5]+(1/y)*x5*x11*D[D[I*f,x8],x7]+(1/y)*x5*x15*D[D[I*f,x13],x3]-(1/y)*x5*x15*D[D[I*f,x15],x5]+(1/y)*x5*x15*D[D[I*f,x12],x7]-3*(1/y)*x6*D[I*f,x6]-(1/y)*x6*x7*D[D[I*f,x14],x2]+(1/y)*x6*x7*D[D[I*f,x11],x4]-(1/y)*x6*x7*D[D[I*f,x7],x6]+(1/y)*x6*x12*D[D[I*f,x10],x1]-(1/y)*x6*x12*D[D[I*f,x12],x6]+(1/y)*x6*x12*D[D[I*f,x9],x8]+(1/y)*x6*x13*D[D[I*f,x10],x4]-(1/y)*x6*x13*D[D[I*f,x13],x6]+(1/y)*x6*x13*D[D[I*f,x14],x9]+(1/y)*x6*x15*D[D[I*f,x10],x2]-(1/y)*x6*x15*D[D[I*f,x15],x6]+(1/y)*x6*x15*D[D[I*f,x11],x9]-3*(1/y)*x7*D[I*f,x7]-(1/y)*x7*x8*D[D[I*f,x14],x3]+(1/y)*x7*x8*D[D[I*f,x11],x5]-(1/y)*x7*x8*D[D[I*f,x8],x7]-(1/y)*x7*x9*D[D[I*f,x13],x2]+(1/y)*x7*x9*D[D[I*f,x15],x4]-(1/y)*x7*x9*D[D[I*f,x9],x7]-(1/y)*x7*x10*D[D[I*f,x10],x7]+(1/y)*x7*x10*D[D[I*f,x13],x11]-(1/y)*x7*x10*D[D[I*f,x15],x14]-(1/y)*x7*x12*D[D[I*f,x13],x3]+(1/y)*x7*x12*D[D[I*f,x15],x5]-(1/y)*x7*x12*D[D[I*f,x12],x7]-3*(1/y)*x8*D[I*f,x8]-(1/y)*x8*x9*D[D[I*f,x10],x1]+(1/y)*x8*x9*D[D[I*f,x12],x6]-(1/y)*x8*x9*D[D[I*f,x9],x8]+(1/y)*x8*x13*D[D[I*f,x10],x5]-(1/y)*x8*x13*D[D[I*f,x13],x8]+(1/y)*x8*x13*D[D[I*f,x14],x12]+(1/y)*x8*x15*D[D[I*f,x10],x3]-(1/y)*x8*x15*D[D[I*f,x15],x8]+(1/y)*x8*x15*D[D[I*f,x12],x11]-3*(1/y)*x9*D[I*f,x9]-(1/y)*x9*x11*D[D[I*f,x10],x2]+(1/y)*x9*x11*D[D[I*f,x15],x6]-(1/y)*x9*x11*D[D[I*f,x11],x9]-(1/y)*x9*x14*D[D[I*f,x10],x4]+(1/y)*x9*x14*D[D[I*f,x13],x6]-(1/y)*x9*x14*D[D[I*f,x14],x9]-3*(1/y)*x10*D[I*f,x10]-3*(1/y)*x11*D[I*f,x11]-(1/y)*x11*x12*D[D[I*f,x10],x3]+(1/y)*x11*x12*D[D[I*f,x15],x8]-(1/y)*x11*x12*D[D[I*f,x12],x11]+(1/y)*x11*x13*D[D[I*f,x10],x7]-(1/y)*x11*x13*D[D[I*f,x13],x11]+(1/y)*x11*x13*D[D[I*f,x15],x14]-3*(1/y)*x12*D[I*f,x12]-(1/y)*x12*x14*D[D[I*f,x10],x5]+(1/y)*x12*x14*D[D[I*f,x13],x8]-(1/y)*x12*x14*D[D[I*f,x14],x12]-3*(1/y)*x13*D[I*f,x13]-3*(1/y)*x14*D[I*f,x14]-(1/y)*x14*x15*D[D[I*f,x10],x7]+(1/y)*x14*x15*D[D[I*f,x13],x11]-(1/y)*x14*x15*D[D[I*f,x15],x14]-3*(1/y)*x15*D[I*f,x15]-9*D[I*f,y]-x0*D[D[I*f,x0],y]+x0*D[d7*D[f,x10],x1]-x0*D[d11*D[f,x13],x1]+x0*D[d14*D[f,x15],x1]-x0*D[d5*D[f,x10],x2]+x0*D[d8*D[f,x13],x2]-x0*D[d12*D[f,x14],x2]+x0*D[d4*D[f,x10],x3]-x0*D[d6*D[f,x13],x3]+x0*D[d9*D[f,x14],x3]-x0*D[d8*D[f,x15],x4]+x0*D[d11*D[f,x12],x4]+x0*D[d6*D[f,x15],x5]-x0*D[d9*D[f,x11],x5]-x0*D[d7*D[f,x12],x6]+x0*D[d8*D[f,x9],x7]-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]-x6*D[D[I*f,x6],y]-x7*D[D[I*f,x7],y]-x8*D[D[I*f,x8],y]-x9*D[D[I*f,x9],y]-x10*D[D[I*f,x10],y]-x11*D[D[I*f,x11],y]-x12*D[D[I*f,x12],y]-x13*D[D[I*f,x13],y]-x14*D[D[I*f,x14],y]-x15*D[D[I*f,x15],y];