The hysteresis behaviour at the L-H-L transitions in tokamak plasma is investigated based on bifurcation concept. The formation of an edge transport barrier (ETB) is modeled via thermal and particle transport equations with the flow shear suppression effect on anomalous transport included. The anomalous transport is modeled based on critical gradients threshold and the flow shear is calculated from the force balance equation, couples the two transport equations leading to a non-linear behaviour. Analytical investigation reveals that the fluxes versus gradients space exhibits bifurcation behaviour with s-curve soft bifurcation type. Apparently, the backward H-L transition occurs at lower values than that of the forward L-H transition, illustrating hysteresis behaviour. The hysteresis properties, i.e. locations of threshold fluxes, gradients and their ratios are analyzed as a function of neoclassical and anomalous transport values and critical gradients. It is found that the minimum heat flux for maintaining H-mode depends on several plasma parameters including the strength of anomalous transport and neoclassical transport. In particular, the hysteresis depth becomes larger when neoclassical transport decreases or anomalous transport increases.