complx.cpp

Go to the documentation of this file.

 
#include <complex>
#include <utility>
#include <iostream>
    
#include <examples.h>
 
typedef std::complex<double> dcomplex;
 
// Return roots of a quadratic equation.
std::pair<dcomplex, dcomplex>
quadratic (dcomplex a, dcomplex b, dcomplex c) {
 
    dcomplex root = std::sqrt(b * b - 4.0 * a * c);
    a = a * 2.0;
    return std::make_pair ((-b + root) / a, (-b - root) / a);
}
 
 
int main () {
    
    dcomplex a(2, 3);
    dcomplex b(4, 5);
    dcomplex c(6, 7);
    
    std::pair<dcomplex, dcomplex> ans = quadratic(a, b, c);
 
    std::cout << "Roots are " << ans.first << " and "
              << ans.second   << std::endl;
 
    return 0;
}