If the force field acting on an artificial Earth satellite is not known a priori with sufficient accuracy to represent its observations on their accuracy level, one may introduce so-called pseudo-stochastic parameters into an orbit determination process, e.g. instantaneous velocity changes at user-defined epochs or piecewise constant accelerations in user-defined adjacent time subintervals or piecewise linear and continuous accelerations in adjacent time subintervals. The procedures, based on standard least-squares, associated with such parameterizations are well established, but they become inefficient (slow) if the number of pseudo-stochastic parameters becomes large. We develop two efficient methods to solve the orbit determination problem in the presence of pseudo-stochastic parameters. The results of the methods are identical to those obtained with conventional least-squares algorithms. The first efficient algorithm also provides the full variance-covariance matrix; the second, even more efficient algorithm, only parts of i
