The variability of memristive devices using multilayer hexagonal boron nitride (h-BN) coupled with Ti and Au electrodes (i.e., Au/Ti/h-BN/Au) is analyzed in depth using different numerical techniques. We extract the reset voltage using three different methods, quantify its cycle-tocycle variability, calculate the charge and flux that allows to minimize the effects of electric noise and the inherent stochasticity of resistive switching, describe the device variability using time series analyses to assess the "memory" effect, and employ a circuit breaker simulator to understand the formation and rupture of the percolation paths that produce the switching. We conclude that the cycle-to-cycle variability of the Au/Ti/h-BN/Au devices presented here is higher than that previously observed in Au/h-BN/Au devices, and hence, they may be useful for data encryption.