k1=(6*a^2-4*a*b+4*b^2);
k2=(3*a^2+5*a*b-5*b^2);
k3=(4*a^2-4*a*b+6*b^2);
k4=(5*a^2-5*a*b-3*b^2);

println(a,b,1,k1,0,0,k1^3,0);
println(a,b,2,k2,0,-1,k2^3,k1^3-k4^3);
println(a,b,3,k3,0,-1,k3^3,0);
println(a,b,4,k4,0,-1,k4^3,k2^3+k3^3);
println(0,0,0,0,1,0,k1^3-k2^3-k3^3-k4^3,0);