clear;clc
a=tril(magic(6),0);
b=1./tril(a,-1)';
b(b==inf)=0;
sym(a+b) % a+b