y1, y2, x1, x2, a, b, g, 
x1-1/3*x1^3-y1+g*(x1-x2)
x2-1/3*x2^3-y2+g*(x2-x1)
-b*y1+a+x1
-b*y2+a+x2
