Starting from the primitive basis functions, you can optimize the radial shape variationally so that the accuracy can be increased. See the details in the Section 'Orbital optimization'.