  function P3()
% Checks out TriMatSol
clc
n = 5;
S = eye(n,n) + tril(randn(n,n),-1);
T = eye(n,n) + triu(randn(n,n),1);
B = randn(n,n);
X = TriMatSol(S,T,B);
R = S*X*T + X - B;
fprintf('\n|| SXT + X - B || = %10.3e\n\n',norm(R,1))
