s, y1, y2, j, 
4*j*y1*y2^2+20*s*y1*y2-3*y1^2*y2-4*j*y1-3*y2
3*j*y1^2*y2+4*y1*y2^2-3*j*y2+4*y1
j^2+1
