function [K]=kzong(k0,dn)
%求解總剛K---輸入變量：層間剛度矩陣k0
%                      結構層數dn
%            返回變量：總剛K
K=zeros(dn);
for i=1:dn-1
    K(i,i)=k0(i)+k0(i+1);
    K(i,i+1)=-k0(i+1);
    K(i+1,i)=-k0(i+1);
end 
K(dn,dn)=k0(dn);