In this paper, the influence of errors which arise in matrix multiplications on theaccuracy of approximate solutions generated by the Block BiCGSTAB method is analyzed.In order to generate high accuracy solutions, a new Block Krylov subspacemethod is also proposed. Some numerical experiments illustrate that high accuracysolutions can be obtained by using the proposed method compared with the BlockBiCGSTAB method.