#include <iostream>
#include <cmath>

using namespace std;

//prototype de la fonction
double fonc(double);    

int main() {
  double x;
  cout << "Entrez une valeur a laquelle calculer les derivees :" << endl;
  cout << "x = ";   cin >> x;
  double h;
  cout << "Entrez une valeur de h :" << endl;
  cout << "h = ";   cin >> h;
  
  double derivee1, derivee2;
  derivee1 = (fonc(x+h) - fonc(x-h)) / (2.*h);
  derivee2 = (fonc(x+h) - 2.*fonc(x) + fonc(x-h)) / (h*h);
  
  cout << "Au premier ordre en h la derivee premiere vaut :" << endl;
  cout << "f'(x) = " << derivee1 << endl;
  cout << "La derivee seconde vaut (attention a h !) : " << endl;
  cout << "f''(x) = " << derivee2 << endl;

  return 0;
}

//definition de la fonction
//dans cet example nous avons choisi un polynome de 2eme degree;
double fonc(double x) {
  double a, b, c;
  a = 10.;
  b = -3.;
  c = 7.;
  
  double y = a*(x*x) + b*x + c;

  return y;
}