如何求解R中的二次方程？

solve_quadratic <- function(a, b, c) {
  discriminant <- b^2 - 4*a*c
  
  if (discriminant > 0) {
    # 两个不同的实根
    root1 <- (-b + sqrt(discriminant)) / (2*a)
    root2 <- (-b - sqrt(discriminant)) / (2*a)
    return(list(root1 = root1, root2 = root2))
  } else if (discriminant == 0) {
    # 一个实根（重根）
    root <- -b / (2*a)
    return(list(root = root))
  } else {
    # 两个共轭复根
    real_part <- -b / (2*a)
    imaginary_part <- sqrt(-discriminant) / (2*a)
    root1 <- complex(real = real_part, imaginary = imaginary_part)
    root2 <- complex(real = real_part, imaginary = -imaginary_part)
    return(list(root1 = root1, root2 = root2))
  }
}

result <- solve_quadratic(a = 1, b = -5, c = 6)
print(result)