can we just make this `dydx = idiff(eq, y, x)`? and use `solve([slope + 1/dydx, eq], [x, y])` below?