def METHOD_NAME(x, y):
    # Temps to make param estimations
    ya = np.abs(y[:])
    phi_0 = cmath.phase(y[0])
    phi_v = 0
    s21min = np.abs(min(ya)/max(ya))
    y_rotated = y*np.exp(-1.j*(phi_v*x+phi_0-phi_v*x[0]))
    y_min = np.argmin(y_rotated.imag)
    y_max = np.argmax(y_rotated.imag)
    # Parameters
    A = (ya[0]+ya[-1])/2.
    f0 = x[ya[:].tolist().index(min(ya[:]))]
    Q = f0/(abs(x[y_min]-x[y_max]))
    Qca = np.abs(Q/np.abs(1-s21min))
    #-440.
    phi_0 = phi_0
    p0i = (A, f0, Q, Qca, phi_0, phi_v)
    return p0i