Computation of Roche limit
Right now, the approximation used to compute the Roche limit is real roche(real dist) const { return dist*std::cbrt(myMass/3); }
.
That is, the star mass is not explicitly taken into account nor G. I suspect it's because they're both set to 1 but I'd like to make them explicit.