Initially the derivation is given of a simple iterative scheme for the calculation of linear best Lp approximations over either a continuum or a discrete point set. The scheme, whilst each iteration is computationally efficient, proves to be divergent. An analysis of the propagation of errors leads to the development of an equally efficient scheme with second order convergence. A further slight modification is shown to ensure convergence.