In this paper, the phase-shifts for neutron-dueteron (n-d) scattering have been determined using the molecular Morse potential as theoretical model of interaction. The Triton (n-d) 2S1/2 ground state initially has been chosen as -7.61 MeV to determine the model parameters using variational Monte-Carlo technique in combination with matrix methods numerical approach to solving the time independent Schrodinger equation (TISE). The obtained potential is incorporated into the phase function equation, which is solved using Runge-Kutta (RK) 4,5 order technique, to calculate the phaseshifts at various lab energies below 15 MeV, for which experimental data is available. The results have been compared with those obtained using another molecular potential named Manning-Rosen (MR) and have been observed to fare better. Finally, the Triton ground state has been chosen as its binding energy (BE), given by -8.481795 MeV, as determined from experimental atomic mass evaluation data and the calculations are repeated. It has been found that these phase-shifts from BE data are slightly better matched with experimental ones as compared to those obtained using -7.61 MeV ground state for Triton (n-d two-body system) modeled using Morse potential.