Day24

// These statements define a class for complex numbers.
// This declaration is stored in complex.h.
#include <iostream>
#include <cmath>
using namespace std;
class complex
{
// Declare function prototypes for public members.
public:
complex();
complex(double a, double b);
void print();
void input();
double magn(complex);
double angle(complex);
complex operator+(complex);
complex operator-(complex);
complex operator*(complex);
complex operator/(complex);
// Declare private members.
double real, imag;
private:

};
//
// Class implementation
//
// These statements define implementation of a complex class.
// This function is the default constructor to initialize a
// complex number that is not given a value.
complex::complex()
{
real = 0;
imag = 0;
}
// This function is the constructor to initialize a complex
// number to a specified value.
complex::complex(double a, double b)
{
real = a;
imag = b;
}
// This function prints a complex number.
void complex::print()
{
if (imag > 0)
cout << real << "+" << imag << "i" << endl;
else
if (imag == 0)
cout << real << endl;
else
cout << real << imag << "i" << endl;
}
// This function reads two values for a complex number.
void complex::input()
{
cin >> real >> imag;
}
// This function defines the sum of complex numbers.
complex complex::operator+(complex c)
{
// Definition of complex addition.
complex temp;
temp.real = c.real + real;
temp.imag = c.imag + imag;
return temp;
}
// This function defines the difference of complex numbers.
complex complex::operator-(complex c)
{
// Definition of complex subtraction.
complex temp;
temp.real = real - c.real;
temp.imag = imag - c.imag;
return temp;
}
// This function defines the product of complex numbers.
complex complex::operator*(complex c)
{
// Definition of complex multiplication.
complex temp;
temp.real = (real*c.real - imag*c.imag);
temp.imag = (imag*c.real + real*c.imag);
return temp;
}





complex complex::operator/(complex c)
{
// Definition of complex division.
complex temp;
temp.real = (real*c.real + imag*c.imag) /
(pow(c.real, 2) + pow(c.imag, 2));
temp.imag = (imag*c.real - real*c.imag) /
(pow(c.real, 2) + pow(c.imag, 2));
return temp;
}

int main(void)
{
// Declare and initialize variables.
complex c1(4, 1), c2(-3, -2), c3;
// Print initial values.
cout << "c1:";
c1.print();
cout << "c2:";
c2.print();
cout << "c3:";
c3.print();
// Compute and print new values.
c3 = c1 + c2;
cout << "c1+c2 = ";
c3.print();
c3 = c1 - c2;
cout << "c1-c2 = ";
c3.print();
c3 = c1*c2;
cout << "c1*c2 = ";
c3.print();
c3 = c1 / c2;
cout << "c1/c2 = ";
c3.print();



// Declare and initialize variables.
double a, b, c, term1, disc;
complex root1, root2;
// Read values for a, b, c from the keyboard.
cout << "Enter real values a, b, c:" << endl;
cin >> a >> b >> c;
// Compute roots of quadratic equation.
term1 = -b / (2 * a);
disc = b*b - 4 * a*c;
if (disc >= 0)
{
root1.real = term1 + sqrt(disc) / (2 * a);
root2.real = term1 - sqrt(disc) / (2 * a);
}
else
{
root1.real = term1;
root1.imag = sqrt(-disc) / (2 * a);
root2.real = term1;
root2.imag = -sqrt(-disc) / (2 * a);
}
// Print roots.
cout << "Roots:" << endl;
if (root1.imag == root2.imag&&root1.real == root2.real)
cout << "double root" << endl;
else
if (root1.imag == 0 && root2.imag == 0)

cout << "real roots" << endl;
else
cout << "complex roots" << endl;
cout.setf(ios::fixed);
cout.precision(2);
root1.print();
root2.print();





system("pause");
// Exit program.
return 0;
}