In this paper we analyze the solution of Borderd Almost Block Diagonal (BABD) linear systems arising from the discretization of BVPs with nonseparated boundary conditions using some well known numerical methods, for example spline collocation at Gaussian points. Several approaches based on cyclic reduction are proposed and compared