If you look at the Lagrangian of a point particle in the gravitational field of a massive star:
<br />
L = \frac{1}{2} \, m \, v^2 - \frac{G \, M \, m}{r}<br />
and you scale M \rightarrow \mu \, M, \mathbf{r} \rightarrow \lambda \, \mathbf{r}, and t \rightarrow \tau \, t, then your Lagrangian changes as:
<br />
L \rightarrow \left( \frac{\lambda}{\tau} \right)^2 \frac{1}{2} \, m \, v^2 - \frac{\mu}{\lambda} \, \frac{G \, M \, m}{r}<br />
If you then impose:
<br />
\left( \frac{\lambda}{\tau} \right)^2 = \frac{\mu}{\lambda} = 1 \Rightarrow \frac{\lambda^3}{\tau^2} = \mu<br />
you see that the Lagrangian gets multiplied by a common factor. This means that if \mathbf{r}(t) was a possible equation of a trajectory, then, so is \lambda \, \mathbf{r}(\tau t) for a system with a gravitational center with a \mu times bigger mass.
If the system loses its mass "slowly", then the trajectory gets modified in the above scaling fashion. There is, however, another relation between the scaling factors. The total mechanical energy of the particle is not conserved, because the loop integral of the gravitational force over one period is non-zero, since it changes with time.
Angular momentum, however, is conserved, because the force is always radially directed. Angular momentum scales as:
<br />
\mathbf{M} = m (\mathbf{r} \times \mathbf{v}) \rightarrow \frac{\lambda^2}{\tau} \, \mathbf{M}<br />
so, we have the further constraint:
<br />
\frac{\lambda^2}{\tau} = 1<br />
These two conditions determine \tau, and \lambda in terms of \mu:
<br />
\tau = \lambda^2 \Rightarrow \frac{\lambda^3}{\lambda^4} = \mu \Rightarrow \lambda = \mu^{-1}, \ \tau = \mu^{-2}<br />
This means that if the mass of the star changes by a fractional amount \epsilon \ll 1, then:
<br />
\epsilon \equiv \frac{\Delta M}{M} = \frac{ M_\mathrm{new} - M_\mathrm{old} }{ M_\mathrm{old} } = \frac{ M_\mathrm{new} }{ M_\mathrm{old} } - 1 = \mu - 1 \Rightarrow \mu = 1 + \epsilon<br />
Then, notice that according to the binomial theorem:
<br />
(1 + \epsilon)^{n} \sim 1 + n \, \epsilon, \ \epsilon \ll 1<br />
Thus:
<br />
\lambda = (1 + \epsilon)^{-1} \sim 1 - \epsilon<br />
and
<br />
\tau = (1 + \epsilon)^{-2} \sim 1 - 2 \, \epsilon<br />
So, we quantify the scaling as, if the mass of the star has a fractional decrease \epsilon, then the linear dimensions of a planet's orbit have a fractional increase \epsilon, and the orbital period has a fractional increase 2 \epsilon.