# 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''''")