In recent years, there are attempts to describe quantization of planetary distance based on time-independent gravitational Schr¨odinger equation, including Rubcic & Rubcic’s method and also Nottale’s Scale Relativity method. Nonetheless, there is no solution yet for time-dependent gravitational Schr ¨odinger equation (TDGSE). In the present paper, a numerical solution of time-dependent gravitational Schrodinger equation is presented, apparently for the first time. This numerical solution leads to gravitational Bohr-radius, as expected. In the subsequent section, we also discuss plausible extension of this gravitational Schr¨odinger equation to include the effect of phion condensate via Gross-Pitaevskii equation, as described recently by Moffat. Alternatively one can consider this condensate from the viewpoint of Bogoliubov de Gennes theory, which can be approximated with coupled time-independent gravitational Schr¨odinger equation. Further observation is of course recommended in order to refute or verify this proposition.