# Author: Vinickw, 2023-12-14

import matplotlib.pyplot as plt
import numpy as np

def f(x):
  return 1/x

x = np.linspace(0.01, 8, 1000)
y = f(x)

fig, ax = plt.subplots(figsize=(6, 6))
ax.plot(x, y)

ax.set_xlim(0, 8)
ax.set_ylim(0, 8)

plt.xticks([1,np.e,np.exp(2)],labels=["1","e","e^2"])
plt.yticks([1,np.e,np.exp(2)],labels=["1","e","e^2"])

plt.fill_between(x, 1, where=(x>=0) & (x<=1),color='blue', alpha=0.2)

plt.fill_between(x, y, where=(x>1) & (x<=np.e),color='red', alpha=0.2)
plt.fill_betweenx(x, y, where=(x>1) & (x<=np.e),color='red', alpha=0.2)

plt.fill_between(x, y, where=(x>np.e) & (x<=np.exp(2)),color='blue', alpha=0.2)
plt.fill_betweenx(x, y, where=(x>np.e) & (x<=np.exp(2)),color='blue', alpha=0.2)

plt.plot(x, y, color="blue")
plt.savefig("xy=1.svg", format='svg')