A method is described which can be used to design a wide class of nonrecursive linear-phase filters including those with arbitrary magnitude specifications and with time domain constraints. The approach is based on formulating the weighted mean-square error between the amplitude responses of the practical and ideal filters as a quadratic function. The filter coefficients are obtained by solving a set of linear equations. This method leads to a lower weighted mean-square error and is computationally more efficient than both the eigenfilter method and the method based on the Remez exchange algorithm. However, the main advantage of our method lies in its computational efficiency. Design examples of many different types of filters are provided.