This method was described by Estabrook et al. in 1985.

     How to run: the value of the parameter '-m' must be 'all' or 'quartest' or 'triplets' to run the quartets and/or triplets method.