..., a polynomial-based partition of unity method is developed, in which locally-constructed least-squares polynomial approximations are smoothly blended to obtain a globally-accurate solution. Comprehensive numerical experiments demonstrate the robustness of the method and confirm its superior accuracy, stability...