  We give sufficient conditions for global existence of positive mild solutions
for the weak coupled system: \begin{eqnarray*} \frac{\partial u_{1}}{\partial
t}
&=&\rho_{1}t^{\rho_{1}-1}\Delta_{\alpha_{1}}u_{1}+t^{\sigma_{1}}u_{2}^{\beta_{1}},\
\ u_{1}\left(0\right) =\varphi_{1}, \\ \frac{\partial u_{2}}{\partial t}
&=&\rho_{2}t^{\rho_{2}-1}\Delta_{\alpha_{2}}u_{2}+t^{\sigma_{2}}u_{1}^{\beta_{2}},\
\ u_{2}\left(0\right) =\varphi_{2}, \end{eqnarray*} where $\Delta_{\alpha_{i}}$
is a fractional Laplacian, $0<\alpha_{i}\leq 2,\ \beta_{i}>1,\ \rho_{i}>0,\
\sigma_{i}>-1\ $are constants and the initial data $\varphi_{i}$ are positive,
bounded and integrable functions.
