#include <iostream>
#include <cmath>
#include <iomanip>

using namespace std;

int main() {   
  const int STEPMAX = 20;
   
  double h, dfNum, dfCalc;
  double x = 1.;
  double delta;

  cout << "          h           ;    log(dfNum - dfAna)  ;        dfNum" << endl << endl;
  for (int step=0; step<STEPMAX; step++) {
    h = pow(10.,-step);

    dfNum = (pow(x+h,2) - pow(x,2)) / h;   // derivee numerique
    dfCalc = 2.*x;   // derivee analytique
    delta = log10(fabs(dfCalc - dfNum));   // difference

    cout << setprecision(20) << fixed << h << " " << delta << " " << dfNum << endl;
  }
  cout << endl;

  return 0;
}