In this paper, a fitting method for arbitrary multivariate scattered data is presented.The method consists of uniforming the set of arbitrary 4dimensional scattered data byusing local weighted least square, evaluating the values of partial derivatives and minedpartial derivatives over the cube-gridded-points only by using those gridded-point data,followed by interpolating smoothly and quickly those gridded-point data. Good resultshave been gotten by computing several known functions.