clear all; close all;

A=[1 2; 2 0];
b=[3; 0];


[L,U]=lu(A)
for j=0:2:6
    b(2)=j;
    x(j/2+1,:)=(U\(L\b)).'
end