Initial gibbs method for orbit determination