R-mode oscillations of rapidly rotating Newtonian stars: A new numerical scheme and its application to the spin evolution of neutron stars