[K0]=MatrixK0(Kape,Kapc,XYZ)

function [K0]=MatrixK0(Kape,Kapc,XYZ)
l1=3232.08;
l2=2740;
l3=2262.74;
l4=1812.07;
l5=1341.64;
d=2.5;
a=pi*(d/2)^2;
Eac=2100;
c=a*Eac*2;
%rigideces cables
Kc1=c/l1;
Kc2=c/l2;
Kc3=c/l3;
Kc4=c/l4;
Kc5=c/l5;
%rigideces tablero (flexurales)
Ehor=250;
Ita=3.12*10^6;
Kt1=Ehor*Ita/100;
Kt2=0.00000001;                            
%rigideces axiales columnas
lc=80;
c2=Ehor/lc;
acs1=18150;
acs2=19800;
acs3=690*110;
acs4=110^2*2;
Ics1=18301250;
Ics2=19965000;
Ics3= 76532500;
Ics4=110^4/6;
KFcs1=c2*Ics1;
KFcs2=c2*Ics2;
KFcs3=c2*Ics3;
KFcs4=c2*Ics4;
KAcs1=c2*acs1;
KAcs2=c2*acs2;
KAcs3=c2*acs3;
KAcs4=c2*acs4;
%matriz de rigidez
KC1=[Kc1;Kc2;Kc3;Kc4;Kc5];
KC2=[Kc5;Kc4;Kc3;Kc2;Kc1];
KK=[Kt2;ones(108,1)*Kt1;Kt2;...%rigideces flexurales tablero (1-->110)
    KC1;KC2;KC1;KC2;...%rigideces axiales cables (111-->130)
    Kape;Kapc;Kapc;Kape;...%rigideces apoyos tablero (131-->134)
    KFcs4*2;ones(5,1)*KFcs4;ones(2,1)*KFcs3;ones(16,1)*KFcs2;ones(9,1)*KFcs1;0.000000001;...%rigideces flexurales columna1 (135-->168)
    KFcs4*2;ones(5,1)*KFcs4;ones(2,1)*KFcs3;ones(16,1)*KFcs2;ones(9,1)*KFcs1;0.000000001;...%rigideces flexurales columna2 (169-->202)
    KAcs4*ones(6,1);KAcs3;KAcs2*ones(16,1);KAcs1*ones(10,1);...%rigideces axiales columna1 ()
    KAcs4*ones(6,1);KAcs3;KAcs2*ones(16,1);KAcs1*ones(10,1)];%rigideces axiales columna2 ()
K0=diag(KK);
return