%Replicate in Matlab the computations performed in Problem 4 of Workshop 1.
%0
%Define:
%
A = eye(6);
f = [1, 1, 1, 1, 1, 2]’;
BarA = A + f*f’
%should give
%
%BarA =
%
% 2 1 1 1 1 2
% 1 2 1 1 1 2
% 1 1 2 1 1 2
% 1 1 1 2 1 2
% 1 1 1 1 2 2
% 2 2 2 2 2 5
%
%Compute the inverse of BarA:
S = 1 + f’ * f;
InvBarA = A – f * 1./S * f’;
%Define:
b = [12, 10, 10, 10, 10, 19]’;
%and compute:
x = InvBarA * b
%which should give
%
%x =
%
% 3.0000
% 1.0000
% 1.0000
% 1.0000
% 1.0000
% 1.0000