The recent development of a meshless method by using radial basis functions will be reported in this talk. Application to both multivariate interpolation and solving partial differential equations have demonstrated the spectral convergence of the