import numpy as np
from scipy.integrate import odeint
import matplotlib.pyplot as plt

# 常微分方程式を定義する
def model(y, t):
    dydt = 1 / (1 + t)
    return dydt

# 初期値を設定する
y0 = 0

# 時間範囲を設定する
t = np.linspace(0, 10, 101)

# 常微分方程式を解く
y = odeint(model, y0, t)

# 結果をプロットする
plt.plot(t, y)
plt.xlabel('t')
plt.ylabel('y')
plt.show()