//polynome de troisieme degre
//exemple de fonction avec des parametres par defaut
#include <iostream>
#include <cmath>     //bibliotheque mathematique

using namespace std;

//declaration de la fonction polynome : a_3 x^3 + a_2 x^2 + a_1 x^1 + a_0
//si a3 n'est pas specifiee, on evaluera un polynome de deuxieme degre
//si a2 et a3 ne sont pas specifiees, on evaluera un polynome de premier degre
double poly(double x, double a0, double a1=0., double a2=0., double a3=0.);

int main() {
  double x = 1.003;
  
  cout << "poly(x,a0=7) = " << poly(x,7) << endl;
  cout << "poly(x,a0=7,a1=6) = " << poly(x,7,6) << endl;
  cout << "poly(x,a0=7,a1=6,a2=5) = " << poly(x,7,6,5) << endl;  
  cout << "poly(x,a0=7,a1=6,a2=5,a3=4) = " << poly(x,7,6,5,4) << endl;

  return 0;
}

//corps de la fonction polynome
double poly(double x, double a0, double a1, double a2, double a3) {
  double p = a0 + (a1 + (a2 + a3 * x) * x) *x;
  //  double p = a0 + a1*x + a2*x*x + a3*x*x*x;
  return p;
}