double harmonic_number ( int n ); double steinerberger_function ( int n, double x ); double steinerberger_integral01 ( int n );