Entendi. Li sobre optim e optimise mas não vi uma forma de isolar a variável. É que na verdade preciso pegar essa expressão:
der.eqimm=expression(exp(b0 + b1/s + b2/x + b3 * (id1/x) * log(ab1) + b4 * (1 - id1/x) +
b5 * (1 - id1/x) * s) * (b4 * (id1/x^2) - (b3 * (id1/x^2) *
log(ab1) + b2/x^2) + b5 * (id1/x^2) * s)/x - exp(b0 + b1/s +
b2/x + b3 * (id1/x) * log(ab1) + b4 * (1 - id1/x) + b5 *
(1 - id1/x) * s)/x^2)
e fazer com que fique assim: x=b4*id1-b3*id1*log(ab1)-b2+b5*id1*s
para isolar o x usou o maple mas queria no R, se puder tentar me ajudar agradeço! Obrigado