xy' = (y^2) * lnx - y