Doppler observation of extrasolar planets through the radial velocity displacement of their host stars can only determine lower limits of planetary masses. We numerically integrate ν Andromedae planetary orbits with various initial conditions of masses and angle variables to investigate which initial configuration produces stable orbits during the timescale of host star's age. According to our preliminary results starting from Lick dataset, ν Andromedae planetary system seems to remain stable over the timescale of its host star's age if sin i > 0.7 where i is the unknown line-of-sight inclination of planetary orbits. In this case we may estimate that the upper limit masses of ν Andromedae planets in our model is about 1/0.7 ∽ 1.43 times larger than its minimum.