X = var('X')
def R(X):
    return numerical_integral(2/((n+1)*pi), 0, X)[0]
n = abs(X)/pi
i = 1
while R(i) > 10^(-3):
    i += 1

I = F(i)+R(i)
I.numerical_approx()