import numpy as np
from ooragan import plot_triptych

# Generate data from theoretical model
def S21(f, Q, Qc, f0):
    return 1 - (Q / Qc) / (1 + 2j * Q * (f - f0) / f0)

Q = 85000
Qc = 100000
f0 = 6.5e9
tau = 0.0000001
a = 1
alpha = 0
phi = 0
frequency = np.linspace(6.4995e9, 6.5005e9, 5001)

fig = plot_triptych(
    frequency,
    S21(frequency, Q, Qc, f0),
)
fig.show()