これまでの関数とは異なり、有理関数をPythonでプロットする際には特別な注意が必要です。未定義点や無限大となる値が存在するため、定義域を分割してエラーを防ぐ必要があります。

1. 関数の定義

有理関数は次のように定義します:

def rational_function(x):
    return 1 / (x - 1)

主な注意点:

• $x = 1$ はゼロ除算を避けるため、計算から除外;
• 関数は2つの定義域（$x = 1$ の左側と右側）に分割。

2. 定義域の分割

ゼロ除算を避けるため、2つの異なるx値のセットを生成:

x_left = np.linspace(-4, 0.99, 250)  # Left of x = 1
x_right = np.linspace(1.01, 4, 250)  # Right of x = 1

値 0.99 および 1.01 により、$x = 1$ を含めず、エラーを防止。

3. 関数のプロット

plt.plot(x_left, y_left, color='blue', linewidth=2, label=r"$f(x) = \frac{1}{x - 1}$")
plt.plot(x_right, y_right, color='blue', linewidth=2)

関数は $x = 1$ でジャンプするため、2つの部分に分けてプロットする必要がある。

4. 漸近線と切片の表示

• 垂直漸近線（$x = 1$）:

plt.axvline(1, color='red', linestyle='--',
            linewidth=1, label="Vertical Asymptote (x=1)")

• 水平漸近線 ($y = 0$):

plt.axhline(0, color='green', linestyle='--', 
            linewidth=1, label="Horizontal Asymptote (y=0)")

• Y切片（$x = 0$）:

y_intercept = rational_function(0)
plt.scatter(0, y_intercept, color='purple', label="Y-Intercept")

5. 方向矢印の追加

関数が無限に伸びることを示すために:

plt.annotate('', xy=(x_right[-1], y_right[-1]), xytext=(x_right[-2], y_right[-2]), arrowprops=dict(arrowstyle='->', color='blue', linewidth=1.5))