# Vedi qui per approfondimenti
# source("http://macosa.dima.unige.it/r.R")
f = function(x) 3/2*x^4-12*x^3+25*x^2-x-45
BF=4; HF=3; Plane(-2,5, -150,150); graph2(f, -2,5, "black")
deriv(f,"x")
# 3/2*(4*x^3) - 12*(3*x^2) + 25*(2*x) - 1
# Se voglio esplicitare il valore dei coefficienti:
x=1; c(3/2*(4*x^3), - 12*(3*x^2), + 25*(2*x), - 1)
#  6 -36  50  -1
df <- function(x) eval(deriv(f,"x"))
deriv2(f,"x")
# 3/2*(4*(3*x^2)) - 12*(3*(2*x)) + 25*2
c( 3/2*(4*(3*x^2)), - 12*(3*(2*x)), + 25*2)
#  18 -72  50
d2f <- function(x) eval(deriv2(f,"x")); graph2(d2f,-3,5, "brown")
deriv3(f,"x")
# 3/2*(4*(3*(2*x))) - 12*(3*2)
c(3/2*(4*(3*(2*x))), - 12*(3*2))
#  36 -72
d3f <- function(x) eval(deriv3(f,"x")); graph2(d3f,-3,5, "magenta")
deriv4(f,"x")
# 3/2*(4*(3*2))
3/2*(4*(3*2))
#  36
d4f <- function(x) eval(deriv4(f,"x")); graph2(d4f,-3,5, "seagreen")
deriv5(f,"x")
#  0
d5f <- function(x) eval(deriv5(f,"x")); graph2(d5f,-3,5, "red")
text(-1.3,80,"f"); text(-1.3,-80,"f'"); text(-0.5,120,"f''")
text(1.5,60,"f'''"); text(-1.5,-15,"f''''")