In this paper we present a new method of numerical evaluation of unknown coefficients of a dynamical system having available information about unknown phase trajectories at some time values. The method consists in the direct integration of given dynamical system with posterior application of a quadrature rules. Using the least square method and possible Constraints we obtain a linear system for determining an unknown coefficient. A numerical example illustrates the method.