# -*- coding: utf-8 -*-
"""
Created on Fri Oct 16 09:43:08 2026

@author: François Coulombeau
"""
import numpy as np
n = 1000000
p = 3
S = 0
for k in range(n,p*n+1):
    S = S+1/k
print(S, "est une approximation de", np.log(p))