We considered the Harris-Nesbet variational method for an
accurate calculation of phase shifts in positron-He+
collisions. Calculations have been carried out for L = 0, 1, 2,
3, 4, 5, 6 partial waves. Very accurate results of phase shift
have been obtained for the first time. The accuracy has been
double-checked by carrying out the calculations with the use of
different large-coupling schemes. We compare our results with
those available in the literature.