/**************************************************************************
 *
 * complx.cpp - Complex Number example program.
 *
 * $Id: //stdlib/dev/examples/stdlib/tutorial/complx.cpp#11 $
 *
 ***************************************************************************
 *
 * Copyright (c) 1994-2005 Quovadx,  Inc., acting through its  Rogue Wave
 * Software division. Licensed under the Apache License, Version 2.0 (the
 * "License");  you may  not use this file except  in compliance with the
 * License.    You    may   obtain   a   copy   of    the   License    at
 * http://www.apache.org/licenses/LICENSE-2.0.    Unless   required    by
 * applicable law  or agreed to  in writing,  software  distributed under
 * the License is distributed on an "AS IS" BASIS,  WITHOUT WARRANTIES OR
 * CONDITIONS OF  ANY KIND, either  express or implied.  See  the License
 * for the specific language governing permissions  and limitations under
 * the License.
 * 
 **************************************************************************/

#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;
}