(* Relation of N=13 for y,z,v,w,v *)

f0[1]=Sum[q^(13*(13/2(p+1/26)^2-1/104))*(-1)^p,{p,-30,30}];
f0[2]=Sum[q^(13*(13/2(p+3/26)^2-1/104))*(-1)^p,{p,-30,30}];
f0[3]=Sum[q^(13*(13/2(p+5/26)^2-1/104))*(-1)^p,{p,-30,30}];
f0[4]=Sum[q^(13*(13/2(p+7/26)^2-1/104))*(-1)^p,{p,-30,30}];
f0[5]=Sum[q^(13*(13/2(p+9/26)^2-1/104))*(-1)^p,{p,-30,30}];
f0[6]=Sum[q^(13*(13/2(p+11/26)^2-1/104))*(-1)^p,{p,-30,30}];
x=f0[1];
y=f0[2];
z=f0[3];
v=f0[4];
w=f0[5];
u=f0[6];
Print["1"]
Series[x^2*z*w-x*y*v^2+z*v*u^2,{q,0,40000}]
Print["2"]
Series[x^2*z*w-y*z^2*v-x*y*w*u+z^2*w^2,{q,0,40000}]
Print["3"]
Series[x*y*z*w-y^2*v^2-x*v^2*u+y*v*w^2,{q,0,40000}]
Print["4"]
Series[z^2*v*u-y*v*w^2+x*w*u^2,{q,0,40000}]
Print["5"]
Series[x*z*v^2-y^2*z*w-y*w*u^2,{q,0,40000}]
Print["6"]
Series[x^3*u-x*z*v^2+x*z*w*u+v*w^2*u,{q,0,40000}]
Print["7"]
Series[x^3*u-z^3*v+v^3*w,{q,0,40000}]
Print["8"]
Series[x^2*y*u-y*z*v^2+z*v*w^2,{q,0,40000}]
Print["9"]
Series[x^3*z-x*y^3+w*u^3,{q,0,40000}]
Print["10"]
Series[x^2*y*u-x*z^2*w+y*z*w*u+w*u^3,{q,0,40000}]
Print["11"]
Series[x*y^2*u-y*z^2*w+z*v^2*u,{q,0,40000}]
Print["12"]
Series[x*y^2*u-x*y*v*w+x^2*u^2+z*w*u^2,{q,0,40000}]
Print["13"]
Series[x^2*z*u-y^2*v*w+z^2*w*u,{q,0,40000}]
Print["14"]
Series[x*v^3-y^3*u-y*w^3,{q,0,40000}]
Print["15"]
Series[x*v^3-y^2*v*w+x*y*u^2-x*v*w*u,{q,0,40000}]
Print["16"]
Series[x^3*v-x^2*z^2+x*y*z*u-z^3*w+v^2*w^2,{q,0,40000}]
Print["17"]
Series[x*y*z*u-x*z*v*w+y*v*w*u,{q,0,40000}]
Print["18"]
Series[x^3*v-y^3*z+w^3*u,{q,0,40000}]
Print["19"]
Series[x*z*v*w-y*v^3-y^2*u^2+v^2*w^2-w^3*u,{q,0,40000}]
Print["20"]
Series[,{q,0,40000}]
Print["21"]
Series[y*v^3-z^3*w-x*u^3,{q,0,40000}]
Print["22"]
Series[x^2*y*v-x*y*z^2+v*w*u^2,{q,0,40000}]
Print["23"]
Series[y^2*z*u-y*z*v*w+z*w^3-v*w*u^2,{q,0,40000}]
Print["24"]
Series[z*w^3-v^3*u-y*u^3,{q,0,40000}]
Print["25"]
Series[x*y^2*v-y^2*z^2-y*z*u^2+z*v*w*u,{q,0,40000}]
Print["26"]
Series[x^2*v*u-x*y*w^2+y*z*u^2,{q,0,40000}]
Print["27"]
Series[x*z^3-y^3*v-z*u^3,{q,0,40000}]
Print["28"]
Series[x^2*z*v-y^3*v-x*y*v*u+y*z^2*u,{q,0,40000}]
Print["29"]
Series[x*v^2*w-y*z^2*u-x*w^2*u,{q,0,40000}]
Print["30"]
Series[x^3*w-y*z^3+v*w^3,{q,0,40000}]
Print["31"]
Series[x*y*z*v-y*z^3-x*z*w^2+y*v^2*w,{q,0,40000}]
Print["32"]
Series[x*z*w^2-y^2*v*u-x*v*u^2,{q,0,40000}]
Print["33"]
Series[x^2*y*w-y^2*z*v+v^2*w*u,{q,0,40000}]
Print["34"]
Series[x*z*v*u-y*z*w^2+v^2*w*u-w^2*u^2,{q,0,40000}]
Print["35"]
Series[x*y^2*w-x*z^2*v+z*w^2*u,{q,0,40000}]