R语言突变级数法的代码

以下是使用R语言实现突变级数法的代码示例：

# 定义突变级数法函数
mutation_series <- function(x, n) {
  terms <- numeric(n)
  for (i in 1:n) {
    terms[i] <- (-1)^(i+1) * x^(2*i-1) / factorial(2*i-1)
  }
  return(sum(terms))
}

# 示例使用
x <- 0.5
n <- 10
result <- mutation_series(x, n)
cat(sprintf("sin(%.2f) ≈ %.6f", x, result))

此代码将计算 $x=0.5$ 的 $\sin(x)$ 近似值，使用了突变级数法的方法，计算前10项。输出结果为：

sin(0.50) ≈ 0.479426

你可以根据需要，修改x和n的值，来进行不同的计算。