The objective of this study was to devise an optimum force system to achieve en-masse retraction of six maxillary anterior teeth in lingual orthodontics (LiO). First, the set of equations was developed based on the mathematical computation to estimate optimum parameters of force system. Then, the computer software based on this mathematical computation was developed for the ease of estimation of force system. The verification of force system obtained with computer software was accomplished by three-dimensional finite element analysis (FEA). In FEA, it was clear that the desired en-masse retraction of six maxillary anterior teeth in LiO was achieved as observed from the vectors of nodal displacements as well as positions of undeformed and deformed models. In this way, mathematically computed optimum force system was verified with FEA. For orthodontists, the developed computer software accurately estimates the required force system according to mathematical computation. It simplifies the task of computation for orthodontists. The orthodontists can easily operate this computer software as it is user-friendly. Further, in-vivo validation of this study is required in future before clinical application.