  function P2()
% Checks out SymSolve
clc
n = 12;
r = 5;
A = hilb(n);
C = rand(n,r);
b = [A C ; C' zeros(r,r)]*ones(n+r,1); f = b(1:n); g = b(n+1:n+r);
x = [A C ; C' zeros(r,r)]\[f;g];
[y z] = SolveSym(A,C,f,g);
fprintf('       via backslash    via SolveSym(A,C,f,g)\n')
fprintf('-----------------------------------------------\n')
disp([x [y;z]])


