y' = p(y)
y'' = p * dp/dy
p * dp/dy * y^3 = 1
p dp = 1/y^3 dy
int p dp = int dy/y^3
1/2 * p^2 = -1/2 * 1/y^2 + C1
p^2 = C1 - 1/y^2