The experimental data from the central composite design runs were utilized for mathematical models far the drilling characteristics containing linear, quadratic and interactive effects of the parameters such as volume fraction of TiC in the composites, drill speed, feed rate and drill diameter. The models were developed via stepwise selection where the insignificant effects were removed using t-test. The models were subjected to optimization of maximizing drill life and satisfying the other constraints.