The theoretical computation of dislocation theory in a given earth model is necessary in the explanation of observations of the co- and post-seismic deformation of earthquakes. Without using reciprocity theorem, we follow the straightforward numerical integration of co-seismic deformation given by Sun et al. to present astraight forward numerical inverse Laplace integration method. This method is used to compute the co- and post-seismic displacement of point dislocations buried in SNREI earth model and is easy to extended to theap plication of geoid and gravity. Comparing with pre-existing method, this method is relatively more straight forward and time-saving, mainly because we sum associated Legendre polynomials and dislocation love numbers before using Riemann-Merlin formula to implement SNIL and use some properties of Green functions in Laplace domain.