Highly accurate calculation of rotating neutron stars