def MillionQi(folder, filename, span):
### Read data from file
result, MeasState = Read_File(folder, filename)
freq, linear, phase, mag = Get_Data(result, span)
n = len(freq)
real = [linear[i] * np.cos(np.deg2rad(phase[i])) for i in range(0, n)]
imag = [linear[i] * np.sin(np.deg2rad(phase[i])) for i in range(0, n)]
complext = np.asarray([real[i] + 1j*imag[i] for i in range(0, n)])
inverset = 1/complext
para_guess = 1,0,1e5,1e5,np.median(freq),0
a, a_err, phi0, phi0_err, Qi, Qi_err, Qc_scaled, Qc_scaled_err, fr, fr_err, aphi0, aphi0_err, result, Fit_result = Fitter.Fit_Million(freq, inverset, para_guess)
fitcurve = a * (1 + Qi/Qc_scaled * np.exp(1j*phi0) * 1/(1 + 1j*2*Qi*(freq/fr-1)))* np.exp(1j*aphi0)
#plt.plot(real, imag,'.')
plt.plot(inverset.real, inverset.imag,'.')
plt.plot(fitcurve.real, fitcurve.imag,'.r')
#plt.plot(freq, real)
#plt.plot(freq, inverset.real, 'r')
#plt.plot(freq, imag)
#plt.plot(freq, inverset.imag, 'r')
plt.show()
print np.rad2deg(phi0), len(inverset), len(complext)
