The paper presents a numerical method for modelling of an unsteady incompressible flow of viscous fluid in a three-dimensional domain. The characteristic feature of the numerical method is nice coupling of the projection method and the operator splitting method into one numerical algorithm. The numerical method was found as a very attractive way of solving equations of the fluid motion since a fast direct method for solving a three-diagonal matrix was used in the computation.