I am trying to obtain the expression for the scale factor $a(t)$ in the Einstein-de Sitter universe, that is, in a universe with a single fluid component, which corresponds to matter. To do that, I start with the Friedmann equation, simplified for the case of a single fluid component, which is the following:
$$\dfrac{\dot{a}}{a}=H_0\sqrt{\Omega_i}\ a^{-\frac{3}{2}(1+w_i)}$$
where $\Omega_i\equiv\Omega_{i,0}=\rho_{i,0}/\rho_{crit,0}$ is the density of the fluid $i$ today divided by the critical density today, considering $i=r,m,\Lambda$ (corresponding respectively to radiation, matter and dark energy). In this case we take $i=m$ (a matter dominated universe), for which $w_i\equiv w_m=0$. I am trying to reproduce the result given in the book Cosmology by Daniel Baumann as the solution to this equation, which is the following:
$$a(t)\propto t^{2/3}$$
Since in this context we assume the condition $a(t_0)=1$, the book mentions later on that, in fact:
$$a(t)=\bigg(\dfrac{t}{t_0}\bigg)^{2/3}$$
However, I see no way to arrive to this conclusion. My calculations are the following:
$$\dfrac{\dot{a}}{a}=H_0 a^{-3/2}\ \ \ \Rightarrow\ \ \ \dot{a}a^{1/2}=H_0\ \ \ \ \Rightarrow\ \ \ \dfrac{da}{dt}a^{1/2}=H_0\ \ \ \Rightarrow\ \ \ \int a^{1/2}da=\int H_0dt\ \ \ $$
$$\Rightarrow\ \ \ \dfrac{a^{3/2}}{3/2}=H_0t+C\ \ \ \Rightarrow a(t)=\bigg(\dfrac{3}{2}H_0t+\dfrac{3}{2}C\bigg)^{2/3}=\bigg(\dfrac{3}{2}H_0t+C_0\bigg)^{2/3}$$
If I impose the condition $a(t_0)=1$, what I obtain is:
$$a(t_0)=1\ \ \ \Rightarrow\ \ \ 1=\bigg(\dfrac{3}{2}H_0t_0+C_0\bigg)^{2/3}\ \ \ \Rightarrow\ \ \ \dfrac{3}{2}H_0t_0+C_0=1\ \ \ \Rightarrow\ \ \ C_0=1-\dfrac{3}{2}H_0t_0$$
Therefore, the final result would be:
$$a(t)=\bigg(\dfrac{3}{2}H_0t-\dfrac{3}{2}H_0t_0+1\bigg)^{2/3}$$
But this is not what the book proposes, and honestly, it looks ugly. I have no idea what is happening Is there any way to obtain what the book states, or maybe the formulas in the book are just approximations used in cosmology? They seem to be neglecting the additive constant, and I am puzzled as to how I can get a multiplicative integration constant instead of an additive one, so that I have something of the form $a(t)=(Ct)^{2/3}$ to later get $C=1/t_0$ by imposing the condition $a(t_0)=1$.
Any help would be greatly appreciated to solve this mystery!