function KM=SUBFunKMatrix(L,n,d1,d2,E)
% n--------將桿分為n段
% L--------為桿長
% P--------鋼管的密度
%d1--------鋼管的內直徑
%d2--------管的外直徑
%E---------彈性模量2.1*10^11
Le=L/n;      
r1=d1/2;      %鋼管的內半徑
r2=d2/2;      %鋼管的外半徑
I=pi*(d2^4-d1^4)/64; 
k=3*E*I/Le^3;

KM=zeros(n,n);
KM(1,1)=2*k;
for i=2:n
   KM(i,i)=2*k;
   KM(i,i-1)=-k;
   KM(i-1,i)=-k;
end
KM(n,n)=k;