The paper describes three algorithms for the construction of initial approximations for the solution of systems of linear algebraic equations by the preconditioned conjugate gradient method. The considered linear systems arise from a finite difference approximation to selfadjoint boundary value problems and the initial approximations are constructed making use of different orthogonal projectors. Algorithms are tested on a model problem which is taken from modelling of semiconductors.