A full-vectorial finite element method is developed to analyze the surface waves propagating at the interface between two media which could be dissipative particularly. The dissipative wave possessing a complex-valued propagation constant can be determined precisely for any given propagation direction and thus the property of losses could be thoroughly analyzed. Besides, by applying a special characteristic of the implicit circular block matrix, we can greatly reduce the computational consumptions in the analysis. By utilizing this method, the Dyakonov-like surface wave (DLSW) at the interface between a dielectric and a metal-dielectric multilayered (MDM) structure is discussed. At first, we consider the case when the involved MDM structure serves as an elliptic medium according to the effective medium approximation (EMA). Different from the conventional Dyakonov surface waves, we find that this kind of DLSW possesses an unexpected leaky property due to an additional hyperbolic-like wave in the MDM structure, resulting in a significant increase of propagation loss compared to the results estimated by a simple effective model based on the EMA. Moreover, such leaky property is found to be sensitive to the period of the MDM structure. Thus, to diminish this non-negligible leaky loss, one can suppress the amplitude of the leaky component by designing the MDM structure with a shorter period. As to the case when MDM structure serves as a hyperbolic medium, we find the calculated results for the DLSW with a small period of the MDM structure will show a great inaccuracy if we ignore the metallic absorption. However, for cases of longer periods, the influence of the metallic absorption for the DLSW becomes slight. In addition, for this DLSW, we find a trade-off between its propagation loss and the field confinement. Its propagation loss is smaller for the longer period of the MDM structure but the field becomes less confined to the propagating interface of the DLSW. Furthermore, as this hyperbolic-like MDM structure with serious nonlocal effect sometimes can support an additional elliptic-like isofrequency contour, we also discuss the DLSW based on this contour. For such DLSW, an apparent leaky property is observed similarly to the case when the MDM structure serves as an elliptic medium. This DLSW propagates with a wider range of propagation direction but suffers from a poor field confinement to the interface it is propagating along.