A new procedure to best-fit earthquake magnitude probability distributions: including an example for Taiwan