else
c=-1/sqrt(1+t^2); s=-c*t; end p=A (i,:); q=A (j,:);
A (i,i)=p(i)*c^2+q (j)*s^2+2*p(j)*s*c; A (j,j)=p(i)*s^2+q (j)*c^2-2*p(j)*s*c; A (i,j)=(q (j)-p(i))*c*s+p(j)*(c^2-s^2); A (j,i)=A (i,j); for t=1:n if t~=i&t~=j
A (i,t)=p(t)*c+q (t)*s; A (t,i)=A (i,t); A (j,t)=q (t)*c-p(t)*s; A (t,j)=A (j,t); end end
p=R(:,i); q=R(:,j); for t=1:n
R(t,i)=p(t)*c+q (t)*s; R(t,j)=-p(t)*s+q (t)*c; end if m>200 %break; end end v=A (1,1); for i=2:n v=[v A (i,i)]; end toc
31