A new variational approach is proposed to calculate scattering phase shifts. In this approach, the variationalparameter is contained not explicitly in the trial function but in the boundary conditions that the scattering wavefunctionshould satisfy. It was found that a high accuracy could be reached. Furthermore, the scattering wavefunction appearsas a polynomial of finite orders in the interior region. This expression is convenient for applications.