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

calcular_raices<-function(a,b,c){
discriminante<-b^2-(4*a*c)
if(discriminante<0) {message("NaN: RAIZ COMPLEJA, NO TIENE SOLUCION REAL")}
else{
x1<- ((-b + sqrt(discriminante))/(2*a))
x2<- ((-b - sqrt(discriminante))/(2*a))
resultados<-list(x1,x2)
names(resultados)<- c("X1","X2")
print(resultados)
}}
# Verificando la función 
ej1<-calcular_raices(2,4,2)
## $X1
## [1] -1
## 
## $X2
## [1] -1
ej2<-calcular_raices(3,5,4)
## NaN: RAIZ COMPLEJA, NO TIENE SOLUCION REAL