Wyznaczanie miejsc zerowych w równaniu kwadratowym w języku C++

//program wyznaczajacy miejsca zerowe 
//rownania kwadratowego
//w postaci f(x)=ax^2 + bx + c 


#include <iostream>
#include <cmath>

using namespace std;

int main (void){

//definiujemy zmienne
double a, b, c, delta, x1, x2;

a=1;
b=4; //dane 
c=1;

//wyznaczanie delty
delta = pow(b,2)-4*a*c;
//wypisanie w konsoli wartosci delty
cout << "delta = " << delta <<endl;
//instrukcja warunkowa
//jezeli delta jest wieksza od 0 to wykonaj instrokcje
if (delta>0){
 
 x1 = (-b-sqrt(delta))/(2*a);
 x2 = (-b+sqrt(delta))/(2*a);
 
 cout << "x1 = " << x1 << endl; 
 cout << "x2 = " << x2 << endl;
}

//jezeli delta jest rowna 0 to ...
//porownujemy za pomoca operatora ==
//jezeli damy jeden znak = to program 
//nie bedzie dzialal prawidlowo
if (delta==0){

 x1 = (-b)/(2*a);
  
 cout << "x1 = " << x1 << endl; 

}
//jezeli delta nie jest dodatnia ani rowna 0 
//to musi byc ujemna

//else oznacz "w przeciwnym razie"
else{
 cout << "nie mam miejsc zerowych"<<endl;
}

return 0;
}///EOF