Numerical modeling of long Josephson junctions in the frame of double sine-gordon equation