2410.sagews

for x in range(1,5):
    for y in range(1,5):
        z=2*x*(y^4)+x^2*y^3-2*x^3*y^2-y^5-x^4*y+2*y;
        if z > 0:
             print x,y,z;

1 1 1 1 2 2 2 3 3