Kinesin-3 and kinesin-1 molecular motors are two families of the kinesin superfamily. It has been experimentally revealed that in monomeric state kinesin-3 is inactive in motility and cargo-mediated dimerization results in superprocessive motion, with an average run length being more than 10-fold longer than that of kinesin-1. In contrast to kinesin-1 showing normally single-exponential distribution of run lengths, dimerized kinesin-3 shows puzzlingly Gaussian distribution of run lengths. Here, based on our proposed model, we studied computationally the dynamics of kinesin-3 and compared with that of kinesin-1, explaining quantitatively the available experimental data and revealing the origin of superprocessivity and Gaussian run length distribution of kinesin-3. Moreover, predicted results are provided on ATP-concentration dependence of run length distribution and force dependence of mean run length and dissociation rate of kinesin-3.