function [KT,K,M]=rigidez(npis,ncol,EY,hpis,mpis,dcol,k);

M=zeros(npis,npis);
K=zeros(npis,npis);
EY=EY*1000000;
for i=1:npis
    
    KT(i)=ncol*(12*EY*k(i))/(hpis(i)^3);
end

cont=npis-1;
K(1,1)=KT(1)+KT(2);
K(1,2)=-1*KT(2);
K(npis,cont)=-1*KT(npis);
K(npis,npis)=KT(npis);
for j=2:cont
    K(j,j-1)=-1*KT(j);
    K(j,j)=KT(j)+KT(j+1);
    K(j,j+1)=-1*KT(j+1);
end
for i=1:npis
    M(i,i)=mpis(i);
end