Metastatic papillary thyroid cancer under treatment with different I-131 protocols: a mathematical model