clear all
%
% Example 02 - random matrices
%
  tic
  m=5
  n=25
  A = rand(m,m,n);
  B = rand(m,m,n);   
  C = rand(m,m,n);
  D = rand(m,m,n);
  E = rand(m,m,n);
  F = zeros(m,n);
  X = zeros(m,n);
  R = zeros(m,n);
%
% Form matrices
%
  for k = 1:m
      C(k,k,1:n) = C(k,k,1:n) + 100.0;
  end
%
% Run example
%
  Xe = ones(m,n);
  F = ResCBPS(A,B,C,D,E,R,Xe);
  toc
  tic
  X = CBPS(A,B,C,D,E,F);
  toc
  R = ResCBPS(A,B,C,D,E,F,X);
  err_abs = max(max(abs(X-Xe)))  
  err_rel = err_abs/max(max(Xe))