complete pivoting

complete

clear

clc

function x = completepivot(A,b)

n = length(b);

Ab = [A b];

for k = 1:n-1

maxval = 0;

p = k;

q = k;

for i = k:n

for j = k:n

if abs(Ab(i,j)) > maxval

maxval = abs(Ab(i,j));

p = i;

q = j;

end

end

end

temp = Ab(k,:);

Ab(k,:) = Ab(p,:);

Ab(p,:) = temp;

temp = Ab(:,k);

Ab(:,k) = Ab(:,q);

Ab(:,q) = temp;

for i = k+1:n

m = Ab(i,k)/Ab(k,k);

Ab(i,k:n+1) = Ab(i,k:n+1) - m*Ab(k,k:n+1);

end

end

x = zeros(n,1);

x(n) = Ab(n,n+1)/Ab(n,n);

for i = n-1:-1:1

sum = 0;

for j = i+1:n

sum = sum + Ab(i,j)*x(j);

end

x(i) = (Ab(i,n+1)-sum)/Ab(i,i);

end

end