In the present paper we address the interaction of pairs of charmed mesons with hidden charm in a finite box. We use the interaction from a recent model based on heavy-quark spin symmetry that predicts molecules of hidden charm in the infinite volume. The energy levels in the box are generated within this model, and from them some synthetic data are generated. These data are then employed to study the inverse problem of getting the energies of the bound states and phase shifts for D (D) over bar or D*(D) over bar*. Different strategies are investigated using the lowest two levels for different values of the box size, and the errors produced are studied. Starting from the upper level, fits to the synthetic data are carried out to determine the scattering length and effective range plus the binding energy of the ground state. A similar strategy using the effective range formula is considered with a simultaneous fit to the two levels-one above and the other one below the threshold. This method turns out to be more efficient than the previous one. Finally, a method based on the fit to the data by means of a potential and a conveniently regularized loop function, turns out to be very efficient and allows us to produce accurate results in the infinite volume starting from levels of the box with errors far larger than the uncertainties obtained in the final results. A regularization method based on Gaussian wave functions turns out to be rather efficient in the analysis and as a byproduct a practical and fast method to calculate the Luscher function with high precision is presented.