In this paper, general model analysis of the form ∑c i φ i (θ) for the fractional loss of lightf(ϑ) exhibited by a close binary system will be considered in the sense of the least-squares criterion. A general recursive method will be established for constructing the normal equations for the most useful functionsφ i (θ), and an economical storage of the method on a digital computer, with its computational steps, will also be given. Moreover, a full recursive computational algorithm for the least-squares approximation will also be established. By means of this algorithm, all the solution vectors, the variance for different orders of fit and the corresponding variance-covariance matrix could be computed once and for all and, moreover, recursively. The economical storage of the algorithm and its computational steps will be given. Finally, some of the practical difficulties encountered in the application of the least-squares criterion will be analysed, and some techniques for detecting and controlling these difficulties are also given. Numerical examples on the Algol system using a Fourier cosine series of the form ∑c i cos [(j - 1) π θ/η] will be given for illustration.