The velocity distribution function (VDF) of the hypothetical Weakly Interacting Massive Particles (WIMPs), currently the most favored candidate for the Dark Matter (DM) in the Galaxy, is determined directly from the circular speed ("rotation") curve data of the Galaxy assuming isotropic VDF. This is done by "inverting" -using Eddington's method -the Navarro-Frenk-White universal density profile of the DM halo of the Galaxy, the parameters of which are determined, by using Markov Chain Monte Carlo (MCMC) technique, from a recently compiled set of observational data on the Galaxy's rotation curve extended to distances well beyond the visible edge of the disk of the Galaxy. The derived most-likely local isotropic VDF strongly differs from the Maxwellian form assumed in the "Standard Halo Model" (SHM) customarily used in the analysis of the results of WIMP direct-detection experiments. A parametrized (non-Maxwellian) form of the derived mostlikely local VDF is given. The astrophysical "g-factor" that determines the effect of the WIMP VDF on the expected event rate in a direct-detection experiment can be lower for the derived most-likely VDF than that for the best Maxwellian fit to it by as much two orders of magnitude at the lowest WIMP mass threshold of a typical experiment.Several experiments worldwide are currently trying to directly detect the hypothetical Weakly Interacting Massive Particles (WIMPs), thought to constitute the Dark Matter (DM) halo of our Galaxy, by looking for nuclear recoil events due to scattering of WIMPs off nuclei of suitably chosen detector materials in low background underground facilities. The rate of nuclear recoil events depends crucially on the local (i.e., solar neighborhood) density and velocity distribution of the WIMPs in the Galaxy [1], which are a priori unknown. Estimates based on a variety of observational data typically yield values for the local density of DM, ρ DM,⊙ , in the range 0.2 -0.4 GeV cm −3 ((0.527 − 1.0) ×10 −2 M ⊙ pc −3 ) [2]. In contrast, not much knowledge directly based on observational data is available on the likely form of the velocity distribution function (VDF) of the WIMPs in the Galaxy. The standard practice is to use what is often referred to as the "Standard Halo Model" (SHM), in which the DM halo of the Galaxy is described as a single-component isothermal sphere [3], for which the VDF is assumed to be isotropic and of Maxwell-Boltzmann (hereafter simply "Maxwellian") form, f (v) ∝ exp(−|v| 2 / v 0 2 ), with a truncation at an assumed value of the local escape speed, and with v 0 = v c,⊙ , the circular rotation velocity at the location of the Sun. Apart from several theoretical issues (see, e.g., [4]) concerning the self-consistency of the SHM * pijush.bhattacharjee@saha.ac.in † soumini.chaudhury@saha.ac.in ‡ susmita.kundu@saha.ac.in § subha@tifr.res.in as a model of a finite-size, finite-mass DM halo of the Galaxy, high resolution cosmological simulations of DM halos [5] give strong indications of significant departure of the VDF from the Maxwellian...