%該程序為由強度屈服系數轉換為延性系數子程序;
function [uc]=u_ru(Ru,T,Tc)
nu=100;
%Ru=2;
To=zeros(1,2); 
To(1)=0.65*Tc*Ru^0.3;%%%初始值假定T>To,即 uc=Ru;
if(T>To(1))
    uc=Ru;        
else%%%T<=To
    for j=1:nu
        tuc=(Ru-1)*To(1)/T+1;       
        To(2)=0.65*Tc*tuc^0.3;
        if(abs((To(2)-To(1))/To(1))<0.001)
            uc=tuc;
        else
            To(1)=To(2);
        end
    end
end