This paper presents an approximation method for simplifying a high-order transfer function to a low-order transfer function. A model reduction is based on minimizing the error function weighted by the numerator polynomial of reduced systems. The proposed methods provide better low frequency fit and a computer aided algorithm that estimates the coefficients vector for the numerator and denominator polynomial on the simplified systems from an overdetermined linear system constructed by frequency responses of the original systems. Two examples are given to illustrate the feasibilities of the suggested schemes.