# 数値微分（前進差分）

# 関数
f <- function(x) {
  return(x ^ 5 - 4 * x ^ 4)
}

# 前進差分
f1 <- function(x, h) {
  return((f(x + h) - f(x)) / h)
}

# 前進差分のグラフ
h <- 1e-1
curve(f1(x, h), xlim = c(-1, 3), ylab = "y")