The method was extended to reduce feed costs (i.e., economic factors) and nitrogen and phosphorus excretions (i.e., environmental factors) by introducing penalty coefficients of nitrogen and phosphorus contents in each ingredient of the diet in the objective function in traditional linear programming.