Mathematical and numerical analysis of a model for anti-angiogenic therapy in metastatic cancers