`def` downSampling(num,den,M,Tc,points):
    omega = np.logspace(np.log10(0.0001),np.log10(0.99*np.pi/(Tc)),num=points,endpoint=True)
    denFreq=np.zeros(omega.size,dtype=complex)
    numFreq=np.zeros(omega.size,dtype=complex) 
    bode=np.zeros(omega.size,dtype=complex) 
    M1=float(M)

    
    for k in range(M):
        
        zeta=np.exp(1.0j*omega*Tc/M)*np.exp(-2.0j*np.pi*k*Tc/M)
        denFreq=denFreq+np.polyval(den, zeta)
        numFreq=numFreq+np.polyval(num, zeta)  
        bode=bode+(numFreq/denFreq)
    bode=bode/M
   
    
    sys=ct.FrequencyResponseData(bode,omega,smooth=True)
    return sys