V. CONCLUSIONS AND DISCUSSIONS In this paper we have investigated the dynamical behavior of the EoS of DE in the region of low redshift in a nearly model-independent way. The redshift in that region is binned and w$_{de}$ in each bin is approximated by a linear expansion of redshift z , and in the large redshift region we set w$_{de}$ to be a constant w$_{L}$ . While fitting the model with some observational data which include SnIa, BAO, CMB and Hubble evolution data, we leave the divided points of bins as free parameters. If the evolution of w$_{de}$ is not monotonous, or is not linear enough in the region under consideration, the best-fitted divided points will represent the turning points, where w$_{de}$ changes its evolving direction significantly. In this way we can explicitly reconstruct w$_{de}$ by using a few bins, and the errors of parameters from observational data will be small due to the small number of bins. First we have tried to find the turning points within the region of redshift z ∈ (0 , 1 . 8), and set w$_{L}$ = − 1 in the region z ∈ (1 . 8 , ∞ ) (Model I). Our results show that the data favor two turning points of w$_{de}$ in z ∈ (0 , 1 . 8), and w$_{de}$ may have an oscillation form [30]. Our results are consistent with those by the UBE method in [27]. Since the main data points are in z ∈ (0 , 1) and our result in Model I shows there may be a turning point around z ∼ 1, to see clearly the dynamical behavior of EoS in that region, we have focused on the region z ∈ (0 , 0 . 9) in Model II. We have found one turning point only in z ∈ (0 , 0 . 9), the reconstructed w$_{de}$ in the best-fitted model is almost the same as that reconstructed in Model I in z ∈ (0 , 0 . 9). We have also obtained the errors of w$_{de}$ at 1 σ and 2 σ in z ∈ (0 , 0 . 9). In both correlated and uncorrelated estimates with a fixed w$_{L}$ = − 1 or a floating constant w$_{L}$ , we found that there is a 2 σ deviation of w$_{de}$ from − 1 around z = 0 . 9. It is interesting to see whether the deviation of EoS from − 1 around z = 0 . 9 is physical, or is caused by some unknown technical causes in fitting. If it is physical, it then clearly shows that DE is dynamical. But in UBE of w$_{de}$ there seems no such distinct deviation around z = 0 . 9, it may be due to the difference between the discontinuity of w$_{de}$ in the piecewise constant case and the continuity in LS case [31]. In [16], where the cubic-spline method is used, there is also no such an explicit deviation around z = 0 . 9, but it is likely due to its set of EoS in the last bin w$_{L}$ = w (1): to fit well with the data of z > 1, w (1) should be much minus, which would suppress the reconstructed w$_{de}$ around z = 1. Of course, it is also possible that such a big deviation around z = 0 . 9 is due to the non-smoothness of 13