Alexandre dos Santos
Engenheiro Florestal, MSc.
Universidade Federal de Lavras
Departamento de Entomologia
Laboratório de Entomologia Florestal
Caixa Postal 3037
37200-000 - Lavras/MG
Fone: +55 (35) 3829-5122

Tomei a liberdade de simplificar um pouco sua expressao (e'
equivalente ao modelo que vc descreve):

txt = expression(frac(1, 1+e^{0.9588-0.0841*x+0.0010*x^2}))
plot(1, type='n')
legend(1, 1, txt, bty='n')

b