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