import numpy as np
import matplotlib.pyplot as plt
from math import pi

N = 300
fr = np.logspace(-3, 3, N)
w = 2 * pi * fr
s = w * 1j

R1 = 1e6
R2 = 1e4
C = 1.6e-6 / 2 / pi

num = R1 * R2 * C * s + R2
den = R1 * R2 * C * s + R1 + R2
H = num / den

plt.figure(1)
plt.subplot(211)
plt.loglog(fr, abs(H))
plt.xlabel('Frequency (Hz)')
plt.ylabel('Magnitude (V/V)')
plt.title('R1 = 1MΩ, R2 = 1kΩ,C = 1.6μF, Rf = 100kΩ')

plt.subplot(212)
plt.semilogx(fr, np.angle(H) * 180 / pi)
plt.xlabel('Frequency (Hz)')
plt.ylabel('Phase (deg)')

plt.show()