Реализовать программу поиска корней квадратного уравнения

#include <iostream>
#include <cmath>
#include <fstream>
using namespace std;
 
int main()
{
float a, b, c, otv, otv1, otv2, D, x;
cout << "ax^2 + bx +c = 0" "\n";
 
 
ifstream file("/home/user37/file1.txt");
if (file.is_open()){
    cout << "0";
}
file >> a >> b >> c;
 
file.close();
 
//cout « " a = ";
//cin » a;
//cout « " b = ";
//cin » b;
//cout « " c = ";
//cin » c;
D = pow(b,2) - (4*a*c);
if (a == 0 and b ==0) {
        cout << "Pri lubom x";
}
if (a == 0) {
        x =(-c)/b;
        cout << x;
}
if (D > 0){
otv1 = ((-b) + sqrt(D))/(2*a);
otv2 = ((-b) - sqrt(D))/(2*a);
cout << "otv1 = " << otv1 << " otv2 = " << otv2;
};
if (D == 0){
otv = (-b)/(2*a);
cout << otv;
};
if (D < 0){
cout << "Net korney";
};
return 0;
}

ifstream file("file1.txt");

if (a == 0) {
      if (b == 0)
            cout << "Pri lubom x";
      else {
           x =(-c)/b;
           cout << x;
      }
}
else {
     D = pow(b,2) - (4*a*c);
     if (D > 0){
           otv1 = ((-b) + sqrt(D))/(2*a);
           otv2 = ((-b) - sqrt(D))/(2*a);
           cout << "otv1 = " << otv1 << " otv2 = " << otv2;
     };
     if (D == 0){
           otv = (-b)/(2*a);
           cout << otv;
     };
     if (D < 0){
           cout << "Net korney";
     };
}