Квадратное уравнение в R

Пример для версий R 2.10.1

Основной средой работы с R является интерактивный интерпретатор, поэтому достаточно определить функцию, решающую уравнение, и вызывать ее с нужными коэффициентами. Используется встроенный тип complex.

quadratic <- function(a, b, c) {
  if (a == 0) {
    print("Not a quadratic equation.")
  } else {
    d <- b * b - 4 * a * c
    if (d == 0) {
      print(paste("x =", -b / 2 / a))
    } else
      if (d > 0) {
        print(paste("x =", (-b + sqrt(d) * c(1, -1)) / 2 / a))
      } else {
        print(paste("x =", (-b + sqrt(as.complex(d)) * c(1, -1)) / 2 / a))
      }
  }
}

Квадратное уравнение в R

Комментарии