A new algorithm called homotopy iteration method based on the homotopy function is studied and improved. By the improved homotopy iteration method, polynomial systems with high order and deficient can be solved fast and efficiently comparing to the original homotopy iteration method. Numerical examples for the ninepoint path synthesis of four-bar linkages show the advantages and efficiency of the improved homotopy iteration method.