An r-coloring of a set S is a function from S to \(\{0,1, \cdots , r-1\}\) . Let \((E_n)\) be the equation \(x_1+x_2+\cdots +x_n=y^2\) . The r-color Rado number \(R_r(E_n)\) of \((E_n)\) is the least integer R, provided it exists, such that every r-coloring of \(\{1,2,\cdots , R\}\) admits a monochromatic solution to \((E_n)\) . In this study, for each \(n,\ r\ge 2\) , we provide a lower bound \(\ell _r\) of \(R_r(E_n)\) . Accordingly, we show that when \(2\le n\le 6\) , \(R_2(E_n)=n\) , and when \(n\ge 7\) , \(R_2(E_n)=\Big \lceil \sqrt{n\lceil { {\sqrt{n}}} \rceil }\Big \rceil \) , which coincides with the lower bound \(\ell _2\) that we have provided. We also show that \(R_3(E_n)=n\) when \(2\le n\le 11\) , and \(R_3(E_{12})=11\) .