This paper is concerned with the problem of determining the time history of position, velocity and thrust acceleration of a space vehicle given only angular observations at discrete times and without prior knowledge of a nominal trajectory. It is shown that the motion, which involves unmeasurable state variables and partially unknown system dynamics, can be mathematically modeled by a system of ordinary differential equations subject to distributed boundary conditions. The resulting system is then solved by iteration with a modified quasilinearization technique. Results obtained from application of the method to a logarithmic spiral orbit will be examined. It will be shown that the method exhibits high precision and fast convergence even with poor estimates of the initial values for the state vector that is used to start the iterative process. ONTINUOUSLY thrusted interplanetary vehicles of the future will require new techniques for their guidance and navigation. At the present time much research has been done on trajectory analysis and optimal control for these missions, yet little work has been done on the critical problem of orbit determination. Jordan1 has shown that it is possible to solve the problem of determining small perturbations in the vehicle's position, velocity, and acceleration by applying filtering techniques to the data if one has prior knowledge of the vehicle's nominal trajectory. Apparently no one, however, has solved the more difficult problem of determining the completely unknown position, velocity, and thrust program of a powered vehicle from the types of observations which are now available; namely telescopic (angular observations) and radar (range and range-rate), without knowledge of a good approximate nominal. In this paper it will be shown that the motion of a thrusted vehicle, which involves unmeasurable state variables and partially unknown system dynamics can be mathematicall y modeled by a system of ordinary differential equations subject to multipoint boundary conditions. Prior knowledge of the vehicle's nominal trajectory is not necessary for a solution. Furthermore, the method of solution is not necessarily restricted to the class of problems where the thrust acceleration magnitude is small compared to the gravitational acceleration. The analytic method of differential approximation and the numerical technique of quasilineariz ation are the basic mathematical tools which will be used to solve the problem. First, the unknown portion of the system dynamics, that which describes the thrust program, is represented by the method of differential approximation which was originally introduced by Bellman 2 for the purpose of state estimation, and extended to the present application by Carpenter.3 Upon application of differential approximation the complete system dynamics can be described by a set of ordi