quadratic <- function(a, b, c) {
  res <- b ^ 2 - 4 * a * c
  if (res >= 0){
    solution_1 = (-b + sqrt(res))/(2 * a)
    solution_2 = (-b - sqrt(res))/(2 * a)
    return(c(solution_1, solution_2))
    }
  
}
print(quadratic(1,-1,0))
## [1] 1 0