from gabarit import analyze_filter

T = 1
square = lambda t: 0.5 if t - t//T*T < T/2 else -0.5

analyze_filter(
    input_signal=square,
    period=T,
    sample_period=0.01,
    recording_periods=4.56,
)