function a = gray_polyfit(x,y)

a = (y - x^2)/(x - x^2);
