The isogeometric collocation method (IGA-C), which is a promising branch of analysis (IGA), can be considered fitting the load function with combination numerical solution and its derivatives. In this study, we develop an iterative method, least-squares progressive-iterative approximation (IG-LSPIA), to solve problem in method. IG-LSPIA starts initial blending function, where control coefficien...