@s3alfisc

It seems fairly straightforward to implement this. Most of fun is done in one function cov in this function. I'm not clear though on what the bootalgo function returns. Is it supposed to be a test statistic or a covariance matrix? If the latter, then we can just pop the function in and we should be done.

Also, don't merge it yet as I think we should merge this in once boot_algo3 and other associated pull requests are merged first