(1-x*y)*S(x,y) = 1 + x*y*(1+y)*dS/dy