f / (df/dy) for f = y^2 + a^2(1-(y-w)^2/b^2) - r^2