import matplotlib.pyplot as plt
import numpy as np
f=[............]
G=[............]
phi=[............]
ftest=np.linspace(0,10000,1000)
freq0=............
def g(freq):
    return 1/np.sqrt(1+freq**2/freq0**2)
def deph(freq):
    return -np.arctan(freq/freq0)
Gtest=g(ftest)
phitest=deph(ftest)

plt.subplot(2,1,1)
plt.plot(ftest,Gtest,'b-')
plt.plot(f,G,'r+')
plt.xscale("log")
plt.xlabel("f(Hz)")
plt.ylabel("G")

plt.subplot(2,1,2)
plt.plot(ftest,phitest,'b-')
plt.plot(f,phi,'k+')
plt.xscale("log")
plt.xlabel("f(Hz)")
plt.ylabel("phi(rad)")
plt.show()