% Test Script P1
% Checks out iSkewSolve
clc
randn('seed',0)
n = 100; b = randn(n,1); c = randn(n-1,1); c(1:2:n-1) = 1000000*c(1:2:n-1);
format long
[x,relerr] = iSkewSolve(b,c);
SnapShot = x(1:10:n)
RelErrEstimate = relerr
A = eye(n,n) + diag(c,1) + diag(-c,-1);
x0 = A\b;
ConditionOfA = cond(A,1)
ActualRelerr = norm(x-x0)/norm(x0)

