I am currently trying to plot some geodesics with maple, in order to do this I'm trying to get Maple to give me an expression for dr/dy (y here is actually a radial variable) from my equation system which is given by

sys := [diff(y(x), x) = sqrt((1-(diff(r(x), x))^2)/(r(x)^4+r(x)^2)), diff(r(x), x) = sqrt((2+1/r(x)^4-1/r(x)^2)/(1/r(x)^4-1)^2)]

Can anyone help me with this?