In this paper, an algorithm with linear iteration form for inverse isoparametric mapping in FEM is built with Taylor expansion technique. Its convergence is proved in theory, and its convergence rate is tested in numeric. The algorithm is efficient and practicable for its simplicity in iteration form, convenience in computer implementation and wide adaptation in element type. The flow chart of the algorithm is presented.