Implement confidence bands

[ManuelNavarroGarcia]

Include the option to compute confidence bands via a non-parametric bootstrap procedure, which entails an estimation of the probability distribution from the data and a sampling from that distribution. This option must be passed through the .fit() method or other function must be defined for this purpose. Moreover, the definition must be flexible enough so that some parameters as the number of repetitions may be tuned by the user.