Talvez as funções abaixo possam ajudar:

# mínimo no intervalo de -100 a 100

optimize(fx, c(-100, 100))

# máximo no intervalo de -100 a 100

optimize(fx, c(-100, 100), maximum=T)

# derivada

deriv(~-6.685+8.227*x-0.841*x^2,"x")

# integral

?integrate

etc