def METHOD_NAME(x,y,**kw):
    '''
    find s initial values for gaussian decaying Ramsey
    '''
    y_avg = mean(y)
    fty = abs(fft(y-y_avg))
    A = max(fty)/len(fty)*2*pi
    hl= int(len(fty)/2)
    f_ind =  ((fty.tolist()).index(max(fty[:hl]))-1)
    f =f_ind/max(x)
    phi = arctan(imag(fty[f_ind])/real(fty[f_ind]))
    flh = abs(fty[:hl]) > (0.45*max(abs(fty[:hl])))
    ind1 = (flh.tolist()).index(True)
    k=0
    while flh[ind1+k]:
        k+=1
    sigm = (1.*k)/max(x)/2
    #print sigm
    T2s = 1./(pi*sigm)
    #print T2s
    return y_avg,A,f,phi,T2s