function [c,s]=givens(A)
%Givens變換得到變換矩陣的參數
if A(2)==0
    c=1;
    s=0;
elseif abs(A(2))>=abs(A(1))
    t=A(1)/A(2);
    s=(1+t^2)^(-1/2);
    c=s*t;
else
    t=A(2)/A(1);
    c=(1+t^2)^(-1/2);
    s=c*t;
end
end