It can be shown that the corresponding (modified) FR equations now read [1]},
\(H^2&=&\frac{8\pi G \rho }{3}\left(1+\gamma ^2 \frac{\rho }{\rho _c}\right) \left(1-\frac{(\gamma ^2+1)\rho }{\Delta ^2\rho _c}\right), ~~~\\\frac{\ddot{a}}{a}&=&-\frac{4\pi G}{3}\left(\rho +3P\right)-\frac{4\pi G P\rho }{\Delta \rho _c} \left[3\left(\gamma ^2+1\right)-2\Delta \right]\nonumber \\&&-\frac{4\pi G \rho ^2}{3\Delta ^2\rho _c}\left[7\gamma ^2-1 +\left(5\gamma ^2-3\right)\left(\Delta -1\right) -\frac{4\gamma ^2\rho }{\rho _c}\right],\nonumber \\\) 
