Solución de la ecuación cuadrática ax2 + bx + c=0

d <- discriminante b2 - 4ac

f1 <- function(a, b, c) {
    if (a == 0 && b == 0 && c == 0) {
        raices <- "Todos los reales"
    } else if (a == 0 && b == 0 && c != 0) {
        raices <- "no hay solución"
    } else if (a == 0 && b != 0 && c != 0) {
        raices <- -c/b
    } else {
        d <- b^2 - 4 * a * c
        if (d > 0) {
            raices <- (-b + c(1, -1) * sqrt(d))/2 * a
        } else if (d == 0) {
            raices <- -b/(2 * a)
        } else {
            raices <- (-b + c(1, -1) * (complex(1, 0, sqrt(-d))/2 * a))
        }
    }
    return(cat("soluciones:", raices))
}


f1(0, 0, 0)
## soluciones: Todos los reales
f1(1, 0, 1)
## soluciones: 0+1i 0-1i
f1(1, 5, 6)
## soluciones: -2 -3
f1(1, 2, 1)
## soluciones: -1