import numpy as np
import sympy as sp
import matplotlib.pyplot as plt
from IPython.display import display
import japanize_matplotlib

n,k=sp.var('n k',integer=True)
ss = sp.summation(1/(k*(k+1)),(k,1,n))
display("部分和は",ss)
lim_s = sp.limit(ss, n, sp.oo)
display("無限級数の収束値は",lim_s)