#include <iostream>
#include <math.h>


double *harmonic_series(int n) {
  double sum = 0;
  double *series = new double[n];

  series[0] = 0;

  for (int i=1; i<n; i++) {
    sum = sum + 1.0/i;
    series[i] = sum;
  }
  
  return series;
}

int main() {
  double *h = harmonic_series(10000);
  for (int i=0; i<10000; i++) 
    cout << h[i] << '\t' 
         << log(i) << '\t' 
         << h[i] - log(i)  << endl;

  delete [] h;
}