A = diag(0:49) + tril(ones(50));
b=ones(50,1);
x=A\b