#1.	Realizar la expansión de la función (a-x) ^6 programar ambas funciones con secuencia (0.995
#hasta 1.005 saltando cada vez 0.001, restar ambas funciones y graficar el error acumulado,
#explicar por qué cree que las líneas del error acumulado resultan de esa manera).

f = function(x) (1-x)^6
ef = function(x) x^6-6*x^5+15*x^4-20* x^3+15* x^2-6* x+1
#a=options(scipen = 999)
x = seq(0.995, 1.005, by=0.001)
y1 = f(x)
y2 = ef(x)
(y1-y2)

x

# grafico


x = seq(0.995, 1.005, by=0.001)
a = x[1]; b = x[length(x)]  #a=x[1]# elige el primer valor de la secuencia en x(0.995) (x = seq(0.995, 1.005, by=0.001))
#b = x[length(x)]: es donde termina la secuencia de x (1.005)
curve((1-x)^6, a, b, col="green"); abline(h=0,v=0, lty=3)
curve(x^6-6*x^5+15*x^4-20* x^3+15* x^
        2-6* x+1, a, b, col="blue",add=T)

x[length(x)]