function y=r(A,a,i,j)
A(j,:)=A(j,:)+a*A(i,:);
y=A;