We solve numerically Eqs. (REF ) by split-operator technique [1]} for trajectories of an atomic white dwarf corresponding to closed eccentric orbits. We consider the Bose-Fermi droplet consisted of \(^{133}\) Cs bosonic and \(^{6}\) Li fermionic atoms. Such mixtures are studied experimentally nowadays [2]}, [3]}, [4]}. The droplet, consisted of 1460 bosonic and 100 fermionic atoms, is initially located at apostron at some distance (about \(7\, r_{\rm wd}\) ) from the artificial black hole, far away from the horizon. The initial numerical parameters are chosen in such a way that we can follow several revolutions of a white dwarf in an eccentric orbit (with periastron at about \(4\, r_{\rm wd}\) ). Such a task is numerically manageable (i.e., is not too intensive with respect to the computational time) provided the white dwarf is not located extremely far from the black hole. This, on the other hand, means that the period of revolution of a white dwarf in our case remains short, in fact it is at most of the order of ten seconds. At each periastron passage a white dwarf is stripped off an approximately (at least within initial few passes) the same mass, which is about \(10^{-3}\)  of the white dwarf mass for a particular orbit discussed below.
