유리수의 자연 로그를 임의의 정밀도 이상의 정밀도로 유리수로 근사

Python

import numpy as np

n=0.03 #input
m=n if n>=1 else 1/n 
#print(m)
p=0.0001 #input


def error(m,n):
    denom=2*m*(2*n+1)*(m+1)**(2*n-1)    
    return ((m+1)**2)*((m-1)**(2*n-1))/denom

def ap_A(n,m):
    s=0
    for k in range(1,n+1):
       r=((m-1)/(m+1))**(2*k-1)
       s+=r/(2*k-1)
       #print(s)
    return 2*s

def n_finder(m,p):
    n=1
    while error(m,n)>p:
        n+=1
        #print(n,error(m,n))
    return n

ap=ap_A(n_finder(m,p),m)
if n<1: ap=-ap
print("An approximate value A s.t. |ln(n)-A|<p is ", ap) # output of result
print(abs(np.log(n)-ap)<p) # check