An algorithm for computing the Arf closure of an algebroid curve with more than one branch