Función que permite calcular las raíces de una ecuación cuadrática

Ecua_cuadrat<-function(a,b,c){
  if((b^2-4*a*c)>=0)
  {x1<-(-b+sqrt(b^2-(4*a*c)))/(2*a)
  x2<-(-b-sqrt(b^2-(4*a*c)))/(2*a)
  resultados<-list(x1,x2)
  print(resultados)
  print(paste("su solucion es real y el valor de x1 =",x1))
  print(paste("su solucion es real y el valor de x2 =",x2))}
  else{"El resultado es solucion compleja"}
}

Ejemplo con soluciones reales

a<-1
b<-(-5)
c<-6
Ecua_cuadrat(a,b,c)
## [[1]]
## [1] 3
## 
## [[2]]
## [1] 2
## 
## [1] "su solucion es real y el valor de x1 = 3"
## [1] "su solucion es real y el valor de x2 = 2"

Ejemplo con soluciones complejas

a<-9
b<-3
c<-10
Ecua_cuadrat(a,b,c)
## [1] "El resultado es solucion compleja"