/*
    Polynomial math:

    Defining and evaluating a polynomium 
*/



#include <iostream>
using namespace std;

// This is the header defining DualNumber class.
#include "Polynomium.hpp"


int main()
{
   cout << "Defining and printing polynomia:\n";
   Polynomium P;
   P[4]=3.0; P[1]=5.0;                // 3x^4 + 5x^1
   cout << "P=" << P << "\n";

   Polynomium Q;
   Q[4]=2.0;  Q[2]=1.0;  Q[0]=10.0;   // 2x^4 + 1x^2 + 10 
   cout << "Q=" << Q << "\n";

   Polynomium R;
   cout << "R=" << R << "\n\n";


   double x=1.0;
   cout << "\nEvaluating polynomia at x=" << x << "\n";
   cout << "P(x)=" << P(x) << "\n";
   cout << "Q(x)=" << Q(x) << "\n";
   cout << "R(x)=" << R(x) << "\n";

   return 0;
}