function [y,count] = Diagonal(D,b)
n = size(b);
count = 0;
for i = 1:n
    count = count+1;
    y(i,1) = b(i)/D(i,i);
end