In summary, we discussed the description of a Bose-Einstein condensate using a wave function, where the functions \phi(x,t) and \rho(x,t) are real. We then explored the probability density and current density of the condensate and wrote down the continuity equation in terms of these functions and their derivatives. Finally, we substituted the wave function into the Schrödinger equation and simplified it, comparing the resulting equation to the continuity equation.
Oh. Your previous expression was actually incorrect. The last 2 terms in the rhs should have contained single derivatives of both rho and phi, and not double derivativs (compare with the expression prior to that).