%================================================================================%
%3D analysis of continuous or solid body using Finite E1ement.Mbthod
%Element kind: 20-Node Iso parametric (by using l4point Gaussian rule)
%Written by: ”Sobhan Rostami & Ali.Mbeinadini"
%NLSC students of structural engineering of Azad university of Kerman-Iran
%E-mail: info@newcivil.com
%================================================================================%
clc
clear all
%Input data:
E=xlsread('Iso20node','input','F9');
v=xlsread('Iso20node','input','E10');
NN=xlsread('Iso20node','input','El2');
Coor=xlsread('Iso20node','input','Cl8:F1000');
for j=1:NN
i=Coor(j,1);
X(i)=Coor(j,2);
Y(i)=Coor(j,3);
Z(i)=Coor(j,4);
end
NE=xlsread('Iso20node','input','El3');
Dir=xlsread('Iso20node','input','M21:AG1000');
NC=xlsread('Iso20node','input','AO9');
Cons=xlsread('Iso20node','input','AJ18:AM1000');
NF=xlsread('Iso20node','input','AX9');
Forc=xlsread('Iso20node','input','AU16:AX1000');
%Definition of properties matrix:
D=E/((l+v)*(1-(2*v)))*[1-v V v 0 0 0;v 1-v v 0 0 0;v v 1-V 0 0 0;
0 0 0 ((1-(2*v))/2) 0 050 0 0 0 ((1*(2*v))/2) 0;0 0 0 0 0 ((1-(2*V))/2)];
%Gaussian quadrature points:
w1=0.335180055401662;
w2=0.886426592797784;
r(1)=-sqrt(19/33); s(1)=-sqrt(19/33); t(l)=sqrt(19/33);
r(2)=-sqrt(19/33); s(2)=-sqrt(19/33); t(2)=-sqrt(19/33);
r(3)=-sqrt(19/33); s(3)=sqrt(19/33); t(3)=-sqrt(19/33);
r(4)=-sqrt(19/33); s(4)=sqrt(l9/33); t(4)=sqrt(l9/33);
r(5)=sqrt(19/33); s(5)=-sqrt(19/33); t(5)=sqrt(19/33);
r(6)=sqrt(19/33); s(6)=-sqrt(19/33); t(6)=-sqrt(19/33);
r(7)=sqrt(19/33); s(7)=sqrt(19/33); t(7)=-sqrt(19/33);
r(8)=sqrt(l9/33); s(8)=sqrt(19/33); t(8)=sqrt(19/33);
r(9)=0; s(9)=-sqrt(38/60); t(9)=0;
r(l0)=0; s(10)=sqrt(38/60); t(10)=0;
r(11)=0; s(11)=0; t(11)=-sqrt(38/60);
r(12)=0; s(12)=0; t(12)=sqrt(38/60);
r(13)=-sqrt(38/60); s(l3)=0; t(13)=0;
r(14)=sqrt(38/60); s(14)=0; t(l4)=0;
%Definitian of strain-displacement matrix 'LB]' for each element:
for n=1:NE
e2(n)=Dir(n,2);
e6(n)=Dir(n,3);
e7(n)=Dir(n,4);
e3(n)=Dir(n,5);
e1(n)=Dir(n,6);
e5(n)=Dir(n,7);
e8(n)=Dir(n,8);
e4(n)=Dir(n,9);
e19(n)=Dir(n,l0);
e14(n)=Dir(n,11);
e20(n)=Dir(n,12);
e10(n)=Dir(n,13);
e18(n)=Dir(n,14);
e16(n)=Dir(n,15);
e17(n)=Dir(n,l6);
e12(n)=Dir(n,17);
e9(n)=Dir(n,18);
e13(n)=Dir(n,l9);
e15(n)=Dir(n,20);
e1l(n)=Dir(n,21);
for i=1:14
Xr(i)=(1/8*(X(e1(n))*(1-s(i))*(1+t(i))*(1-t(i)+s(i)+(2*r(i)))+...
X(e2(n))*(1-s(i))*(1-t(i))*(1+t(i)+s(i)+(2*r(i)))+...
X(e3(n))*(1+s(i))*(1-t(i))*(1+t(i)-s(i)+(2*r(i)))+...
X(e4(n))*(l+s(i))*(1+t(i))*(1-t(i)-s(i)+(2*r(i)))-...
X(e5(n))*(l-s(i))*(1+t(i))*(1+t(i)+s(i)-(2*r(i)))-...
X(e6(n))*(1-s(i))*(1-t(i))*(1+t(i)+s(i)-(2*r(i)))-...
X(e7(n))*(1+s(i))*(1-t(i))*(1+t(i)-s(i)-(2*r(i)))-...
X(e8(n))*(1+s(i))*(1+t(i))*(1-t(i)-s(i)-(2*r(i))))+...
1/4*(-X(e9(n))*(1-s(i))*(1-t(i)^2)-X(e10(n))*(1-s(i)^2)*(l-t(i))-...
X(e11(n))*(1+s(i))*(1-t(i)^2)-X(e12(n))*(1-s(i)^2)*(1+t(i))+...
X(e13(n))*(l-s(i))*(1-t(i)^2)+X(e14(n))*(1-s(i)^2)*(1-t(i))+...
X(e15(n))*(1+s(i))*(1-t(i)^2)+X(e16(n))*(1-s(i)^2)*(1+t(i))-...
X(e17(n))*(2*r(i))*(1+s(i))*(1+t(i))-...
X(e18(n))*(2*r(i))*(1-s(i))*(1+t(i))-...
X(e19(n))*(2*r(i))*(1-s(i))*(1-t(i))-...
X(e20(n))*(2*r(i))*(1+s(i))*(1-t(i))));
Xs(i)=(1/8*(X(e1(n))*(1-r(i))*(1+t(i))*(1-t(i)+r(i)+(2*s(i)))+...
X(e2(n))*(l-r(i))*(1-t(i))*(1+t(i)+r(i)+(2*s(i)))-...
X(e3(n))*(1-r(i))*(1-t(i))*(l+t(i)+r(i)-(2*s(i)))-...
X(e4(n))*(1-r(i))*(1+t(i))*(1-t(i)+r(i)-(2*s(i)))+...
X(e5(n))*(1+r(i))*(1+t(i))*(1-t(i)+r(i)+(2*s(i)))+...
X(e6(n))*(1+r(i))*(1-t(i))*(1+t(i)-r(i)+(2*s(i)))-...
X(e7(n))*(1+r(i))*(1-t(i))*(1+t(i)-r(i)-(2*s(i)))-...
X(e8(n))*(1+r(i))*(1+t(i))*(1-t(i)-r(i)-(2*s(i))))+...
1/4*(-X(e9(n))*(l-r(i))*(1-t(i)^2)-...
X(e10(n))*(2*s(i))*(1-r(i))*(1-t(i))+...
X(e11(n))*(1-r(i))*(1-t(i)^2)-...
X(e12(n))*(2*s(i))*(1-r(i))*(1+t(i))-...
X(e13(n))*(1+r(i))*(1-t(i)^2)-...
X(e14(n))*(2*s(i))*(1+r(i))*(1-t(i))+...
X(e15(n))*(1+r(i))*(1-t(i)^2)-...
X(e16(n))*(2*s(i))*(1+r(i))*(1+t(i))+...
X(e17(n))*(1-r(i)^2)*(1+t(i))-...
X(e18(n))*(1-r(i)^2)*(1+t(i))-...
X(e19(n))*(1-r(i)^2)*(l-t(i))+...
X(e20(n))*(1-r(i)^2)*(l-t(i))));
Xt(i)=(1/8*(-X(e1(n))*(l-r(i))*(1-s(i))*(1+s(i)+r(i)-(2*t(i)))+...
X(e2(n))*(1-r(i))*(1-s(i))*(1+s(i)+r(i)+(2*t(i)))+...
X(e3(n))*(1-r(i))*(1+s(i))*(1-s(i)+r(i)+(2*t(i)))-...
X(e4(n))*(1-r(i))*(1+s(i))*(1-s(i)+r(i)-(2*t(i)))-...
X(e5(n))*(1+r(i))*(1-s(i))*(1+s(i)-r(i)-(2*t(i)))+...
X(e6(n))*(1+r(i))*(1-s(i))*(1+s(i)-r(i)+(2*t(i)))+...
X(e7(n))*(l+r(i))*(1+s(i))*(1-s(i)-r(i)+(2*t(i)))-...
X(e8(n))*(1+r(i))*(1+s(i))*(1-s(i)-r(i)-(2*t(i))))+...
1/4*(-X(e9(n))*(2*t(i))*(1-r(i))*(1-s(i))-...
X(e10(n))*(1-r(i))*(1-s(i)^2)-...
X(e11(n))*(2*t(i))*(1-r(i))*(1+s(i))+...
X(e12(n))*(1-r(i))*(1-s(i)^2)-...
X(e13(n))*(2*t(i))*(1+r(i))*(1-s(i))-...
X(e14(n))*(1+r(i))*(1-s(i)^2)-...
X(e15(n))*(2*t(i))*(1+r(i))*(1+s(i))+...
X(e16(n))*(1+r(i))*(1-s(i)^2)+...
X(e17(n))*(l-r(i)^2)*(1+s(i))+...
X(e18(n))*(l-r(i)^2)*(1-s(i))-...
X(e19(n))*(1-r(i)^2)*(1-s(i))-...
X(e20(n))*(1-r(i)^2)*(1+s(i))));
Yr(i)=(1/8*(Y(e1(n))*(1-s(i))*(1+t(i))*(1-t(i)+s(i)+(2*r(i)))+...
Y(e2(n))*(1-s(i))*(1-t(i))*(1+t(i)+s(i)+(2*r(i)))+...
Y(e3(n))*(1+s(i))*(1-t(i))*(1+t(i)-s(i)+(2*r(i)))+...
Y(e4(n))*(1+s(i))*(1+t(i))*(1-t(i)-s(i)+(2*r(i)))-...
Y(e5(n))*(1-s(i))*(1+t(i))*(1-t(i)+s(i)-(2*r(i)))-...
Y(e6(n))*(1-s(i))*(1-t(i))*(1+t(i)+s(i)-(2*r(i)))-...
Y(e7(n))*(1+s(i))*(1-t(i))*(1+t(i)-s(i)-(2*r(i)))-...
Y(e8(n))*(1+s(i))*(1+t(i))*(1-t(i)-s(i)-(2*r(i))))+...
1/4*(-Y(e9(n))*(l-s(i))*(1-t(i)^2)-Y(e10(n))*(1-s(i)^2)*(l-t(i))-...
Y(e11(n))*(1+s(i))*(1-t(i)^2)-Y(e12(n))*(1-s(i)^2)*(1+t(i))+...
Y(e13(n))*(1-s(i))*(1-t(i)^2)+Y(e14(n))*(1-s(i)^2)*(1-t(i))+...
Y(e15(n))*(l+s(i))*(1-t(i)^2)+Y(e16(n))*(1-s(i)^2)*(l+t(i))-...
Y(e17(n))*(2*r(i))*(1+s(i))*(1+t(i))-...
Y(el8(n))*(2*r(i))*(1-s(i))*(l+t(i))-...
Y(e19(n))*(2*r(i))*(1-s(i))*(1-t(i))-...
Y(e20(n))*(2*r(i))*(1+s(i))*(1-t(i))));
Ys(i)=(1/8*(Y(e1(n))*(1-r(i))*(1+t(i))*(1-t(i)+r(i)+(2*s(i)))+...
Y(e2(n))*(1-r(i))*(1-t(i))*(1+t(i)+r(i)+(2*s(i)))-...
Y(e3(n))*(1-r(i))*(1-t(i))*(1+t(i)+r(i)-(2*s(i)))-...
Y(e4(n))*(1-r(i))*(1+t(i))*(1-t(i)+r(i)-(2*s(i)))+...
Y(e5(n))*(1+r(i))*(1+t(i))*(1-t(i)-r(i)+(2*s(i)))+...
Y(e6(n))*(l+r(i))*(1-t(i))*(1+t(i)-r(i)+(2*s(i)))-...
Y(e7(n))*(1+r(i))*(1-t(i))*(1+t(i)-r(i)-(2*s(i)))-...
Y(e8(n))*(1+r(i))*(1+t(i))*(1-t(i)-r(i)-(2*s(i))))+...
1/4*(-Y(e9(n))*(l-r(i))*(1-t(i)^2)-...
Y(el0(n))*(2*s(i))*(1-r(i))*(l-t(i))+...
Y(e11(n))*(1-r(i))*(1-t(i)^2)-...
Y(e12(n))*(2*s(i))*(1-r(i))*(1+t(i))-...
Y(e13(n))*(1+r(i))*(1-t(i)^2)-...
Y(e14(n))*(2*s(i))*(1+r(i))*(1-t(i))+...
Y(e15(n))*(l+r(i))*(1-t(i)^2)-...
Y(e16(n))*(2*s(i))*(1+r(i))*(1+t(i))+...
Y(e17(n))*(1-r(i)^2)*(1+t(i))-...
Y(e18(n))*(1-r(i)^2)*(1+t(i))-...
Y(e19(n))*(1-r(i)^2)*(1-t(i))+...
Y(e20(n))*(l-r(i)^2)*(1-t(i))));
Yt(i)=(1/8*(-Y(e1(n))*(1-r(i))*(1-s(i))*(1+s(i)+r(i)-(2*t(i)))+...
Y(e2(n))*(1-r(i))*(1-s(i))*(1+s(i)+r(i)+(2*t(i)))+...
Y(e3(n))*(1-r(i))*(1+s(i))*(1-s(i)+r(i)+(2*t(i)))-...
Y(e4(n))*(1-r(i))*(1+s(i))*(1-s(i)+r(i)-(2*t(i)))-...
Y(e5(n))*(1+r(i))*(1-s(i))*(1+s(i)-r(i)-(2*t(i)))+...
Y(e6(n))*(1+r(i))*(1-s(i))*(1+s(i)-r(i)+(2*t(i)))+...
Y(e7(n))*(1+r(i))*(1+s(i))*(1-s(i)-r(i)+(2*t(i)))-...
Y(e8(n))*(l+r(i))*(1+s(i))*(1-s(i)-r(i)-(2*t(i))))+...
1/4*(-Y(e9(n))*(2*t(i))*(1-r(i))*(1-s(i))-...
Y(e10(n))*(1-r(i))*(1-s(i)^2)-...
Y(e11(n))*(2*t(i))*(1-r(i))*(1+s(i))+...
Y(e12(n))*(1-r(i))*(1-s(i)^2)-...
Y(e13(n))*(2*t(i))*(1+r(i))*(1-s(i))-...
Y(e14(n))*(1+r(i))*(1-s(i)^2)-...
Y(e15(n))*(2*t(i))*(1+r(i))*(1+s(i))+...
Y(e16(n))*(l+r(i))*(1-s(i)^2)+...
Y(e17(n))*(1-r(i)^2)*(1+s(i))+...
Y(e18(n))*(1-r(i)^2)*(1-s(i))-...
Y(e19(n))*(1-r(i)^2)*(1-s(i))-...
Y(e20(n))*(1-r(i)^2)*(1+s(i))));
Zr(i)=(1/8*(Z(e1(n))*(1-s(i))*(1+t(i))*(1-t(i)+s(i)+(2*r(i)))+...
Z(e2(n))*(1-s(i))*(1-t(i))*(1+t(i)+s(i)+(2*r(i)))+...
Z(e3(n))*(1+s(i))*(1-t(i))*(1+t(i)-s(i)+(2*r(i)))+...
Z(e4(n))*(l+s(i))*(1+t(i))*(1-t(i)-s(i)+(2*r(i)))-...
Z(e5(n))*(1-s(i))*(1+t(i))*(1-t(i)+s(i)-(2*r(i)))-...
Z(e6(n))*(1-s(i))*(1-t(i))*(1+t(i)+s(i)-(2*r(i)))-...
Z(e7(n))*(1+s(i))*(1-t(i))*(1+t(i)-s(i)-(2*r(i)))-...
Z(e8(n))*(1+s(i))*(1+t(i))*(1-t(i)-s(i)-(2*r(i))))+...
1/4*(-Z(e9(n))*(1-s(i))*(1-t(i)^2)-Z(e10(n))*(1-s(i)^2)*(1-t(i))-...
Z(e1l(n))*(1+s(i))*(1-t(i)^2)-Z(e12(n))*(1-s(i)^2)*(1+t(i))+...
Z(e13(n))*(1-s(i))*(1-t(i)^2)+Z(e14(n))*(1-s(i)^2)*(1*t(i))+...
Z(el5(n))*(1+s(i))*(1-t(i)^2)+Z(e16(n))*(l-s(i)^2)*(1+t(i))-...
Z(e17(n))*(2*r(i))*(1+s(i))*(1+t(i))-...
Z(e18(n))*(2*r(i))*(1-s(i))*(1+t(i))-...
z(e19(n))*(2*r(i))*(1-s(i))*(1-t(i))-...
Z(e20(n))*(2*r(i))*(1+s(i))*(1-t(i))));
Zs(i)=(1/8*(Z(e1(n))*(1-r(i))*(l+t(i))*(1-t(i)+r(i)+(2*s(i)))+...
Z(e2(n))*(l-r(i))*(1-t(i))*(1+t(i)+r(i)+(2*s(i)))-...
Z(e3(n))*(l-r(i))*(1-t(i))*(1+t(i)+r(i)-(2*s(i)))-...
z(e4(n))*(1-r(i))*(1+t(i))*(1-t(i)+r(i)-(2*s(i)))+...
Z(e5(n))*(1+r(i))*(l+t(i))*(l-t(i)-r(i)+(2*s(i)))+...
z(e6(n))*(1+r(i))*(1-t(i))*(1+t(i)-r(i)+(2*s(i)))-...
Z(e7(n))*(1+r(i))*(1-t(i))*(1+t(i)-r(i)-(2*s(i)))-...
Z(e8(n))*(1+r(i))*(l+t(i))*(1-t(i)-r(i)-(2*s(i))))+...
1/4*(-Z(e9(n))*(1-r(i))*(1-t(i)^2)-...
Z(e10(n))*(2*s(i))*(1-r(i))*(1-t(i))+...
Z(e11(n))*(l-r(i))*(1-t(i)^2)-...
Z(e12(n))*(2*s(i))*(1-r(i))*(1+t(i))-...
Z(e13(n))*(1+r(i))*(1-t(i)^2)-...
Z(e14(n))*(2*s(i))*(1+r(i))*(1-t(i))+...
Z(e15(n))*(1+r(i))*(1-t(i)^2)-...
Z(e16(n))*(2*s(i))*(1+r(i))*(1+t(i))+...
Z(e17(n))*(1-r(i)^2)*(1+t(i))-...
Z(el8(n))*(1-r(i)^2)*(1+t(i))-...
Z(e19(n))*(1-r(i)^2)*(1-t(i))+...
Z(e20(n))*(1-r(i)^2)*(1-t(i))));
Zt(i)=(1/8*(-z(e1(n))*(1-r(i))*(1-s(i))*(1+s(i)+r(i)-(2*t(i)))+...
Z(e2(n))*(1-r(i))*(1-s(i))*(1+s(i)+r(i)+(2*t(i)))+...
Z(e3(n))*(1-r(i))*(1+s(i))*(1-s(i)+r(i)+(2*t(i)))-...
Z(e4(n))*(1-r(i))*(1+s(i))*(1-s(i)+r(i)-(2*t(i)))-...
Z(e5(n))*(1+r(i))*(1-s(i))*(1+s(i)-r(i)-(2*t(i)))+...
Z(e6(n))*(1+r(i))*(1-s(i))*(1+s(i)-r(i)+(2*t(i)))+...
Z(e7(n))*(1+r(i))*(1+s(i))*(l-s(i)-r(i)+(2*t(i)))-...
Z(e8(n))*(1+r(i))*(1+s(i))*(1-s(i)-r(i)-(2*t(i))))+...
1/4*(-Z(e9(n))*(2*t(i))*(1-r(i))*(1-s(i))-...
Z(el0(n))*(1-r(i))*(1-s(i)^2)-...
Z(e11(n))*(2*t(i))*(1-r(i))*(1+s(i))+...
Z(e12(n))*(1-r(i))*(1-s(i)^2)-...
Z(e13(n))*(2*t(i))*(1+r(i))*(1-s(i))-...
Z(el4(n))*(1+r(i))*(1-s(i)^2)-...
Z(el5(n))*(2*t(i))*(1+r(i))*(1+S(i))+...
Z(e16(n))*(1+r(i))*(1-s(i)^2)+...
Z(e17(n))*(1-r(i)^2)*(1+s(i))+...
Z(e18(n))*(1-r(i)^2)*(1-s(i))-...
Z(e19(n))*(1-r(i)^2)*(l-s(i))-...
Z(e20(n))*(1-r(i)^2)*(1+s(i))));
N1r(i)=1/8*(1-s(i))*(1+t(i))*(1-t(i)+s(i)+2*r(i));
N2r(i)=1/8*(1-s(i))*(1-t(i))*(1+t(i)+s(i)+2*r(i));
N3r(i)=1/8*(1+s(i))*(1-t(i))*(1+t(i)-s(i)+2*r(i));
N4r(i)=1/8*(1+s(i))*(1+t(i))*(1-t(i)-s(i)+2*r(i));
N5r(i)=-1/8*(1-s(i))*(1+t(i))*(1-t(i)+s(i)-2*r(i));
N6r(i)=-1/8*(1-s(i))*(1-t(i))*(1+t(i)+s(i)-2*r(i));
N7r(i)=-1/8*(1+s(i))*(1-t(i))*(1+t(i)-s(i)-2*r(i));
N8r(i)=-1/8*(1+s(i))*(1+t(i))*(1-t(i)-s(i)-2*r(i));
N9r(i)=-l/4*(l-s(i))*(1-t(i)^2);
N10r(i)=-1/4*(1-s(i)^2)*(1-t(i));
Nl1r(i)=-1/4*(1+s(i))*(1-t(i)^2);
N12r(i)=-1/4*(1-s(i)^2)*(1+t(i));
N13r(i)=1/4*(l-s(i))*(l-t(i)^2);
N14r(i)=1/4*(1-s(i)^2)*(l-t(i));
Nl5r(i)=1/4*(l+s(i))*(1-t(i)^2);
N16r(i)=1/4*(1-s(i)^2)*(1+t(i));
Nl7r(i)=-l/4*(2*r(i))*(l+s(i))*(1+t(i));
Nl8r(i)=-1/4*(2*r(i))*(l-s(i))*(1+t(i));
N19r(i)=-1/4*(2*r(i))*(1-s(i))*(1-t(i));
N20r(i)=-1/4*(2*r(i))*(l+s(i))*(1-t(i));


N1s(i)=1/8*(l-r(i))*(l+t(i))*(1-t(i)+r(i)+2*s(i));
N2s(i)=1/8*(1-r(i))*(1-t(i))*(1+t(i)+r(i)+2*s(i));
N3s(i)=-1/8*(l-r(i))*(1-t(i))*(1+t(i)+r(i)-2*s(i));
N4s(i)=-1/8*(l-r(i))*(1+t(i))*(1-t(i)+r(i)-2*s(i));
N5s(i)=1/8*(1+r(i))*(1+t(i))*(1-t(i)-r(i)+2*s(i));
N6s(i)=1/8*(l+r(i))*(1-t(i))*(1+t(i)-r(i)+2*s(i));
N7s(i)=-1/8*(1+r(i))*(l-t(i))*(1+t(i)-r(i)-2*s(i));
N8s(i)=-1/8*(l+r(i))*(1+t(i))*(1-t(i)-r(i)-2*s(i));
N9s(i)=-1/4*(1-r(i))*(1-t(i)^2);
N10s(i)=-1/4*((2*s(i))*(l-r(i))*(1-t(i)));
Nl1s(i)=1/4*(l-r(i))*(1-t(i)^2);
N12s(i)=-1/4*(2*s(i))*(1-r(i))*(1+t(i));
Nl3s(i)=-1/4*(l+r(i))*(1-t(i)^2);
N14s(i)=-1/4*(2*s(i))*(1+r(i))*(1-t(i));
N15s(i)=l/4*(1+r(i))*(1-t(i)^2);
Nl6s(i)=-1/4*(2*s(i))*(1+r(i))*(1+t(i));
N17s(i)=1/4*(1-r(i)^2)*(1+t(i));
N18s(i)=-1/4*(1-r(i)^2)*(l+t(i));
N19s(i)=-1/4*(1-r(i)^2)*(1-t(i));
N20s(i)=1/4*(1-r(i)^2)*(l-t(i));


Nlt(i)=-1/8*(l-r(i))*(l-s(i))*(1+s(i)+r(i)-2*t(i));
N2t(i)=1/8*(1-r(i))*(1-s(i))*(1+s(i)+r(i)+2*t(i));
N3t(i)=1/8*(l-r(i))*(1+s(i))*(1-s(i)+r(i)+2*t(i));
N4t(i)=-1/8*(1-r(i))*(1+s(i))*(1-s(i)+r(i)-2*t(i));
N5t(i)=-1/8*(1+r(i))*(1-s(i))*(1+s(i)-r(i)-2*t(i));
N6t(i)=1/8*(l+r(i))*(1-s(i))*(1+s(i)-r(i)+2*t(i));
N7t(i)=1/8*(1+r(i))*(1+s(i))*(1-s(i)-r(i)+2*t(i));
N8t(i)=-1/8*(1+r(i))*(1+s(i))*(1-s(i)-r(i)-2*t(i));
N9t(i)=-1/4*(2*t(i))*(1-r(i))*(1-s(i));
N10t(i)=-1/4*(1-r(i))*(1-s(i)^2);
Nl1t(i)=-1/4*(2*t(i))*(1-r(i))*(1+s(i));
N12t(i)=1/4*(1-r(i))*(1-s(i)^2);
Nl3t(i)=-l/4*(2*t(i))*(1+r(i))*(1-s(i));
N14t(i)=-1/4*(l+r(i))*(l-s(i)^2);
N15t(i)=-1/4*(2*t(i))*(1+r(i))*(1+s(i));
N16t(i)=1/4*(l+r(i))*(1-s(i)^2);
N17t(i)=1/4*(1-r(i)^2)*(1+s(i));
Nl8t(i)=1/4*(1-r(i)^2)*(l-s(i));
N19t(i)=-1/4*(1-r(i)^2)*(1-s(i));
N20t(i)=-1/4*(1-r(i)^2)*(1+s(i));

J(i)=Xr(i)*(Ys(i)*Zt(i)-Zs(i)*Yt(i))-Yr(i)*(Xs(i)*Zt(i)-Zs(i)*Xt(i))+...
 Zr(i)*(Xs(i)*Yt(i)-Ys(i)*Xt(i));
b11=N1r(i)*(Ys(i)*Zt(i)-Zs(i)*Yt(i))-Yr(i)*(Nls(i)*Zt(i)-Zs(i)*N1t(i))+...
 Zr(i)*(Nls(i)*Yt(i)-Ys(i)*Nlt(i));
b12=Xr(i)*(N1s(i)*Zt(i)-Zs(i)*N1t(i))-Nlr(i)*(Xs(i)*Zt(i)-Zs(i)*Xt(i))+...
 Zr(i)*(Xs(i)*N1t(i)-N1s(i)*Xt(i));
bl3=Xr(i)*(Ys(i)*N1t(i)-Nls(i)*Yt(i))-Yr(i)*(Xs(i)*N1t(i)-N1s(i)*Xt(i))+...
N1r(i)*(Xs(i)*Yt(i)-Ys(i)*Xt(i));
b21=N2r(i)*(Ys(i)*Zt(i)-Zs(i)*Yt(i))-Yr(i)*(N2s(i)*Zt(i)-Zs(i)*N2t(i))+...
Zr(i)*(N2s(i)*Yt(i)-Ys(i)*N2t(i));
b22=Xr(i)*(N2s(i)*Zt(i)-Zs(i)*N2t(i))-N2r(i)*(Xs(i)*Zt(i)-Zs(i)*Xt(i))+...
Zr(i)*(XS(i)*N2t(i)-N2S(i)*Xt(i));
b23=Xr(i)*(Ys(i)*N2t(i)-N2s(i)*Yt(i))-Yr(i)*(Xs(i)*N2t(i)-N2s(i)*Xt(i))+...
N2r(i)*(Xs(i)*Yt(i)-Ys(i)*Xt(i));
b31=N3r(i)*(Ys(i)*Zt(i)-Zs(i)*Yt(i))-Yr(i)*(N3s(i)*Zt(i)-Zs(i)*N3t(i))+...
Zr(i)*(N3s(i)*Yt(i)-Ys(i)*N3t(i));
b32=Xr(i)*(N3s(i)*Zt(i)-Zs(i)*N3t(i))-N3r(i)*(Xs(i)*Zt(i)-Zs(i)*Xt(i))+...
Zr(i)*(Xs(i)*N3t(i)-N3s(i)*Xt(i));
b33=Xr(i)*(Ys(i)*N3t(i)-N3s(i)*Yt(i))-Yr(i)*(Xs(i)*N3t(i)-N3s(i)*Xt(i))+...
N3r(i)*(Xs(i)*Yt(i)-Ys(i)*Xt(i));
b41=N4r(i)*(Ys(i)*Zt(i)-Zs(i)*Yt(i))-Yr(i)*(N4s(i)*Zt(i)-Zs(i)*N4t(i))+...
Zr(i)*(N4s(i)*Yt(i)-Ys(i)*N4t(i));
b42=Xr(i)*(N4s(i)*Zt(i)-Zs(i)*N4t(i))-N4r(i)*(Xs(i)*Zt(i)-Zs(i)*Xt(i))+...
Zr(i)*(Xs(i)*N4t(i)-N4s(i)*Xt(i));
b43=Xr(i)*(Ys(i)*N4t(i)-N4s(i)*Yt(i))-Yr(i)*(Xs(i)*N4t(i)-N4s(i)*Xt(i))+...
N4r(i)*(Xs(i)*Yt(i)-Ys(i)*Xt(i));
b5l=N5r(i)*(Ys(i)*Zt(i)-Zs(i)*Yt(i))-Yr(i)*(N5s(i)*Zt(i)-Zs(i)*N5t(i))+...
Zr(i)*(N55(i)*Yt(i)-Ys(i)*N5t(i));
b52=Xr(i)*(N5s(i)*Zt(i)-Zs(i)*N5t(i))-N5r(i)*(XS(i)*Zt(i)-Zs(i)*Xt(i))+...
Zr(i)*(Xs(i)*N5t(i)-N5s(i)*Xt(i));
b53=Xr(i)*(Ys(i)*N5t(i)-N5s(i)*Yt(i))-Yr(i)*(Xs(i)*N5t(i)-N5s(i)*Xt(i))+...
N5r(i)*(Xs(i)*Yt(i)-Ys(i)*Xt(i));
b6l=N6r(i)*(Ys(i)*Zt(i)-Zs(i)*Yt(i))-Yr(i)*(N6s(i)*Zt(i)-Zs(i)*N6t(i))+...
Zr(i)*(N65(i)*Yt(i)-Ys(i)*N6t(i));
b62=Xr(i)*(N6s(i)*Zt(i)-Zs(i)*N6t(i))-N6r(i)*(Xs(i)*Zt(i)-Zs(i)*Xt(i))+...
Zr(i)*(Xs(i)*N6t(i)-N6s(i)*Xt(i));
b63=Xr(i)*(Ys(i)*N6t(i)-N6s(i)*Yt(i))-Yr(i)*(Xs(i)*N6t(i)-N6s(i)*Xt(i))+...
N6r(i)*(Xs(i)*Yt(i)-Ys(i)*Xt(i));
b71=N7r(i)*(Ys(i)*Zt(i)-Zs(i)*Yt(i))-Yr(i)*(N7s(i)*Zt(i)-Zs(i)*N7t(i))+...
Zr(i)*(N7s(i)*Yt(i)-Y5(i)*N7t(i));
b72=Xr(i)*(N7s(i)*Zt(i)-Zs(i)*N7t(i))-N7r(i)*(Xs(i)*Zt(i)-Zs(i)*Xt(i))+...
Zr(i)*(Xs(i)*N7t(i)-N7s(i)*Xt(i));
b73=Xr(i)*(Ys(i)*N7t(i)-N7s(i)*Yt(i))-Yr(i)*(Xs(i)*N7t(i)-N7s(i)*Xt(i))+...
N7r(i)*(Xs(i)*Yt(i)-Ys(i)*Xt(i));
b81=N8r(i)*(Ys(i)*Zt(i)-Zs(i)*Yt(i))-Yr(i)*(N8s(i)*Zt(i)-Zs(i)*N8t(i))+...
Zr(i)*(N8s(i)*Yt(i)-Ys(i)*N8t(i));
b82=Xr(i)*(N8s(i)*Zt(i)-Zs(i)*N8t(i))-N8r(i)*(Xs(i)*Zt(i)-Zs(i)*Xt(i))+...
Zr(i)*(Xs(i)*N8t(i)-N8s(i)*Xt(i));
b83=Xr(i)*(Ys(i)*N8t(i)-N8s(i)*Yt(i))-Yr(i)*(Xs(i)*N8t(i)-N8s(i)*Xt(i))+...
N8r(i)*(Xs(i)*Yt(i)-Ys(i)*Xt(i));
b91=N9r(i)*(Ys(i)*Zt(i)-Zs(i)*Yt(i))-Yr(i)*(N9s(i)*Zt(i)-Zs(i)*N9t(i))+...
Zr(i)*(N9S(i)*Yt(i)-Ys(i)*N9t(i));
b92=Xr(i)*(N9s(i)*Zt(i)-Zs(i)*N9t(i))-N9r(i)*(Xs(i)*Zt(i)-Zs(i)*Xt(i))+...
Zr(i)*(Xs(i)*N9t(i)-N9s(i)*Xt(i));
b93=Xr(i)*(Ys(i)*N9t(i)-N9S(i)*Yt(i))-Yr(i)*(Xs(i)*N9t(i)-N9s(i)*Xt(i))+...
N9r(i)*(Xs(i)*Yt(i)-Ys(i)*Xt(i));
b101=Nl0r(i)*(Ys(i)*Zt(i)-Zs(i)*Yt(i))-Yr(i)*(N10s(i)*Zt(i)-Zs(i)*N10t(i))+...
Zr(i)*(N10s(i)*Yt(i)-Ys(i)*N10t(i));
b102=Xr(i)*(N10s(i)*Zt(i)-Zs(i)*N10t(i))-Nl0r(i)*(Xs(i)*Zt(i)-Zs(i)*Xt(i))+...
Zr(i)*(Xs(i)*N10t(i)-Nl0s(i)*Xt(i));
b103=Xr(i)*(Ys(i)*N10t(i)-N10s(i)*Yt(i))-Yr(i)*(Xs(i)*N10t(i)-N10s(i)*Xt(i))+...
N10r(i)*(Xs(i)*Yt(i)-Ys(i)*Xt(i));
b1l1=N11r(i)*(Ys(i)*Zt(i)-Zs(i)*Yt(i))-Yr(i)*(N11s(i)*Zt(i)-Zs(i)*N11t(i))+...
Zr(i)*(Nlls(i)*Yt(i)-Ys(i)*N11t(i));
b112=Xr(i)*(N11s(i)*Zt(i)-Zs(i)*N11t(i))-N11r(i)*(Xs(i)*Zt(i)-Zs(i)*Xt(i))+...
Zr(i)*(Xs(i)*N11t(i)-N1ls(i)*Xt(i));
b113=Xr(i)*(Ys(i)*N11t(i)-N11s(i)*Yt(i))-Yr(i)*(Xs(i)*N11t(i)-N11s(i)*Xt(i))+...
N11r(i)*(Xs(i)*Yt(i)-Ys(i)*Xt(i));
bl2l=N12r(i)*(Ys(i)*Zt(i)-Zs(i)*Yt(i))-Yr(i)*(Nl2s(i)*Zt(i)-Zs(i)*N12t(i))+...
Zr(i)*(Nl2s(i)*Yt(i)-Ys(i)*N12t(i));
b122=Xr(i)*(Nl2s(i)*Zt(i)-Zs(i)*N12t(i))-N12r(i)*(Xs(i)*Zt(i)-Zs(i)*Xt(i))+...
Zr(i)*(Xs(i)*N12t(i)-Nl2s(i)*Xt(i));
bl23=Xr(i)*(Ys(i)*N12t(i)-N12s(i)*Yt(i))-Yr(i)*(Xs(i)*N12t(i)-Nl2s(i)*Xt(i))+...
Nl2r(i)*(Xs(i)*Yt(i)-Ys(i)*Xt(i));
b131=N13r(i)*(Ys(i)*Zt(i)-Zs(i)*Yt(i))-Yr(i)*(N13s(i)*Zt(i)-Zs(i)*N13t(i))+...
Zr(i)*(N13s(i)*Yt(i)-Ys(i)*N13t(i));
bl32=Xr(i)*(Nl3s(i)*Zt(i)-Zs(i)*Nl3t(i))-Nl3r(i)*(Xs(i)*Zt(i)-Zs(i)*Xt(i))+...
Zr(i)*(Xs(i)*N13t(i)-N13s(i)*Xt(i));
bl33=Xr(i)*(Ys(i)*N13t(i)-N13s(i)*Yt(i))-Yr(i)*(Xs(i)*N13t(i)-Nl3s(i)*Xt(i))+...
N13r(i)*(Xs(i)*Yt(i)-Ys(i)*Xt(i));
bl41=N14r(i)*(Ys(i)*Zt(i)-Zs(i)*Yt(i))-Yr(i)*(Nl4s(i)*Zt(i)-Zs(i)*N14t(i))+...
Zr(i)*(Nl4s(i)*Yt(i)-Ys(i)*N14t(i));
bl42=Xr(i)*(Nl4s(i)*Zt(i)-Zs(i)*Nl4t(i))-N14r(i)*(Xs(i)*Zt(i)-Zs(i)*Xt(i))+...
Zr(i)*(Xs(i)*N14t(i)-N14s(i)*Xt(i));
bl43=Xr(i)*(Ys(i)*N14t(i)-N14s(i)*Yt(i))-Yr(i)*(Xs(i)*N14t(i)-Nl4s(i)*Xt(i))+...
Nl4r(i)*(Xs(i)*Yt(i)-Ys(i)*Xt(i));
bl51=N15r(i)*(Ys(i)*Zt(i)-Zs(i)*Yt(i))-Yr(i)*(N15s(i)*Zt(i)-Zs(i)*N15t(i))+...
Zr(i)*(N15s(i)*Yt(i)-Ys(i)*N15t(i));
bl52=Xr(i)*(N15s(i)*Zt(i)-Zs(i)*Nl5t(i))-N15r(i)*(Xs(i)*Zt(i)-Zs(i)*Xt(i))+...
Zr(i)*(Xs(i)*N15t(i)-Nl5s(i)*Xt(i));
b153=Xr(i)*(Ys(i)*N15t(i)*N15s(i)*Yt(i))-Yr(i)*(Xs(i)*N15t(i)-N15s(i)*Xt(i))+...
Nl5r(i)*(Xs(i)*Yt(i)-Ys(i)*Xt(i));
bl61=N16r(i)*(Ys(i)*Zt(i)-Zs(i)*Yt(i))-Yr(i)*(Nl6s(i)*Zt(i)-Zs(i)*N16t(i))+...
Zr(i)*(Nl6s(i)*Yt(i)-Ys(i)*N16t(i));
b162=Xr(i)*(N16s(i)*Zt(i)-Zs(i)*N16t(i))-N16r(i)*(Xs(i)*Zt(i)-Zs(i)*Xt(i))+...
Zr(i)*(Xs(i)*N16t(i)-Nl6s(i)*Xt(i));
b163=Xr(i)*(Ys(i)*N16t(i)-N16s(i)*Yt(i))-Yr(i)*(Xs(i)*N16t(i)-N16s(i)*Xt(i))+...
N16r(i)*(Xs(i)*Yt(i)-Ys(i)*Xt(i));
b171=N17r(i)*(Ys(i)*Zt(i)-Zs(i)*Yt(i))-Yr(i)*(Nl7s(i)*Zt(i)-Zs(i)*Nl7t(i))+...
Zr(i)*(Nl7s(i)*Yt(i)-Ys(i)*Nl7t(i));
b172=Xr(i)*(N17s(i)*Zt(i)-Zs(i)*N17t(i))-N17r(i)*(Xs(i)*Zt(i)-Zs(i)*Xt(i))+...
Zr(i)*(Xs(i)*N17t(i)-N17s(i)*Xt(i));
bl73=Xr(i)*(Ys(i)*N17t(i)-N17s(i)*Yt(i))-Yr(i)*(Xs(i)*N17t(i)-Nl7S(i)*Xt(i))+...
N17r(i)*(Xs(i)*Yt(i)-Ys(i)*Xt(i));
b181=N18r(i)*(Ys(i)*Zt(i)-Zs(i)*Yt(i))-Yr(i)*(N18s(i)*Zt(i)-Zs(i)*N18t(i))+...
Zr(i)*(N18s(i)*Yt(i)-Ys(i)*N18t(i));
bl82=Xr(i)*(Nl8s(i)*Zt(i)-Zs(i)*N18t(i))-N18r(i)*(Xs(i)*Zt(i)-Zs(i)*Xt(i))+...
Zr(i)*(Xs(i)*N18t(i)-N18s(i)*Xt(i));
bl83=Xr(i)*(Ys(i)*Nl8t(i)-N18s(i)*Yt(i))-Yr(i)*(Xs(i)*N18t(i)-Nl8s(i)*Xt(i))+...
N18r(i)*(Xs(i)*Yt(i)-Ys(i)*Xt(i));
b191=N19r(i)*(Ys(i)*Zt(i)-Zs(i)*Yt(i))-Yr(i)*(N19s(i)*Zt(i)-Zs(i)*N19t(i))+...
Zr(i)*(N19s(i)*Yt(i)-Ys(i)*N19t(i));
bl92=Xr(i)*(Nl9s(i)*Zt(i)-Zs(i)*N19t(i))-N19r(i)*(Xs(i)*Zt(i)-Zs(i)*Xt(i))+...
Zr(i)*(Xs(i)*N19t(i)-N19s(i)*Xt(i));
bl93=Xr(i)*(Ys(i)*N19t(i)-N19S(i)*Yt(i))-Yr(i)*(Xs(i)*N19t(i)-Nl9s(i)*Xt(i))+...
N19r(i)*(Xs(i)*Yt(i)-Ys(i)*Xt(i));
bl91=N19r(i)*(Ys(i)*Zt(i)-Zs(i)*Yt(i))-Yr(i)*(Nl9s(i)*Zt(i)-Zs(i)*N19t(i))+...
Zr(i)*(N19s(i)*Yt(i)-Ys(i)*N19t(i));
bl92=Xr(i)*(Nl9s(i)*Zt(i)-Zs(i)*Nl9t(i))-N19r(i)*(Xs(i)*Zt(i)-Zs(i)*Xt(i))+...
Zr(i)*(Xs(i)*N19t(i)-N19s(i)*Xt(i));
b193=Xr(i)*(Ys(i)*N19t(i)-N19s(i)*Yt(i))-Yr(i)*(Xs(i)*N19t(i)-N19s(i)*Xt(i))+...
N19r(i)*(Xs(i)*Yt(i)-Ys(i)*Xt(i));
b201=N20r(i)*(Ys(i)*Zt(i)-Zs(i)*Yt(i))-Yr(i)*(N20s(i)*Zt(i)-Zs(i)*N20t(i))+...
Zr(i)*(N20s(i)*Yt(i)-Ys(i)*N20t(i));
b202=Xr(i)*(N20s(i)*Zt(i)-Zs(i)*N20t(i))-N20r(i)*(Xs(i)*Zt(i)-Zs(i)*Xt(i))+...
Zr(i)*(Xs(i)*N20t(i)-N20s(i)*Xt(i));
b203=Xr(i)*(Ys(i)*N20t(i)-N20s(i)*Yt(i))-Yr(i)*(Xs(i)*N20t(i)-N20s(i)*Xt(i))+...
N20r(i)*(Xs(i)*Yt(i)-Ys(i)*Xt(i));
B(:,:,i)=1/J(i)*[bll 0 0 b2l 0 0 b3l 0 0 b41 0 0 b51 0 0 b6l 0 0 b7l 0 0 ...
b81 0 0 b91 0 0 b101 0 0 blll 0 0 b121 0 0 bl3l 0 0 b141 ...
0 0 b151 0 0 b161 0 0 b171 0 0 b181 0 0 b191 0 0 b201 0 0 ...
0 b12 0 0 b22 0 0 b32 0 0 b42 0 0 b52 0 0 b62 0 0 b72 0 0 ...
b82 0 0 b92 0 0 b102 0 0 b112 0 0 b122 0 0 bl32 0 0 bl42 ...
0 0 b152 0 0 b162 0 0 b172 0 0 bl82 0 0 b192 0 0 b202 0 ...
0 0 b13 0 0 b23 0 0 b33 0 0 b43 0 0 b53 0 0 b63 0 0 b73 0 ...
0 b83 0 0 b93 0 0 bl03 0 0 b113 0 0 b123 0 0 b133 0 0 b143 ...
0 0 b153 0 0 b163 0 0 b173 0 0 b183 0 0 b193 0 0 b203 ...
b12 bll 0 b22 b2l 0 b32 b3l 0 b42 b41 0 b52 b51 0 b62 b6l 0 b72 b7l 0 b82 ...
b81 0 b92 b91 0 b102 bl01 0 b112 blll 0 b122 b121 0 bl32 bl3l 0 bl42 b141 ...
0 b152 b151 0 b162 b161 0 b172 b171 0 bl82 b18l 0 b192 b191 0 b202 b201 0 ...
0 b13 b12 0 b23 b22 0 b33 b32 0 b43 b42 0 b53 b52 0 b63 b62 0 b73 b72 0 b83 ...
b82 0 b93 b92 0 bl03 b102 0 b113 b112 0 b123 b122 0 b133 bl32 0 b143 bl42 ...
0 b153 b152 0 b163 b162 0 b173 b172 0 b183 bl82 0 b193 b192 0 b203 b202 ...
b13 0 bll b23 0 b2l b33 0 b3l b43 0 b41 b53 0 b51 b63 0 b6l b73 0 b7l b83 0 ...
b81 b93 0 b9l bl03 0 b10l bl13 0 blll b123 0 b12l b133 0 bl3l b143 0 b141 ...
b153 0 b15l b163 0 b161 b173 0 b171 b183 0 b181 b193 0 b191 b203 0 b201];
end
Bl(:,:,n)=B(:,:,1); B2(:,:,n)=B(:,:,2); B3(:,:,n)=B(:,:,3);
B4(:,:,n)=B(:,:,4); B5(:,:,n)=B(:,:,5); B6(:,:,n)=B(:,:,6);
B7(:,:,n)=B(:,:,7); B8(:,:,n)=B(:,:,8); B9(:,:,n)=B(:,:,9);
B10(:,:,n)=B(:,:,l0); Bl1(:,:,n)=B(:,:,l1); B12(:,:,n)=B(:,:,12);
B13(:,:,n)=B(:,:,l3); B14(:,:,n)=B(:,:,14);

Jl(n)=J(1); J2(n)=J(2); J3(n)=J(3); J4(n)=J(4); J5(n)=J(5);
J6(n)=J(6); J7(n)=J(7); J8(n)=J(8); J9(n)=J(9); J10(n)=J(l0);
J11(n)=J(11); J12(n)=J(12); Jl3(n)=J(13); Jl4(n)=J(14);

%Definition of stiffhess(loca1) matrix '[k]' for each element:
k(:,:,n)=(w1)*Bl(:,:,n)'*D*Bl(:,:,n)*Jl(n)+...
(w1)*B2(:,:,n)'*D*B2(:,:,n)*J2(n)+...
(w1)*B3(:,:,n)'*D*B3(:,:,n)*J3(n)+...
(w1)*B4(:,:,n)'*D*B4(:,:,n)*J4(n)+...
(w1)*B5(:,:,n)'*D*B5(:,:,n)*J5(n)+...
(w1)*B6(:,:,n)'*D*B6(:,:,n)*J6(n)+...
(w1)*B7(:,:,n)'*D*B7(:,:,n)*J7(n)+...
(w1)*B8(:,:,n)'*D*B8(:,:,n)*J8(n)+...
(w2)*B9(:,:,n)'*D*B9(:,:,n)*J9(n)+...
(w2)*B10(:,:,n)'*D*B10(:,:,n)*J10(n)+...
(w2)*Bll(:,:,n)'*D*Bll(:,:,n)*Jl1(n)+...
(w2)*B12(:,:,n)'*D*B12(:,:,n)*J12(n)+...
(w2)*Bl3(:,:,n)'*D*Bl3(:,:,n)*Jl3(n)+...
(w2)*B14(:,:,n)'*D*B14(:,:,n)*J14(n);
end
%Definition of structure stiffhess(g1abal) matrix(plant of submatrixas):
K=zeros(3*NN,3*NN);
for n=1:NE
K(3*e1(n)-2:3*e1(n),3*e1(n)-2:3*e1(n))=...
K(3*e1(n)-2:3*e1(n),3*e1(n)-2:3*el(n))+k(1:3,1:3,n);
K(3*e1(n)-2:3*e1(n),3*e2(n)-2:3*e2(n))=...
K(3*e1(n)-2:3*e1(n),3*e2(n)-2:3*e2(n))+k(l:3,4:6,n);
K(3*el(n)-2:3*el(n),3*e3(n)-2:3*e3(n))=...
K(3*e1(n)-2:3*e1(n),3*e3(n)-2:3*e3(n))+k(1:3,7:9,n);
K(3*el(n)-2:3*el(n),3*e4(n)-2:3*e4(n))=...
K(3*e1(n)-2:3*e1(n),3*e4(n)-2:3*e4(n))+k(1:3,10:12,n);
K(3*el(n)-2:3*el(n),3*e5(n)-2:3*e5(n))=...
K(3*e1(n)-2:3*e1(n),3*e5(n)-2:3*e5(n))+k(l:3,13:l5,n);
K(3*el(n)-2:3*el(n),3*e6(n)-2:3*e6(n))=...
K(3*e1(n)-2:3*e1(n),3*e6(n)-2:3*e6(n))+k(1:3,16:18,n);
K(3*el(n)-2:3*el(n),3*e7(n)-2:3*e7(n))=...
K(3*e1(n)-2:3*e1(n),3*e7(n)-2:3*e7(n))+k(1:3,19:21,n);
K(3*e1(n)-2:3*e1(n),3*e8(n)-2:3*e8(n))=...
K(3*e1(n)-2:3*e1(n),3*e8(n)-2:3*e8(n))+k(1:3,22:24,n);
K(3*e1(n)-2:3*e1(n),3*e9(n)-2:3*e9(n))=...
K(3*e1(n)-2:3*e1(n),3*e9(n)-2:3*e9(n))+k(1:3,25:27,n);
K(3*e1(n)-2:3*e1(n),3*e10(n)-2:3*e10(n))=...
K(3*e1(n)-2:3*e1(n),3*el0(n)-2:3*e10(n))+k(1:3,28:30,n);
K(3*el(n)-2:3*el(n),3*e11(n)-2:3*e11(n))=...
K(3*e1(n)-2:3*e1(n),3*e11(n)-2:3*e11(n))+k(1:3,31:33,n);
K(3*el(n)-2:3*e1(n),3*e12(n)-2:3*e12(n))=...
K(3*e1(n)-2:3*e1(n),3*e12(n)-2:3*e12(n))+k(1:3,34:36,n);
K(3*e1(n)-2:3*e1(n),3*e13(n)-2:3*e13(n))=...
K(3*e1(n)-2:3*e1(n),3*e13(n)-2:3*e13(n))+k(l:3,37:39,n);
K(3*e1(n)-2:3*e1(n),3*e14(n)-2:3*e14(n))=...
K(3*e1(n)-2:3*e1(n),3*e14(n)-2:3*el4(n))+k(1:3,40:42,n);
K(3*e1(n)-2:3*e1(n),3*e15(n)-2:3*e15(n))=...
K(3*e1(n)-2:3*e1(n),3*el5(n)-2:3*e15(n))+k(1:3,43:45,n);
K(3*el(n)-2:3*el(n),3*e16(n)-2:3*e16(n))=...
K(3*e1(n)-2:3*e1(n),3*e16(n)-2:3*e16(n))+k(1:3,46:48,n);
K(3*e1(n)-2:3*e1(n),3*e17(n)-2:3*e17(n))=...
K(3*e1(n)-2:3*e1(n),3*e17(n)-2:3*e17(n))+k(1:3,49:51,n);
K(3*e1(n)-2:3*e1(n),3*e18(n)-2:3*e18(n))=...
K(3*e1(n)-2:3*e1(n),3*e18(n)-2:3*e18(n))+k(1:3,52:54,n);
K(3*e1(n)-2:3*e1(n),3*e19(n)-2:3*e19(n))=...
K(3*e1(n)-2:3*e1(n),3*e19(n)-2:3*el9(n))+k(1:3,55:57,n);
K(3*el(n)-2:3*e1(n),3*e20(n)-2:3*e20(n))=...
K(3*e1(n)-2:3*e1(n),3*e20(n)-2:3*e20(n))+k(1:3,58:60,n);
K(3*e2(n)-2:3*e2(n),3*e1(n)-2:3*e1(n))=...
K(3*e2(n)-2:3*e2(n),3*el(n)-2:3*el(n))+k(4:6,1:3,n);
K(3*e2(n)-2:3*e2(n),3*e2(n)-2:3*e2(n))=...
K(3*e2(n)-2:3*e2(n),3*e2(n)-2:3*e2(n))+k(4:6,4:6,n);
K(3*e2(n)-2:3*e2(n),3*e3(n)-2:3*e3(n))=...
K(3*e2(n)-2:3*e2(n),3*e3(n)-2:3*e3(n))+k(4:6,7:9,n);
K(3*e2(n)-2:3*e2(n),3*e4(n)-2:3*e4(n))=...
K(3*e2(n)-2:3*e2(n),3*e4(n)-2:3*e4(n))+k(4:6,10:12,n);
K(3*e2(n)-2:3*e2(n),3*e5(n)-2:3*e5(n))=...
K(3*e2(n)-2:3*e2(n),3*e5(n)-2:3*e5(n))+k(4:6,13:15,n);
K(3*e2(n)-2:3*e2(n),3*e6(n)-2:3*e6(n))=...
K(3*e2(n)-2:3*e2(n),3*e6(n)-2:3*e6(n))+k(4:6,16:18,n);
K(3*e2(n)-2:3*e2(n),3*e7(n)-2:3*e7(n))=...
K(3*e2(n)-2:3*e2(n),3*e7(n)-2:3*e7(n))+k(4:6,19:21,n);
K(3*e2(n)-2:3*e2(n),3*e8(n)-2:3*e8(n))=...
K(3*e2(n)-2:3*e2(n),3*eB(n)-2:3*e8(n))+k(4:6,22:24,n);
K(3*e2(n)-2:3*e2(n),3*e9(n)-2:3*e9(n))=...
K(3*e2(n)-2:3*e2(n),3*e9(n)-2:3*e9(n))+k(4:6,25:27,n);
K(3*e2(n)-2:3*e2(n),3*e10(n)-2:3*e10(n))=...
K(3*e2(n)-2:3*e2(n),3*e10(n)-2:3*e10(n))+k(4:6,28:30,n);
K(3*e2(n)-2:3*e2(n),3*e11(n)-2:3*e11(n))=...
K(3*e2(n)-2:3*e2(n),3*e1l(n)-2:3*ell(n))+k(4:6,31:33,n);
K(3*e2(n)-2:3*e2(n),3*e12(n)-2:3*e12(n))=...
K(3*e2(n)-2:3*e2(n),3*el2(n)-2:3*e12(n))+k(4:6,34:36,n);
K(3*e2(n)-2:3*e2(n),3*e13(n)-2:3*e13(n))=...
K(3*e2(n)-2:3*e2(n),3*el3(n)-2:3*e13(n))+k(4:6,37:39,n);
K(3*e2(n)-2:3*e2(n),3*e14(n)-2:3*e14(n))=...
K(3*e2(n)-2:3*e2(n),3*e14(n)-2:3*e14(n))+k(4:6,40:42,n);
K(3*e2(n)-2:3*e2(n),3*e15(n)-2:3*e15(n))=...
K(3*e2(n)-2:3*e2(n),3*e15(n)-2:3*e15(n))+k(4:6,43:45,n);
K(3*e2(n)-2:3*e2(n),3*e16(n)-2:3*e16(n))=...
K(3*e2(n)-2:3*e2(n),3*e16(n)-2:3*el6(n))+k(4:6,46:48,n);
K(3*e2(n)-2:3*e2(n),3*e17(n)-2:3*e17(n))=...
K(3*e2(n)-2:3*e2(n),3*el7(n)-2:3*e17(n))+k(4:6,49:51,n);
K(3*e2(n)-2:3*e2(n),3*e18(n)-2:3*e18(n))=...
K(3*e2(n)-2:3*e2(n),3*e18(n)-2:3*el8(n))+k(4:6,52:54,n);
K(3*e2(n)-2:3*e2(n),3*e19(n)-2:3*e19(n))=...
K(3*e2(n)-2:3*e2(n),3*e19(n)-2:3*e19(n))+k(4:6,55:57,n);
K(3*e2(n)-2:3*e2(n),3*e20(n)-2:3*e20(n))=...
K(3*e2(n)-2:3*e2(n),3*e20(n)-2:3*e20(n))+k(4:6,58:60,n);
K(3*e3(n)-2:3*e3(n),3*e1(n)-2:3*e1(n))=...
K(3*e3(n)-2:3*e3(n),3*e1(n)-2:3*e1(n))+k(7:9,1:3,n);
K(3*e3(n)-2:3*e3(n),3*e2(n)-2:3*e2(n))=...
K(3*e3(n)-2:3*e3(n),3*e2(n)-2:3*e2(n))+k(7:9,4:6,n);
K(3*e3(n)-2:3*e3(n),3*e3(n)-2:3*e3(n))=...
K(3*e3(n)-2:3*e3(n),3*e3(n)-2:3*e3(n))+k(7:9,7:9,n);
K(3*e3(n)-2:3*e3(n),3*e4(n)-2:3*e4(n))=...
K(3*e3(n)-2:3*e3(n),3*e4(n)-2:3*e4(n))+k(7:9,10:12,n);
K(3*e3(n)-2:3*e3(n),3*e5(n)-2:3*e5(n))=...
K(3*e3(n)-2:3*e3(n),3*e5(n)-2:3*e5(n))+k(7:9,13:l5,n);
K(3*e3(n)-2:3*e3(n),3*e6(n)-2:3*e6(n))=...
K(3*e3(n)-2:3*e3(n),3*e6(n)-2:3*e6(n))+k(7:9,16:18,n);
K(3*e3(n)-2:3*e3(n),3*e7(n)-2:3*e7(n))=...
K(3*e3(n)-2:3*e3(n),3*e7(n)-2:3*e7(n))+k(7:9,19:21,n);
K(3*e3(n)-2:3*e3(n),3*e8(n)-2:3*e8(n))=...
K(3*e3(n)-2:3*e3(n),3*e8(n)-2:3*e8(n))+k(7:9,22:24,n);
K(3*e3(n)-2:3*e3(n),3*e9(n)-2:3*e9(n))=...
K(3*e3(n)-2:3*e3(n),3*e9(n)-2:3*e9(n))+k(7:9,25:27,n);
K(3*e3(n)-2:3*e3(n),3*e10(n)-2:3*e10(n))=...
K(3*e3(n)-2:3*e3(n),3*el0(n)-2:3*e10(n))+k(7:9,28:30,n);
K(3*e3(n)-2:3*e3(n),3*e11(n)-2:3*e11(n))=...
K(3*e3(n)-2:3*e3(n),3*e11(n)-2:3*e11(n))+k(7:9,31:33,n);
K(3*e3(n)-2:3*e3(n),3*e12(n)-2:3*e12(n))=...
K(3*e3(n)-2:3*e3(n),3*e12(n)-2:3*e12(n))+k(7:9,34:36,n);
K(3*e3(n)-2:3*e3(n),3*e13(n)-2:3*e13(n))=...
K(3*e3(n)-2:3*e3(n),3*e13(n)-2:3*e13(n))+k(7:9,37:39,n);
K(3*e3(n)-2:3*e3(n),3*e14(n)-2:3*e14(n))=...
K(3*e3(n)-2:3*e3(n),3*e14(n)-2:3*el4(n))+k(7:9,40:42,n);
K(3*e3(n)-2:3*e3(n),3*e15(n)-2:3*e15(n))=...
K(3*e3(n)-2:3*e3(n),3*e15(n)-2:3*e15(n))+k(7:9,43:45,n);
K(3*e3(n)-2:3*e3(n),3*e16(n)-2:3*e16(n))=...
K(3*e3(n)-2:3*e3(n),3*e16(n)-2:3*e16(n))+k(7:9,46:48,n);
K(3*e3(n)-2:3*e3(n),3*e17(n)-2:3*e17(n))=...
K(3*e3(n)-2:3*e3(n),3*el7(n)-2:3*e17(n))+k(7:9,49:51,n);
K(3*e3(n)-2:3*e3(n),3*e18(n)-2:3*e18(n))=...
K(3*e3(n)-2:3*e3(n),3*e18(n)-2:3*e18(n))+k(7:9,52:54,n);
K(3*e3(n)-2:3*e3(n),3*e19(n)-2:3*e19(n))=...
K(3*e3(n)-2:3*e3(n),3*e19(n)-2:3*el9(n))+k(7:9,55:57,n);
K(3*e3(n)-2:3*e3(n),3*e20(n)-2:3*e20(n))=...
K(3*e3(n)-2:3*e3(n),3*e20(n)-2:3*e20(n))+k(7:9,58:60,n);
K(3*e4(n)-2:3*e4(n),3*e1(n)-2:3*e1(n))=...
K(3*e4(n)-2:3*e4(n),3*e1(n)-2:3*el(n))+k(10:12,1:3,n);
K(3*e4(n)-2:3*e4(n),3*e2(n)-2:3*e2(n))=...
K(3*e4(n)-2:3*e4(n),3*e2(n)-2:3*e2(n))+k(10:12,4:6,n);
K(3*e4(n)-2:3*e4(n),3*e3(n)-2:3*e3(n))=...
K(3*e4(n)-2:3*e4(n),3*e3(n)-2:3*e3(n))+k(10:12,7:9,n);
K(3*e4(n)-2:3*e4(n),3*e4(n)-2:3*e4(n))=...
K(3*e4(n)-2:3*e4(n),3*e4(n)-2:3*e4(n))+k(10:12,10:12,n);
K(3*e4(n)-2:3*e4(n),3*e5(n)*2:3*e5(n))=...
K(3*e4(n)-2:3*e4(n),3*e5(n)-2:3*e5(n))+k(10:12,13:l5,n);
K(3*e4(n)-2:3*e4(n),3*e6(n)-2:3*e6(n))=...
K(3*e4(n)-2:3*e4(n),3*e6(n)-2:3*e6(n))+k(10:12,16:l8,n);
K(3*e4(n)-2:3*e4(n),3*e7(n)-2:3*e7(n))=...
K(3*e4(n)-2:3*e4(n),3*e7(n)-2:3*e7(n))+k(10:12,19:21,n);
K(3*e4(n)-2:3*e4(n),3*e8(n)-2:3*e8(n))=...
K(3*e4(n)-2:3*e4(n),3*e8(n)-2:3*e8(n))+k(10:12,22:24,n);
K(3*e4(n)-2:3*e4(n),3*e9(n)-2:3*e9(n))=...
K(3*e4(n)-2:3*e4(n),3*e9(n)-2:3*e9(n))+k(10:12,25:27,n);
K(3*e4(n)-2:3*e4(n),3*e10(n)-2:3*e10(n))=...
K(3*e4(n)-2:3*e4(n),3*e10(n)-2:3*el0(n))+k(10:12,28:30,n);
K(3*e4(n)-2:3*e4(n),3*e11(n)-2:3*e11(n))=...
K(3*e4(n)-2:3*e4(n),3*e11(n)-2:3*e11(n))+k(10:12,31:33,n);
K(3*e4(n)-2:3*e4(n),3*e12(n)-2:3*e12(n))=...
K(3*e4(n)-2:3*e4(n),3*e12(n)-2:3*e12(n))+k(10:12,34:36,n);
K(3*e4(n)-2:3*e4(n),3*e13(n)-2:3*e13(n))=...
K(3*e4(n)-2:3*e4(n),3*el3(n)-2:3*e13(n))+k(l0:12,37:39,n);
K(3*e4(n)-2:3*e4(n),3*e14(n)-2:3*e14(n))=...
K(3*e4(n)-2:3*e4(n),3*e14(n)-2:3*e14(n))+k(10:12,40:42,n);
K(3*e4(n)-2:3*e4(n),3*e15(n)-2:3*e15(n))=...
K(3*e4(n)-2:3*e4(n),3*e15(n)-2:3*e15(n))+k(10:12,43:45,n);
K(3*e4(n)-2:3*e4(n),3*e16(n)-2:3*e16(n))=...
K(3*e4(n)-2:3*e4(n),3*e16(n)-2:3*e16(n))+k(10:12,46:48,n);
K(3*e4(n)-2:3*e4(n),3*e17(n)-2:3*e17(n))=...
K(3*e4(n)-2:3*e4(n),3*e17(n)-2:3*e17(n))+k(10:12,49:51,n);
K(3*e4(n)-2:3*e4(n),3*e1B(n)-2:3*e18(n))=...
K(3*e4(n)-2:3*e4(n),3*e18(n)-2:3*e18(n))+k(10:12,52:54,n);
K(3*e4(n)-2:3*e4(n),3*e19(n)-2:3*e19(n))=...
K(3*e4(n)-2:3*e4(n),3*el9(n)-2:3*e19(n))+k(10:l2,55:57,n);
K(3*e4(n)-2:3*e4(n),3*e20(n)-2:3*e20(n))=...
K(3*e4(n)-2:3*e4(n),3*e20(n)-2:3*e20(n))+k(10:12,58:60,n);
K(3*e5(n)-2:3*e5(n),3*e1(n)-2:3*e1(n))=...
K(3*e5(n)-2:3*e5(n),3*e1(n)-2:3*el(n))+k(13:15,1:3,n);
K(3*e5(n)-2:3*e5(n),3*e2(n)-2:3*e2(n))=...
K(3*e5(n)-2:3*e5(n),3*e2(n)-2:3*e2(n))+k(l3:15,4:6,n);
K(3*e5(n)-2:3*e5(n),3*e3(n)-2:3*e3(n))=...
K(3*e5(n)-2:3*e5(n),3*e3(n)-2:3*e3(n))+k(13:15,7:9,n);
K(3*e5(n)-2:3*e5(n),3*e4(n)-2:3*e4(n))=...
K(3*e5(n)-2:3*e5(n),3*e4(n)-2:3*e4(n))+k(13:15,10:l2,n);
K(3*e5(n)-2:3*e5(n),3*e5(n)-2:3*e5(n))=...
K(3*e5(n)-2:3*e5(n),3*e5(n)-2:3*e5(n))+k(13:15,13:l5,n);
K(3*e5(n)-2:3*e5(n),3*e6(n)-2:3*e6(n))=...
K(3*e5(n)-2:3*e5(n),3*e6(n)-2:3*e6(n))+k(13:15,16:l8,n);
K(3*e5(n)-2:3*e5(n),3*e7(n)-2:3*e7(n))=...
K(3*e5(n)-2:3*e5(n),3*e7(n)-2:3*e7(n))+k(13:15,19:21,n);
K(3*e5(n)-2:3*e5(n),3*e8(n)-2:3*e8(n))=...
K(3*e5(n)-2:3*e5(n),3*eB(n)-2:3*e8(n))+k(13:15,22:24,n);
K(3*e5(n)-2:3*e5(n),3*e9(n)-2:3*e9(n))=...
K(3*e5(n)-2:3*e5(n),3*e9(n)-2:3*e9(n))+k(13:15,25:27,n);
K(3*e5(n)-2:3*e5(n),3*e10(n)-2:3*e10(n))=...
K(3*e5(n)-2:3*e5(n),3*e10(n)-2:3*e10(n))+k(13:15,28:30,n);
K(3*e5(n)-2:3*e5(n),3*e11(n)-2:3*e11(n))=...
K(3*e5(n)-2:3*e5(n),3*e11(n)-2:3*e1l(n))+k(13:15,31:33,n);
K(3*e5(n)-2:3*e5(n),3*e12(n)-2:3*e12(n))=...
K(3*e5(n)-2:3*e5(n),3*el2(n)-2:3*e12(n))+k(13:l5,34:36,n);
K(3*e5(n)-2:3*e5(n),3*e13(n)-2:3*e13(n))=...
K(3*e5(n)-2:3*e5(n),3*e13(n)-2:3*e13(n))+k(13:15,37:39,n);
K(3*e5(n)-2:3*e5(n),3*e14(n)-2:3*e14(n))=...
K(3*e5(n)-2:3*e5(n),3*e14(n)-2:3*e14(n))+k(13:15,40:42,n);
K(3*e5(n)-2:3*e5(n),3*e15(n)-2:3*e15(n))=...
K(3*e5(n)-2:3*e5(n),3*e15(n)-2:3*e15(n))+k(13:15,43:45,n);
K(3*e5(n)-2:3*e5(n),3*e16(n)-2:3*e16(n))=...
K(3*e5(n)-2:3*e5(n),3*e16(n)-2:3*e16(n))+k(13:15,46:48,n);
K(3*e5(n)-2:3*e5(n),3*e17(n)-2:3*e17(n))=...
K(3*e5(n)-2:3*e5(n),3*el7(n)-2:3*e17(n))+k(13:l5,49:51,n);
K(3*e5(n)-2:3*e5(n),3*e18(n)-2:3*e18(n))=...
K(3*e5(n)-2:3*e5(n),3*e18(n)-2:3*e18(n))+k(13:15,52:54,n);
K(3*e5(n)-2:3*e5(n),3*e19(n)-2:3*e19(n))=...
K(3*e5(n)-2:3*e5(n),3*e19(n)-2:3*e19(n))+k(13:15,55:57,n);
K(3*e5(n)-2:3*e5(n),3*e20(n)-2:3*e20(n))=...
K(3*e5(n)-2:3*e5(n),3*e20(n)-2:3*e20(n))+k(13:15,58:60,n);
K(3*e6(n)-2:3*e6(n),3*e1(n)-2:3*e1(n))=...
K(3*e6(n)-2:3*e6(n),3*e1(n)-2:3*e1(n))+k(16:18,l:3,n);
K(3*e6(n)-2:3*e6(n),3*e2(n)-2:3*e2(n))=...
K(3*e6(n)-2:3*e6(n),3*e2(n)-2:3*e2(n))+k(16:18,4:6,n);
K(3*e6(n)-2:3*e6(n),3*e3(n)-2:3*e3(n))=...
K(3*e6(n)-2:3*e6(n),3*e3(n)-2:3*e3(n))+k(16:18,7:9,n);
K(3*e6(n)-2:3*e6(n),3*e4(n)-2:3*e4(n))=...
K(3*e6(n)-2:3*e6(n),3*e4(n)-2:3*e4(n))+k(16:18,10:l2,n);
K(3*e6(n)-2:3*e6(n),3*e5(n)-2:3*e5(n))=...
K(3*e6(n)-2:3*e6(n),3*e5(n)-2:3*e5(n))+k(16:18,13:l5,n);
K(3*e6(n)-2:3*e6(n),3*e6(n)-2:3*e6(n))=...
K(3*e6(n)-2:3*e6(n),3*e6(n)-2:3*e6(n))+k(l6:18,l6:18,n);
K(3*e6(n)-2:3*e6(n),3*e7(n)-2:3*e7(n))=...
K(3*e6(n)-2:3*e6(n),3*e7(n)-2:3*e7(n))+k(16:18,19:21,n);
K(3*e6(n)-2:3*e6(n),3*e8(n)-2:3*e8(n))=...
K(3*e6(n)-2:3*e6(n),3*e8(n)-2:3*e8(n))+k(16:18,22:24,n);
K(3*e6(n)-2:3*e6(n),3*e9(n)*2:3*e9(n))=...
K(3*e6(n)-2:3*e6(n),3*e9(n)-2:3*e9(n))+k(16:18,25:27,n);
K(3*e6(n)-2:3*e6(n),3*e10(n)-2:3*e10(n))=...
K(3*e6(n)-2:3*e6(n),3*e10(n)-2:3*e10(n))+k(16:18,28:30,n);
K(3*e6(n)-2:3*e6(n),3*e11(n)-2:3*e11(n))=...
K(3*e6(n)-2:3*e6(n),3*e1l(n)-2:3*el1(n))+k(16:l8,31:33,n);
K(3*e6(n)-2:3*e6(n),3*e12(n)-2:3*e12(n))=...
K(3*e6(n)-2:3*e6(n),3*e12(n)-2:3*e12(n))+k(16:18,34:36,n);
K(3*e6(n)-2:3*e6(n),3*e13(n)-2:3*e13(n))=...
K(3*e6(n)-2:3*e6(n),3*e13(n)-2:3*e13(n))+k(16:18,37:39,n);
K(3*e6(n)-2:3*e6(n),3*e14(n)-2:3*e14(n))=...
K(3*e6(n)-2:3*e6(n),3*e14(n)-2:3*e14(n))+k(l6:18,40:42,n);
K(3*e6(n)-2:3*e6(n),3*e15(n)-2:3*e15(n))=...
K(3*e6(n)-2:3*e6(n),3*el5(n)-2:3*e15(n))+k(16:18,43:45,n);
K(3*e6(n)-2:3*e6(n),3*e16(n)-2:3*e16(n))=...
K(3*e6(n)-2:3*e6(n),3*e16(n)-2:3*el6(n))+k(16:l8,46:48,n);
K(3*e6(n)-2:3*e6(n),3*e17(n)-2:3*e17(n))=...
K(3*e6(n)-2:3*e6(n),3*e17(n)-2:3*e17(n))+k(16:18,49:51,n);
K(3*e6(n)-2:3*e6(n),3*e18(n)-2:3*e18(n))=...
K(3*e6(n)-2:3*e6(n),3*e18(n)-2:3*e18(n))+k(16:18,52:54,n);
K(3*e6(n)-2:3*e6(n),3*e19(n)-2:3*e19(n))=...
K(3*e6(n)-2:3*e6(n),3*e19(n)-2:3*e19(n))+k(16:18,55:57,n);
K(3*e6(n)-2:3*e6(n),3*e20(n)-2:3*e20(n))=...
K(3*e6(n)-2:3*e6(n),3*e20(n)-2:3*e20(n))+k(16:l8,58:60,n);
K(3*e7(n)-2:3*e7(n),3*e1(n)-2:3*e1(n))=...
K(3*e7(n)-2:3*e7(n),3*e1(n)-2:3*e1(n))+k(19:21,1:3,n);
K(3*e7(n)-2:3*e7(n),3*e2(n)-2:3*e2(n))=...
K(3*e7(n)-2:3*e7(n),3*e2(n)-2:3*e2(n))+k(19:21,4:6,n);
K(3*e7(n)-2:3*e7(n),3*e3(n)-2:3*e3(n))=...
K(3*e7(n)-2:3*e7(n),3*e3(n)-2:3*e3(n))+k(19:21,7:9,n);
K(3*e7(n)-2:3*e7(n),3*e4(n)-2:3*e4(n))=...
K(3*e7(n)-2:3*e7(n),3*e4(n)-2:3*e4(n))+k(19:21,10:l2,n);
K(3*e7(n)-2:3*e7(n),3*e5(n)-2:3*e5(n))=...
K(3*e7(n)-2:3*e7(n),3*e5(n)-2:3*e5(n))+k(19:21,13:15,n);
K(3*e7(n)-2:3*e7(n),3*e6(n)-2:3*e6(n))=...
K(3*e7(n)-2:3*e7(n),3*e6(n)-2:3*e6(n))+k(19:21,16:l8,n);
K(3*e7(n)-2:3*e7(n),3*e7(n)-2:3*e7(n))=...
K(3*e7(n)-2:3*e7(n),3*e7(n)-2:3*e7(n))+k(19:21,19:21,n);
K(3*e7(n)-2:3*e7(n),3*e8(n)-2:3*e8(n))=...
K(3*e7(n)-2:3*e7(n),3*e8(n)-2:3*e8(n))+k(l9:21,22:24,n);
K(3*e7(n)-2:3*e7(n),3*e9(n)-2:3*e9(n))=...
K(3*e7(n)-2:3*e7(n),3*e9(n)-2:3*e9(n))+k(19:21,25:27,n);
K(3*e7(n)-2:3*e7(n),3*e10(n)-2:3*e10(n))=...
K(3*e7(n)-2:3*e7(n),3*e10(n)-2:3*e10(n))+k(19:21,28:30,n);
K(3*e7(n)-2:3*e7(n),3*e11(n)-2:3*e11(n))=...
K(3*e7(n)-2:3*e7(n),3*e11(n)-2:3*e11(n))+k(19:21,31:33,n);
K(3*e7(n)-2:3*e7(n),3*e12(n)-2:3*e12(n))=...
K(3*e7(n)-2:3*e7(n),3*e12(n)-2:3*e12(n))+k(19:21,34:36,n);
K(3*e7(n)-2:3*e7(n),3*e13(n)-2:3*e13(n))=...
K(3*e7(n)-2:3*e7(n),3*el3(n)-2:3*e13(n))+k(19:21,37:39,n);
K(3*e7(n)-2:3*e7(n),3*e14(n)-2:3*e14(n))=...
K(3*e7(n)-2:3*e7(n),3*e14(n)-2:3*el4(n))+k(19:21,40:42,n);
K(3*e7(n)-2:3*e7(n),3*e15(n)-2:3*e15(n))=...
K(3*e7(n)-2:3*e7(n),3*e15(n)-2:3*e15(n))+k(19:21,43:45,n);
K(3*e7(n)-2:3*e7(n),3*e16(n)-2:3*e16(n))=...
K(3*e7(n)-2:3*e7(n),3*e16(n)-2:3*e16(n))+k(19:21,46:48,n);
K(3*e7(n)-2:3*e7(n),3*e17(n)-2:3*e17(n))=...
K(3*e7(n)-2:3*e7(n),3*e17(n)-2:3*e17(n))+k(19:21,49:51,n);
K(3*e7(n)-2:3*e7(n),3*e18(n)-2:3*e18(n))=...
K(3*e7(n)-2:3*e7(n),3*e18(n)-2:3*e18(n))+k(19:21,52:54,n);
K(3*e7(n)-2:3*e7(n),3*e19(n)-2:3*e19(n))=...
K(3*e7(n)-2:3*e7(n),3*el9(n)-2:3*e19(n))+k(19:21,55:57,n);
K(3*e7(n)-2:3*e7(n),3*e20(n)-2:3*e20(n))=...
K(3*e7(n)-2:3*e7(n),3*e20(n)-2:3*e20(n))+k(19:21,58:60,n);
K(3*e8(n)-2:3*e8(n),3*e1(n)-2:3*e1(n))=...
K(3*e8(n)-2:3*e8(n),3*el(n)-2:3*e1(n))+k(22:24,1:3,n);
K(3*e8(n)-2:3*e8(n),3*e2(n)-2:3*e2(n))=...
K(3*e8(n)-2:3*e8(n),3*e2(n)-2:3*e2(n))+k(22:24,4:6,n);
K(3*e8(n)-2:3*e8(n),3*e3(n)-2:3*e3(n))=...
K(3*e8(n)-2:3*e8(n),3*e3(n)-2:3*e3(n))+k(22:24,7:9,n);
K(3*e8(n)-2:3*e8(n),3*e4(n)-2:3*e4(n))=...
K(3*e8(n)-2:3*e8(n),3*e4(n)-2:3*e4(n))+k(22:24,10:l2,n);
K(3*e8(n)-2:3*e8(n),3*e5(n)-2:3*e5(n))=...
K(3*e8(n)-2:3*e8(n),3*e5(n)-2:3*e5(n))+k(22:24,13:15,n);
K(3*e8(n)-2:3*e8(n),3*e6(n)-2:3*e6(n))=...
K(3*e8(n)-2:3*e8(n),3*e6(n)-2:3*e6(n))+k(22:24,l6:18,n);
K(3*e8(n)-2:3*e8(n),3*e7(n)-2:3*e7(n))=...
K(3*e8(n)-2:3*e8(n),3*e7(n)-2:3*e7(n))+k(22:24,19:21,n);
K(3*e8(n)-2:3*e8(n),3*e8(n)-2:3*e8(n))=...
K(3*e8(n)-2:3*e8(n),3*e8(n)-2:3*e8(n))+k(22:24,22:24,n);
K(3*e8(n)-2:3*e8(n),3*e9(n)-2:3*e9(n))=...
K(3*e8(n)-2:3*e8(n),3*e9(n)-2:3*e9(n))+k(22:24,25:27,n);
K(3*e8(n)-2:3*e8(n),3*e10(n)-2:3*e10(n))=...
K(3*e8(n)-2:3*e8(n),3*e10(n)-2:3*e10(n))+k(22:24,28:30,n);
K(3*e8(n)-2:3*e8(n),3*e11(n)-2:3*e11(n))=...
K(3*e8(n)-2:3*e8(n),3*ell(n)-2:3*e11(n))+k(22:24,31:33,n);
K(3*e8(n)-2:3*e8(n),3*e12(n)-2:3*e12(n))=...
K(3*e8(n)-2:3*e8(n),3*e12(n)-2:3*e12(n))+k(22:24,34:36,n);
K(3*e8(n)-2:3*e8(n),3*e13(n)-2:3*e13(n))=...
K(3*e8(n)-2:3*e8(n),3*el3(n)-2:3*e13(n))+k(22:24,37:39,n);
K(3*e8(n)-2:3*e8(n),3*e14(n)-2:3*e14(n))=...
K(3*e8(n)-2:3*e8(n),3*e14(n)-2:3*e14(n))+k(22:24,40:42,n);
K(3*e8(n)-2:3*e8(n),3*e15(n)-2:3*e15(n))=...
K(3*e8(n)-2:3*e8(n),3*e15(n)-2:3*e15(n))+k(22:24,43:45,n);
K(3*e8(n)-2:3*e8(n),3*e16(n)-2:3*e16(n))=...
K(3*e8(n)-2:3*e8(n),3*el6(n)-2:3*e16(n))+k(22:24,46:48,n);
K(3*e8(n)-2:3*e8(n),3*e17(n)-2:3*e17(n))=...
K(3*e8(n)-2:3*e8(n),3*e17(n)-2:3*e17(n))+k(22:24,49:51,n);
K(3*e8(n)-2:3*e8(n),3*e18(n)-2:3*e18(n))=...
K(3*e8(n)-2:3*e8(n),3*e18(n)-2:3*e18(n))+k(22:24,52:54,n);
K(3*e8(n)-2:3*e8(n),3*e19(n)-2:3*e19(n))=...
K(3*e8(n)-2:3*e8(n),3*e19(n)-2:3*e19(n))+k(22:24,55:57,n);
K(3*e8(n)-2:3*e8(n),3*e20(n)-2:3*e20(n))=...
K(3*e8(n)-2:3*e8(n),3*e20(n)-2:3*e20(n))+k(22:24,58:60,n);
K(3*e9(n)-2:3*e9(n),3*e1(n)-2:3*e1(n))=...
K(3*e9(n)-2:3*e9(n),3*e1(n)-2:3*e1(n))+k(25:27,1:3,n);
K(3*e9(n)-2:3*e9(n),3*e2(n)-2:3*e2(n))=...
K(3*e9(n)-2:3*e9(n),3*e2(n)-2:3*e2(n))+k(25:27,4:6,n);
K(3*e9(n)-2:3*e9(n),3*e3(n)-2:3*e3(n))=...
K(3*e9(n)-2:3*e9(n),3*e3(n)-2:3*e3(n))+k(25:27,7:9,n);
K(3*e9(n)-2:3*e9(n),3*e4(n)-2:3*e4(n))=...
K(3*e9(n)-2:3*e9(n),3*e4(n)-2:3*e4(n))+k(25:27,10:12,n);
K(3*e9(n)-2:3*e9(n),3*e5(n)-2:3*e5(n))=...
K(3*e9(n)-2:3*e9(n),3*e5(n)-2:3*e5(n))+k(25:27,13:15,n);
K(3*e9(n)-2:3*e9(n),3*e6(n)-2:3*e6(n))=...
K(3*e9(n)-2:3*e9(n),3*e6(n)-2:3*e6(n))+k(25:27,16:18,n);
K(3*e9(n)-2:3*e9(n),3*e7(n)-2:3*e7(n))=...
K(3*e9(n)-2:3*e9(n),3*e7(n)-2:3*e7(n))+k(25:27,19:21,n);
K(3*e9(n)-2:3*e9(n),3*e8(n)-2:3*e8(n))=...
K(3*e9(n)-2:3*e9(n),3*e8(n)-2:3*e8(n))+k(25:27,22:24,n);
K(3*e9(n)-2:3*e9(n),3*e9(n)-2:3*e9(n))=...
K(3*e9(n)-2:3*e9(n),3*e9(n)-2:3*e9(n))+k(25:27,25:27,n);
K(3*e9(n)-2:3*e9(n),3*e10(n)-2:3*e10(n))=...
K(3*e9(n)-2:3*e9(n),3*e10(n)-2:3*e10(n))+k(25:27,28:30,n);
K(3*e9(n)-2:3*e9(n),3*e11(n)-2:3*e11(n))=...
K(3*e9(n)-2:3*e9(n),3*el1(n)-2:3*e11(n))+k(25:27,31:33,n);
K(3*e9(n)-2:3*e9(n),3*e12(n)-2:3*e12(n))=...
K(3*e9(n)-2:3*e9(n),3*e12(n)-2:3*e12(n))+k(25:27,34:36,n);
K(3*e9(n)-2:3*e9(n),3*e13(n)-2:3*e13(n))=...
K(3*e9(n)-2:3*e9(n),3*e13(n)-2:3*e13(n))+k(25:27,37:39,n);
K(3*e9(n)-2:3*e9(n),3*e14(n)-2:3*e14(n))=...
K(3*e9(n)-2:3*e9(n),3*e14(n)-2:3*e14(n))+k(25:27,40:42,n);
K(3*e9(n)-2:3*e9(n),3*e15(n)-2:3*e15(n))=...
K(3*e9(n)-2:3*e9(n),3*e15(n)-2:3*e15(n))+k(25:27,43:45,n);
K(3*e9(n)-2:3*e9(n),3*e16(n)-2:3*e16(n))=...
K(3*e9(n)-2:3*e9(n),3*el6(n)-2:3*e16(n))+k(25:27,46:48,n);
K(3*e9(n)-2:3*e9(n),3*e17(n)-2:3*e17(n))=...
K(3*e9(n)-2:3*e9(n),3*e17(n)-2:3*e17(n))+k(25:27,49:51,n);
K(3*e9(n)-2:3*e9(n),3*e18(n)-2:3*e18(n))=...
K(3*e9(n)-2:3*e9(n),3*e18(n)-2:3*e18(n))+k(25:27,52:54,n);
K(3*e9(n)-2:3*e9(n),3*e19(n)-2:3*e19(n))=...
K(3*e9(n)-2:3*e9(n),3*e19(n)-2:3*e19(n))+k(25:27,55:57,n);
K(3*e9(n)-2:3*e9(n),3*e20(n)-2:3*e20(n))=...
K(3*e9(n)-2:3*e9(n),3*e20(n)-2:3*e20(n))+k(25:27,58:60,n);
K(3*e10(n)-2:3*e10(n),3*e1(n)-2:3*e1(n))=...
K(3*e10(n)-2:3*e10(n),3*el(n)-2:3*e1(n))+k(28:30,l:3,n);
K(3*e10(n)-2:3*e10(n),3*e2(n)-2:3*e2(n))=...
K(3*e10(n)-2:3*e10(n),3*e2(n)-2:3*e2(n))+k(28:30,4:6,n);
K(3*el0(n)-2:3*el0(n),3*e3(n)-2:3*e3(n))=...
K(3*e10(n)-2:3*e10(n),3*e3(n)-2:3*e3(n))+k(28:30,7:9,n);
K(3*e10(n)-2:3*el0(n),3*e4(n)-2:3*e4(n))=...
K(3*e10(n)-2:3*e10(n),3*e4(n)-2:3*e4(n))+k(28:30,10:12,n);
K(3*el0(n)-2:3*el0(n),3*e5(n)-2:3*e5(n))=...
K(3*e10(n)-2:3*e10(n),3*e5(n)-2:3*e5(n))+k(28:30,l3:15,n);
K(3*e10(n)-2:3*e10(n),3*e6(n)-2:3*e6(n))=...
K(3*el0(n)-2:3*e10(n),3*e6(n)-2:3*e6(n))+k(28:30,l6:18,n);
K(3*e10(n)-2:3*e10(n),3*e7(n)-2:3*e7(n))=...
K(3*e10(n)-2:3*e10(n),3*e7(n)-2:3*e7(n))+k(28:30,l9:21,n);
K(3*e10(n)-2:3*el0(n),3*e8(n)-2:3*e8(n))=...
K(3*e10(n)-2:3*el0(n),3*e8(n)-2:3*e8(n))+k(28:30,22:24,n);
K(3*e10(n)-2:3*el0(n),3*e9(n)-2:3*e9(n))=...
K(3*e10(n)-2:3*e10(n),3*e9(n)-2:3*e9(n))+k(28:30,25:27,n);
K(3*el0(n)-2:3*e10(n),3*e10(n)-2:3*e10(n))=...
K(3*e10(n)-2:3*e10(n),3*e10(n)-2:3*e10(n))+k(28:30,28:30,n);
K(3*e10(n)-2:3*e10(n),3*el1(n)-2:3*e11(n))=...
K(3*e10(n)-2:3*e10(n),3*e1l(n)-2:3*ell(n))+k(28:30,31:33,n);
K(3*e10(n)-2:3*e10(n),3*e12(n)-2:3*e12(n))=...
K(3*e10(n)-2:3*e10(n),3*el2(n)-2:3*e12(n))+k(28:30,34:36,n);
K(3*e10(n)-2:3*e10(n),3*e13(n)-2:3*e13(n))=...
K(3*e10(n)-2:3*e10(n),3*e13(n)-2:3*el3(n))+k(28:30,37:39,n);
K(3*e10(n)-2:3*el0(n),3*e14(n)-2:3*e14(n))=...
K(3*e10(n)-2:3*e10(n),3*e14(n)-2:3*e14(n))+k(28:30,40:42,n);
K(3*e10(n)-2:3*e10(n),3*e15(n)-2:3*e15(n))=...
K(3*e10(n)-2:3*e10(n),3*e15(n)-2:3*e15(n))+k(28:30,43:45,n);
K(3*e10(n)-2:3*e10(n),3*e16(n)-2:3*e16(n))=...
K(3*e10(n)-2:3*e10(n),3*el6(n)-2:3*el6(n))+k(28:30,46:48,n);
K(3*e10(n)-2:3*e10(n),3*e17(n)-2:3*e17(n))=...
K(3*e10(n)-2:3*e10(n),3*el7(n)-2:3*e17(n))+k(28:30,49:51,n);
K(3*e10(n)-2:3*e10(n),3*e18(n)-2:3*e18(n))=...
K(3*e10(n)-2:3*el0(n),3*e18(n)-2:3*e18(n))+k(28:30,52:54,n);
K(3*e10(n)-2:3*el0(n),3*e19(n)-2:3*e19(n))=...
K(3*e10(n)-2:3*e10(n),3*e19(n)-2:3*e19(n))+k(28:30,55:57,n);
K(3*e10(n)-2:3*el0(n),3*e20(n)-2:3*e20(n))=...
K(3*e10(n)-2:3*e10(n),3*e20(n)-2:3*e20(n))+k(28:30,58:60,n);
K(3*e11(n)-2:3*ell(n),3*e1(n)-2:3*el(n))=...
K(3*e11(n)-2:3*e11(n),3*el(n)-2:3*e1(n))+k(31:33,l:3,n);
K(3*e11(n)-2:3*ell(n),3*e2(n)-2:3*e2(n))=...
K(3*e11(n)-2:3*e11(n),3*e2(n)-2:3*e2(n))+k(31:33,4:6,n);
K(3*el1(n)-2:3*e11(n),3*e3(n)-2:3*e3(n))=...
K(3*e11(n)-2:3*e11(n),3*e3(n)-2:3*e3(n))+k(31:33,7:9,n);
K(3*e11(n)-2:3*e11(n),3*e4(n)-2:3*e4(n))=...
K(3*el1(n)-2:3*e11(n),3*e4(n)-2:3*e4(n))+k(31:33,l0:12,n);
K(3*e11(n)-2:3*e11(n),3*e5(n)-2:3*e5(n))=...
K(3*e11(n)-2:3*e11(n),3*e5(n)-2:3*e5(n))+k(31:33,13:15,n);
K(3*el1(n)-2:3*ell(n),3*e6(n)-2:3*e6(n))=...
K(3*e11(n)-2:3*e11(n),3*e6(n)-2:3*e6(n))+k(31:33,16:l8,n);
K(3*el1(n)-2:3*ell(n),3*e7(n)-2:3*e7(n))=...
K(3*e11(n)-2:3*e11(n),3*e7(n)-2:3*e7(n))+k(31:33,19:21,n);
K(3*e11(n)-2:3*e1l(n),3*e8(n)-2:3*e8(n))=...
K(3*e11(n)-2:3*e11(n),3*e8(n)-2:3*e8(n))+k(31:33,22:24,n);
K(3*el1(n)-2:3*e11(n),3*e9(n)-2:3*e9(n))=...
K(3*e11(n)-2:3*e11(n),3*e9(n)-2:3*e9(n))+k(31:33,25:27,n);
K(3*e11(n)-2:3*e11(n),3*e10(n)-2:3*e10(n))=...
K(3*e11(n)-2:3*e11(n),3*el0(n)-2:3*e10(n))+k(31:33,28:30,n);
K(3*e11(n)-2:3*el1(n),3*e11(n)-2:3*e11(n))=...
K(3*e11(n)-2:3*e11(n),3*e1l(n)-2:3*e11(n))+k(31:33,31:33,n);
K(3*e11(n)-2:3*el1(n),3*e12(n)-2:3*e12(n))=...
K(3*e11(n)-2:3*e11(n),3*e12(n)-2:3*e12(n))+k(31:33,34:36,n);
K(3*el1(n)-2:3*e11(n),3*e13(n)-2:3*e13(n))=...
K(3*e11(n)-2:3*e11(n),3*e13(n)-2:3*el3(n))+k(31:33,37:39,n);
K(3*e11(n)-2:3*e11(n),3*e14(n)*2:3*e14(n))=...
K(3*e11(n)-2:3*e11(n),3*el4(n)-2:3*el4(n))+k(31:33,40:42,n);
K(3*e11(n)-2:3*e11(n),3*e15(n)-2:3*e15(n))=...
K(3*e11(n)-2:3*e11(n),3*el5(n)-2:3*e15(n))+k(31:33,43:45,n);
K(3*e11(n)-2:3*el1(n),3*e16(n)-2:3*e16(n))=...
K(3*e11(n)-2:3*e11(n),3*e16(n)-2:3*e16(n))+k(31:33,46:48,n);
K(3*e11(n)-2:3*el1(n),3*e17(n)-2:3*e17(n))=...
K(3*e11(n)-2:3*e11(n),3*e17(n)-2:3*e17(n))+k(31:33,49:51,n);
K(3*el1(n)-2:3*e11(n),3*el8(n)-2:3*e18(n))=...
K(3*e11(n)-2:3*e11(n),3*el8(n)-2:3*e18(n))+k(31:33,52:54,n);
K(3*el1(n)-2:3*e11(n),3*e19(n)-2:3*e19(n))=...
K(3*e11(n)-2:3*e11(n),3*el9(n)-2:3*el9(n))+k(31:33,55:57,n);
K(3*e11(n)-2:3*e11(n),3*e20(n)-2:3*e20(n))=...
K(3*e11(n)-2:3*e11(n),3*e20(n)-2:3*e20(n))+k(31:33,58:60,n);
K(3*e12(n)-2:3*e12(n),3*el(n)-2:3*e1(n))=...
K(3*e12(n)-2:3*e12(n),3*e1(n)-2:3*e1(n))+k(34:36,l:3,n);
K(3*e12(n)-2:3*e12(n),3*e2(n)-2:3*e2(n))=...
K(3*e12(n)-2:3*e12(n),3*e2(n)-2:3*e2(n))+k(34:36,4:6,n);
K(3*e12(n)-2:3*e12(n),3*e3(n)-2:3*e3(n))=...
K(3*e12(n)-2:3*e12(n),3*e3(n)-2:3*e3(n))+k(34:36,7:9,n);
K(3*e12(n)-2:3*el2(n),3*e4(n)-2:3*e4(n))=...
K(3*e12(n)-2:3*e12(n),3*e4(n)-2:3*e4(n))+k(34:36,10:12,n);
K(3*el2(n)-2:3*el2(n),3*e5(n)-2:3*e5(n))=...
K(3*e12(n)-2:3*e12(n),3*e5(n)-2:3*e5(n))+k(34:36,13:15,n);
K(3*e12(n)-2:3*e12(n),3*e6(n)-2:3*e6(n))=...
K(3*e12(n)-2:3*e12(n),3*e6(n)-2:3*e6(n))+k(34:36,16:18,n);
K(3*e12(n)-2:3*e12(n),3*e7(n)-2:3*e7(n))=...
K(3*e12(n)-2:3*e12(n),3*e7(n)-2:3*e7(n))+k(34:36,19:21,n);
K(3*e12(n)-2:3*e12(n),3*e8(n)-2:3*e8(n))=...
K(3*e12(n)-2:3*e12(n),3*e8(n)-2:3*e8(n))+k(34:36,22:24,n);
K(3*el2(n)-2:3*el2(n),3*e9(n)-2:3*e9(n))=...
K(3*e12(n)-2:3*e12(n),3*e9(n)-2:3*e9(n))+k(34:36,25:27,n);
K(3*el2(n)-2:3*e12(n),3*e10(n)-2:3*e10(n))=...
K(3*e12(n)-2:3*e12(n),3*e10(n)-2:3*e10(n))+k(34:36,28:30,n);
K(3*el2(n)-2:3*e12(n),3*e11(n)-2:3*e11(n))=...
K(3*e12(n)-2:3*e12(n),3*e11(n)-2:3*e11(n))+k(34:36,31:33,n);
K(3*el2(n)-2:3*e12(n),3*e12(n)-2:3*e12(n))=...
K(3*e12(n)-2:3*e12(n),3*e12(n)-2:3*e12(n))+k(34:36,34:36,n);
K(3*el2(n)-2:3*e12(n),3*e13(n)-2:3*e13(n))=...
K(3*e12(n)-2:3*e12(n),3*e13(n)-2:3*el3(n))+k(34:36,37:39,n);
K(3*e12(n)-2:3*e12(n),3*el4(n)-2:3*e14(n))=...
K(3*e12(n)-2:3*e12(n),3*el4(n)-2:3*e14(n))+k(34:36,40:42,n);
K(3*e12(n)-2:3*e12(n),3*e15(n)-2:3*e15(n))=...
K(3*e12(n)-2:3*e12(n),3*e15(n)-2:3*e15(n))+k(34:36,43:45,n);
K(3*el2(n)-2:3*e12(n),3*e16(n)-2:3*e16(n))=...
K(3*e12(n)-2:3*e12(n),3*e16(n)-2:3*e16(n))+k(34:36,46:48,n);
K(3*e12(n)-2:3*e12(n),3*e17(n)-2:3*e17(n))=...
K(3*e12(n)-2:3*e12(n),3*e17(n)-2:3*e17(n))+k(34:36,49:51,n);
K(3*el2(n)-2:3*e12(n),3*el8(n)-2:3*e18(n))=...
K(3*e12(n)-2:3*e12(n),3*el8(n)-2:3*el8(n))+k(34:36,52:54,n);
K(3*e12(n)-2:3*e12(n),3*e19(n)-2:3*e19(n))=...
K(3*e12(n)-2:3*e12(n),3*el9(n)-2:3*e19(n))+k(34:36,55:57,n);
K(3*e12(n)-2:3*e12(n),3*e20(n)-2:3*e20(n))=...
K(3*e12(n)-2:3*e12(n),3*e20(n)-2:3*e20(n))+k(34:36,58:60,n);
K(3*e13(n)-2:3*e13(n),3*e1(n)-2:3*e1(n))=...
K(3*e13(n)-2:3*e13(n),3*el(n)-2:3*e1(n))+k(37:39,l:3,n);
K(3*e13(n)-2:3*e13(n),3*e2(n)-2:3*e2(n))=...
K(3*e13(n)-2:3*e13(n),3*e2(n)-2:3*e2(n))+k(37:39,4:6,n);
K(3*e13(n)-2:3*el3(n),3*e3(n)-2:3*e3(n))=...
K(3*e13(n)-2:3*el3(n),3*e3(n)-2:3*e3(n))+k(37:39,7:9,n);
K(3*e13(n)-2:3*el3(n),3*e4(n)-2:3*e4(n))=...
K(3*e13(n)-2:3*e13(n),3*e4(n)-2:3*e4(n))+k(37:39,10:12,n);
K(3*el3(n)-2:3*e13(n),3*e5(n)-2:3*e5(n))=...
K(3*e13(n)-2:3*e13(n),3*e5(n)-2:3*e5(n))+k(37:39,l3:15,n);
K(3*e13(n)-2:3*e13(n),3*e6(n)-2:3*e6(n))=...
K(3*el3(n)-2:3*e13(n),3*e6(n)-2:3*e6(n))+k(37:39,l6:18,n);
K(3*e13(n)-2:3*e13(n),3*e7(n)-2:3*e7(n))=...
K(3*e13(n)-2:3*e13(n),3*e7(n)-2:3*e7(n))+k(37:39,19:21,n);
K(3*el3(n)-2:3*el3(n),3*e8(n)-2:3*e8(n))=...
K(3*e13(n)-2:3*e13(n),3*e8(n)-2:3*e8(n))+k(37:39,22:24,n);
K(3*el3(n)-2:3*el3(n),3*e9(n)-2:3*e9(n))=...
K(3*e13(n)-2:3*e13(n),3*e9(n)-2:3*e9(n))+k(37:39,25:27,n);
K(3*el3(n)-2:3*e13(n),3*e10(n)-2:3*e10(n))=...
K(3*e13(n)-2:3*e13(n),3*el0(n)-2:3*el0(n))+k(37:39,28:30,n);
K(3*e13(n)-2:3*e13(n),3*e11(n)-2:3*e11(n))=...
K(3*e13(n)-2:3*e13(n),3*ell(n)-2:3*ell(n))+k(37:39,31:33,n);
K(3*e13(n)-2:3*e13(n),3*e12(n)-2:3*e12(n))=...
K(3*e13(n)-2:3*e13(n),3*el2(n)-2:3*e12(n))+k(37:39,34:36,n);
K(3*e13(n)-2:3*e13(n),3*e13(n)-2:3*e13(n))=...
K(3*e13(n)-2:3*el3(n),3*el3(n)-2:3*e13(n))+k(37:39,37:39,n);
K(3*e13(n)-2:3*el3(n),3*e14(n)-2:3*e14(n))=...
K(3*e13(n)-2:3*e13(n),3*e14(n)-2:3*e14(n))+k(37:39,40:42,n);
K(3*el3(n)-2:3*e13(n),3*e15(n)-2:3*e15(n))=...
K(3*e13(n)-2:3*e13(n),3*e15(n)-2:3*e15(n))+k(37:39,43:45,n);
K(3*e13(n)-2:3*e13(n),3*e16(n)*2:3*e16(n))=...
K(3*e13(n)-2:3*e13(n),3*el6(n)-2:3*el6(n))+k(37:39,46:48,n);
K(3*e13(n)-2:3*e13(n),3*e17(n)-2:3*e17(n))=...
K(3*e13(n)-2:3*e13(n),3*el7(n)-2:3*e17(n))+k(37:39,49:51,n);
K(3*e13(n)-2:3*e13(n),3*e1B(n)-2:3*e18(n))=...
K(3*e13(n)-2:3*e13(n),3*e18(n)-2:3*e18(n))+k(37:39,52:54,n);
K(3*e13(n)-2:3*e13(n),3*e19(n)-2:3*e19(n))=...
K(3*e13(n)-2:3*e13(n),3*e19(n)-2:3*e19(n))+k(37:39,55:57,n);
K(3*el3(n)-2:3*e13(n),3*e20(n)-2:3*e20(n))=...
K(3*e13(n)-2:3*e13(n),3*e20(n)-2:3*e20(n))+k(37:39,58:60,n);
K(3*el4(n)-2:3*el4(n),3*e1(n)-2:3*e1(n))=...
K(3*e14(n)-2:3*e14(n),3*e1(n)-2:3*e1(n))+k(40:42,1:3,n);
K(3*el4(n)-2:3*el4(n),3*e2(n)-2:3*e2(n))=...
K(3*e14(n)-2:3*e14(n),3*e2(n)-2:3*e2(n))+k(40:42,4:6,n);
K(3*e14(n)-2:3*el4(n),3*e3(n)-2:3*e3(n))=...
K(3*e14(n)-2:3*e14(n),3*e3(n)-2:3*e3(n))+k(40:42,7:9,n);
K(3*el4(n)-2:3*el4(n),3*e4(n)-2:3*e4(n))=...
K(3*e14(n)-2:3*e14(n),3*e4(n)-2:3*e4(n))+k(40:42,10:12,n);
K(3*el4(n)-2:3*el4(n),3*e5(n)-2:3*e5(n))=...
K(3*el4(n)-2:3*e14(n),3*e5(n)-2:3*e5(n))+k(40:42,l3:15,n);
K(3*e14(n)-2:3*e14(n),3*e6(n)-2:3*e6(n))=...
K(3*el4(n)-2:3*e14(n),3*e6(n)-2:3*e6(n))+k(40:42,l6:18,n);
K(3*e14(n)-2:3*e14(n),3*e7(n)-2:3*e7(n))=...
K(3*e14(n)-2:3*e14(n),3*e7(n)-2:3*e7(n))+k(40:42,l9:21,n);
K(3*e14(n)-2:3*el4(n),3*e8(n)-2:3*e8(n))=...
K(3*e14(n)-2:3*e14(n),3*e8(n)-2:3*e8(n))+k(40:42,22:24,n);
K(3*el4(n)-2:3*el4(n),3*e9(n)-2:3*e9(n))=...
K(3*e14(n)-2:3*e14(n),3*e9(n)-2:3*e9(n))+k(40:42,25:27,n);
K(3*el4(n)-2:3*el4(n),3*el0(n)-2:3*e10(n))=...
K(3*e14(n)-2:3*e14(n),3*el0(n)-2:3*el0(n))+k(40:42,28:30,n);
K(3*e14(n)-2:3*e14(n),3*ell(n)-2:3*e11(n))=...
K(3*e14(n)-2:3*e14(n),3*e1l(n)-2:3*el1(n))+k(40:42,31:33,n);
K(3*e14(n)-2:3*e14(n),3*e12(n)-2:3*e12(n))=...
K(3*e14(n)-2:3*e14(n),3*e12(n)-2:3*e12(n))+k(40:42,34:36,n);
K(3*e14(n)-2:3*el4(n),3*e13(n)-2:3*e13(n))=...
K(3*e14(n)-2:3*e14(n),3*el3(n)-2:3*e13(n))+k(40:42,37:39,n);
K(3*el4(n)-2:3*el4(n),3*e14(n)-2:3*e14(n))=...
K(3*e14(n)-2:3*e14(n),3*e14(n)-2:3*e14(n))+k(40:42,40:42,n);
K(3*el4(n)-2:3*e14(n),3*el5(n)-2:3*e15(n))=...
K(3*e14(n)-2:3*e14(n),3*el5(n)-2:3*el5(n))+k(40:42,43:45,n);
K(3*el4(n)-2:3*e14(n),3*el6(n)-2:3*e16(n))=...
K(3*e14(n)-2:3*e14(n),3*el6(n)-2:3*e16(n))+k(40:42,46:48,n);
K(3*e14(n)-2:3*e14(n),3*e17(n)-2:3*e17(n))=...
K(3*e14(n)-2:3*e14(n),3*e17(n)-2:3*e17(n))+k(40:42,49:51,n);
K(3*el4(n)-2:3*el4(n),3*e18(n)-2:3*e18(n))=...
K(3*e14(n)-2:3*e14(n),3*e18(n)-2:3*e18(n))+k(40:42,52:54,n);
K(3*e14(n)-2:3*el4(n),3*e19(n)-2:3*e19(n))=...
K(3*e14(n)-2:3*e14(n),3*e19(n)-2:3*e19(n))+k(40:42,55:57,n);
K(3*el4(n)-2:3*e14(n),3*e20(n)-2:3*e20(n))=...
K(3*e14(n)-2:3*e14(n),3*e20(n)-2:3*e20(n))+k(40:42,58:60,n);
K(3*el5(n)-2:3*el5(n),3*el(n)-2:3*e1(n))=...
K(3*e15(n)-2:3*e15(n),3*e1(n)-2:3*e1(n))+k(43:45,1:3,n);
K(3*el5(n)-2:3*el5(n),3*e2(n)-2:3*e2(n))=...
K(3*e15(n)-2:3*e15(n),3*e2(n)-2:3*e2(n))+k(43:45,4:6,n);
K(3*e15(n)-2:3*e15(n),3*e3(n)-2:3*e3(n))=...
K(3*e15(n)-2:3*e15(n),3*e3(n)-2:3*e3(n))+k(43:45,7:9,n);
K(3*e15(n)-2:3*e15(n),3*e4(n)-2:3*e4(n))=...
K(3*e15(n)-2:3*e15(n),3*e4(n)-2:3*e4(n))+k(43:45,l0:12,n);
K(3*e15(n)-2:3*el5(n),3*e5(n)-2:3*e5(n))=...
K(3*e15(n)-2:3*el5(n),3*e5(n)-2:3*e5(n))+k(43:45,l3:15,n);
K(3*e15(n)-2:3*el5(n),3*e6(n)-2:3*e6(n))=...
K(3*e15(n)-2:3*e15(n),3*e6(n)-2:3*e6(n))+k(43:45,16:18,n);
K(3*el5(n)-2:3*e15(n),3*e7(n)-2:3*e7(n))=...
K(3*e15(n)-2:3*e15(n),3*e7(n)-2:3*e7(n))+k(43:45,19:21,n);
K(3*e15(n)-2:3*e15(n),3*e8(n)-2:3*e8(n))=...
K(3*el5(n)-2:3*e15(n),3*e8(n)-2:3*e8(n))+k(43:45,22:24,n);
K(3*e15(n)-2:3*e15(n),3*e9(n)-2:3*e9(n))=...
K(3*e15(n)-2:3*e15(n),3*e9(n)-2:3*e9(n))+k(43:45,25:27,n);
K(3*e15(n)-2:3*e15(n),3*e10(n)-2:3*e10(n))=...
K(3*e15(n)-2:3*e15(n),3*e10(n)-2:3*e10(n))+k(43:45,28:30,n);
K(3*el5(n)-2:3*el5(n),3*e11(n)-2:3*e11(n))=...
K(3*e15(n)-2:3*e15(n),3*e11(n)-2:3*e11(n))+k(43:45,31:33,n);
K(3*e15(n)-2:3*e15(n),3*e12(n)-2:3*e12(n))=...
K(3*e15(n)-2:3*e15(n),3*e12(n)-2:3*e12(n))+k(43:45,34:36,n);
K(3*el5(n)-2:3*e15(n),3*el3(n)-2:3*e13(n))=...
K(3*e15(n)-2:3*e15(n),3*el3(n)-2:3*el3(n))+k(43:45,37:39,n);
K(3*e15(n)-2:3*e15(n),3*el4(n)-2:3*e14(n))=...
K(3*e15(n)-2:3*e15(n),3*e14(n)-2:3*e14(n))+k(43:45,40:42,n);
K(3*e15(n)-2:3*e15(n),3*e15(n)-2:3*e15(n))=...
K(3*e15(n)-2:3*el5(n),3*e15(n)-2:3*e15(n))+k(43:45,43:45,n);
K(3*e15(n)-2:3*el5(n),3*e16(n)-2:3*e16(n))=...
K(3*e15(n)-2:3*e15(n),3*el6(n)-2:3*e16(n))+k(43:45,46:48,n);
K(3*el5(n)-2:3*el5(n),3*e17(n)-2:3*e17(n))=...
K(3*e15(n)-2:3*e15(n),3*e17(n)-2:3*e17(n))+k(43:45,49:51,n);
K(3*el5(n)-2:3*el5(n),3*e18(n)-2:3*e18(n))=...
K(3*e15(n)-2:3*e15(n),3*el8(n)-2:3*el8(n))+k(43:45,52:54,n);
K(3*e15(n)-2:3*e15(n),3*el9(n)-2:3*e19(n))=...
K(3*e15(n)-2:3*e15(n),3*el9(n)-2:3*el9(n))+k(43:45,55:57,n);
K(3*e15(n)-2:3*e15(n),3*e20(n)-2:3*e20(n))=...
K(3*e15(n)-2:3*e15(n),3*e20(n)-2:3*e20(n))+k(43:45,58:60,n);
K(3*e16(n)-2:3*e16(n),3*el(n)-2:3*e1(n))=...
K(3*e16(n)-2:3*e16(n),3*el(n)-2:3*e1(n))+k(46:48,l:3,n);
K(3*e16(n)-2:3*e16(n),3*e2(n)-2:3*e2(n))=...
K(3*e16(n)-2:3*e16(n),3*e2(n)-2:3*e2(n))+k(46:48,4:6,n);
K(3*el6(n)-2:3*el6(n),3*e3(n)-2:3*e3(n))=...
K(3*e16(n)-2:3*el6(n),3*e3(n)-2:3*e3(n))+k(46:48,7:9,n);
K(3*el6(n)-2:3*el6(n),3*e4(n)-2:3*e4(n))=...
K(3*e16(n)-2:3*e16(n),3*e4(n)-2:3*e4(n))+k(46:48,10:12,n);
K(3*e16(n)-2:3*e16(n),3*e5(n)-2:3*e5(n))=...
K(3*el6(n)-2:3*e16(n),3*e5(n)-2:3*e5(n))+k(46:48,l3:15,n);
K(3*e16(n)-2:3*e16(n),3*e6(n)-2:3*e6(n))=...
K(3*e16(n)-2:3*e16(n),3*e6(n)-2:3*e6(n))+k(46:48,l6:18,n);
K(3*e16(n)-2:3*e16(n),3*e7(n)-2:3*e7(n))=...
K(3*e16(n)-2:3*e16(n),3*e7(n)-2:3*e7(n))+k(46:48,l9:21,n);
K(3*el6(n)-2:3*el6(n),3*e8(n)-2:3*e8(n))=...
K(3*e16(n)-2:3*el6(n),3*e8(n)-2:3*e8(n))+k(46:48,22:24,n);
K(3*e16(n)-2:3*el6(n),3*e9(n)-2:3*e9(n))=...
K(3*e16(n)-2:3*e16(n),3*e9(n)-2:3*e9(n))+k(46:48,25:27,n);
K(3*el6(n)-2:3*e16(n),3*e10(n)-2:3*e10(n))=...
K(3*e16(n)-2:3*e16(n),3*e10(n)-2:3*el0(n))+k(46:48,28:30,n);
K(3*e16(n)-2:3*e16(n),3*e1l(n)-2:3*e11(n))=...
K(3*e16(n)-2:3*e16(n),3*el1(n)-2:3*ell(n))+k(46:48,31:33,n);
K(3*e16(n)-2:3*e16(n),3*e12(n)-2:3*e12(n))=...
K(3*e16(n)-2:3*e16(n),3*e12(n)-2:3*e12(n))+k(46:48,34:36,n);
K(3*el6(n)-2:3*el6(n),3*e13(n)-2:3*e13(n))=...
K(3*e16(n)-2:3*e16(n),3*e13(n)-2:3*e13(n))+k(46:48,37:39,n);
K(3*el6(n)-2:3*el6(n),3*e14(n)-2:3*e14(n))=...
K(3*e16(n)-2:3*e16(n),3*e14(n)-2:3*e14(n))+k(46:48,40:42,n);
K(3*el6(n)-2:3*e16(n),3*el5(n)-2:3*e15(n))=...
K(3*e16(n)-2:3*e16(n),3*el5(n)-2:3*e15(n))+k(46:48,43:45,n);
K(3*e16(n)-2:3*e16(n),3*el6(n)-2:3*e16(n))=...
K(3*e16(n)-2:3*e16(n),3*el6(n)-2:3*el6(n))+k(46:48,46:48,n);
K(3*e16(n)-2:3*e16(n),3*e17(n)-2:3*e17(n))=...
K(3*e16(n)-2:3*e16(n),3*e17(n)-2:3*e17(n))+k(46:48,49:51,n);
K(3*e16(n)-2:3*el6(n),3*e18(n)-2:3*e18(n))=...
K(3*e16(n)-2:3*e16(n),3*e18(n)-2:3*e18(n))+k(46:48,52:54,n);
K(3*el6(n)-2:3*el6(n),3*e19(n)-2:3*e19(n))=...
K(3*e16(n)-2:3*e16(n),3*e19(n)-2:3*e19(n))+k(46:48,55:57,n);
K(3*el6(n)-2:3*e16(n),3*e20(n)-2:3*e20(n))=...
K(3*e16(n)-2:3*e16(n),3*e20(n)-2:3*e20(n))+k(46:48,58:60,n);
K(3*el7(n)-2:3*el7(n),3*e1(n)-2:3*e1(n))=...
K(3*e17(n)-2:3*e17(n),3*e1(n)-2:3*e1(n))+k(49:51,1:3,n);
K(3*e17(n)-2:3*e17(n),3*e2(n)-2:3*e2(n))=...
K(3*e17(n)-2:3*e17(n),3*e2(n)-2:3*e2(n))+k(49:51,4:6,n);
K(3*e17(n)-2:3*e17(n),3*e3(n)-2:3*e3(n))=...
K(3*e17(n)-2:3*e17(n),3*e3(n)-2:3*e3(n))+k(49:51,7:9,n);
K(3*e17(n)-2:3*e17(n),3*e4(n)-2:3*e4(n))=...
K(3*e17(n)-2:3*e17(n),3*e4(n)-2:3*e4(n))+k(49:51,10:12,n);
K(3*e17(n)-2:3*e17(n),3*e5(n)-2:3*e5(n))=...
K(3*e17(n)-2:3*e17(n),3*e5(n)-2:3*e5(n))+k(49:51,13:15,n);
K(3*el7(n)-2:3*el7(n),3*e6(n)-2:3*e6(n))=...
K(3*e17(n)-2:3*e17(n),3*e6(n)-2:3*e6(n))+k(49:51,16:18,n);
K(3*el7(n)-2:3*el7(n),3*e7(n)-2:3*e7(n))=...
K(3*e17(n)-2:3*e17(n),3*e7(n)-2:3*e7(n))+k(49:51,19:21,n);
K(3*e17(n)-2:3*el7(n),3*e8(n)-2:3*e8(n))=...
K(3*e17(n)-2:3*e17(n),3*e8(n)-2:3*e8(n))+k(49:51,22:24,n);
K(3*el7(n)-2:3*el7(n),3*e9(n)-2:3*e9(n))=...
K(3*e17(n)-2:3*e17(n),3*e9(n)-2:3*e9(n))+k(49:51,25:27,n);
K(3*el7(n)-2:3*e17(n),3*e10(n)-2:3*e10(n))=...
K(3*e17(n)-2:3*e17(n),3*e10(n)-2:3*el0(n))+k(49:51,28:30,n);
K(3*e17(n)-2:3*e17(n),3*ell(n)-2:3*e11(n))=...
K(3*e17(n)-2:3*e17(n),3*ell(n)-2:3*e1l(n))+k(49:51,31:33,n);
K(3*e17(n)-2:3*e17(n),3*e12(n)-2:3*e12(n))=...
K(3*e17(n)-2:3*e17(n),3*e12(n)-2:3*e12(n))+k(49:51,34:36,n);
K(3*e17(n)-2:3*el7(n),3*e13(n)-2:3*e13(n))=...
K(3*e17(n)-2:3*e17(n),3*e13(n)-2:3*e13(n))+k(49:51,37:39,n);
K(3*e17(n)-2:3*el7(n),3*e14(n)-2:3*e14(n))=...
K(3*e17(n)-2:3*e17(n),3*e14(n)-2:3*e14(n))+k(49:51,40:42,n);
K(3*el7(n)-2:3*e17(n),3*el5(n)-2:3*e15(n))=...
K(3*e17(n)-2:3*e17(n),3*el5(n)-2:3*el5(n))+k(49:51,43:45,n);
K(3*e17(n)-2:3*e17(n),3*e16(n)-2:3*e16(n))=...
K(3*e17(n)-2:3*e17(n),3*e16(n)-2:3*e16(n))+k(49:51,46:48,n);
K(3*e17(n)-2:3*e17(n),3*e17(n)-2:3*e17(n))=...
K(3*e17(n)-2:3*e17(n),3*el7(n)-2:3*e17(n))+k(49:51,49:51,n);
K(3*e17(n)-2:3*el7(n),3*e18(n)-2:3*e18(n))=...
K(3*e17(n)-2:3*e17(n),3*e18(n)-2:3*e18(n))+k(49:51,52:54,n);
K(3*e17(n)-2:3*el7(n),3*e19(n)-2:3*e19(n))=...
K(3*e17(n)-2:3*e17(n),3*e19(n)-2:3*e19(n))+k(49:51,55:57,n);
K(3*el7(n)-2:3*e17(n),3*e20(n)-2:3*e20(n))=...
K(3*e17(n)-2:3*e17(n),3*e20(n)-2:3*e20(n))+k(49:51,58:60,n);
K(3*e18(n)-2:3*el8(n),3*el(n)-2:3*e1(n))=...
K(3*e18(n)-2:3*e18(n),3*e1(n)-2:3*e1(n))+k(52:54,1:3,n);
K(3*el8(n)-2:3*e18(n),3*e2(n)-2:3*e2(n))=...
K(3*e18(n)-2:3*e18(n),3*e2(n)-2:3*e2(n))+k(52:54,4:6,n);
K(3*e18(n)-2:3*e18(n),3*e3(n)-2:3*e3(n))=...
K(3*e18(n)-2:3*e18(n),3*e3(n)-2:3*e3(n))+k(52:54,7:9,n);
K(3*e18(n)-2:3*e18(n),3*e4(n)-2:3*e4(n))=...
K(3*e18(n)-2:3*e18(n),3*e4(n)-2:3*e4(n))+k(52:54,10:12,n);
K(3*el8(n)-2:3*el8(n),3*e5(n)-2:3*e5(n))=...
K(3*e18(n)-2:3*e18(n),3*e5(n)-2:3*e5(n))+k(52:54,13:15,n);
K(3*el8(n)-2:3*el8(n),3*e6(n)-2:3*e6(n))=...
K(3*e18(n)-2:3*e18(n),3*e6(n)-2:3*e6(n))+k(52:54,16:18,n);
K(3*e18(n)-2:3*e18(n),3*e7(n)-2:3*e7(n))=...
K(3*e18(n)-2:3*e18(n),3*e7(n)-2:3*e7(n))+k(52:54,l9:21,n);
K(3*el8(n)-2:3*e18(n),3*e8(n)-2:3*e8(n))=...
K(3*e18(n)-2:3*e18(n),3*e8(n)-2:3*e8(n))+k(52:54,22:24,n);
K(3*e18(n)-2:3*e18(n),3*e9(n)-2:3*e9(n))=...
K(3*e18(n)-2:3*e18(n),3*e9(n)-2:3*e9(n))+k(52:54,25:27,n);
K(3*e18(n)-2:3*el8(n),3*e10(n)-2:3*e10(n))=...
K(3*e18(n)-2:3*e18(n),3*e10(n)-2:3*e10(n))+k(52:54,28:30,n);
K(3*e18(n)-2:3*el8(n),3*e11(n)-2:3*e11(n))=...
K(3*e18(n)-2:3*e18(n),3*e11(n)-2:3*e11(n))+k(52:54,31:33,n);
K(3*el8(n)-2:3*e18(n),3*e12(n)-2:3*e12(n))=...
K(3*e18(n)-2:3*e18(n),3*e12(n)-2:3*e12(n))+k(52:54,34:36,n);
K(3*e18(n)-2:3*e18(n),3*e13(n)*2:3*e13(n))=...
K(3*e18(n)-2:3*e18(n),3*el3(n)-2:3*el3(n))+k(52:54,37:39,n);
K(3*e18(n)-2:3*e18(n),3*e14(n)-2:3*e14(n))=...
K(3*e18(n)-2:3*e18(n),3*el4(n)-2:3*e14(n))+k(52:54,40:42,n);
K(3*e18(n)-2:3*e18(n),3*e15(n)-2:3*e15(n))=...
K(3*e18(n)-2:3*e18(n),3*e15(n)-2:3*e15(n))+k(52:54,43:45,n);
K(3*e18(n)-2:3*el8(n),3*e16(n)-2:3*e16(n))=...
K(3*e18(n)-2:3*e18(n),3*e16(n)-2:3*e16(n))+k(52:54,46:48,n);
K(3*el8(n)-2:3*e18(n),3*e17(n)-2:3*e17(n))=...
K(3*e18(n)-2:3*e18(n),3*e17(n)-2:3*e17(n))+k(52:54,49:51,n);
K(3*e18(n)-2:3*e18(n),3*e18(n)-2:3*e18(n))=...
K(3*e18(n)-2:3*e18(n),3*el8(n)-2:3*el8(n))+k(52:54,52:54,n);
K(3*e18(n)-2:3*e18(n),3*e19(n)-2:3*e19(n))=...
K(3*e18(n)-2:3*e18(n),3*el9(n)-2:3*e19(n))+k(52:54,55:57,n);
K(3*e18(n)-2:3*e18(n),3*e20(n)-2:3*e20(n))=...
K(3*e18(n)-2:3*e18(n),3*e20(n)-2:3*e20(n))+k(52:54,58:60,n);
K(3*e19(n)-2:3*el9(n),3*e1(n)-2:3*e1(n))=...
K(3*e19(n)-2:3*e19(n),3*e1(n)-2:3*el(n))+k(55:57,1:3,n);
K(3*el9(n)-2:3*el9(n),3*e2(n)-2:3*e2(n))=...
K(3*e19(n)-2:3*e19(n),3*e2(n)-2:3*e2(n))+k(55:57,4:6,n);
K(3*e19(n)-2:3*el9(n),3*e3(n)-2:3*e3(n))=...
K(3*e19(n)-2:3*e19(n),3*e3(n)-2:3*e3(n))+k(55:57,7:9,n);
K(3*e19(n)-2:3*e19(n),3*e4(n)-2:3*e4(n))=...
K(3*el9(n)-2:3*e19(n),3*e4(n)-2:3*e4(n))+k(55:57,l0:12,n);
K(3*e19(n)-2:3*e19(n),3*e5(n)-2:3*e5(n))=...
K(3*e19(n)-2:3*e19(n),3*e5(n)-2:3*e5(n))+k(55:57,l3:15,n);
K(3*e19(n)-2:3*el9(n),3*e6(n)-2:3*e6(n))=...
K(3*e19(n)-2:3*e19(n),3*e6(n)-2:3*e6(n))+k(55:57,l6:18,n);
K(3*e19(n)-2:3*el9(n),3*e7(n)-2:3*e7(n))=...
K(3*e19(n)-2:3*e19(n),3*e7(n)-2:3*e7(n))+k(55:57,19:21,n);
K(3*el9(n)-2:3*el9(n),3*e8(n)-2:3*e8(n))=...
K(3*el9(n)-2:3*e19(n),3*e8(n)-2:3*e8(n))+k(55:57,22:24,n);
K(3*e19(n)-2:3*e19(n),3*e9(n)-2:3*e9(n))=...
K(3*e19(n)-2:3*e19(n),3*e9(n)-2:3*e9(n))+k(55:57,25:27,n);
K(3*e19(n)-2:3*e19(n),3*e10(n)-2:3*e10(n))=...
K(3*e19(n)-2:3*e19(n),3*e10(n)-2:3*e10(n))+k(55:57,28:30,n);
K(3*e19(n)-2:3*el9(n),3*e11(n)-2:3*e11(n))=...
K(3*e19(n)-2:3*e19(n),3*e11(n)-2:3*e11(n))+k(55:57,31:33,n);
K(3*e19(n)-2:3*el9(n),3*e12(n)-2:3*e12(n))=...
K(3*e19(n)-2:3*e19(n),3*e12(n)-2:3*e12(n))+k(55:57,34:36,n);
K(3*el9(n)-2:3*e19(n),3*e13(n)-2:3*e13(n))=...
K(3*e19(n)-2:3*e19(n),3*e13(n)-2:3*el3(n))+k(55:57,37:39,n);
K(3*e19(n)-2:3*e19(n),3*el4(n)-2:3*e14(n))=...
K(3*e19(n)-2:3*e19(n),3*el4(n)-2:3*e14(n))+k(55:57,40:42,n);
K(3*e19(n)-2:3*e19(n),3*e15(n)-2:3*e15(n))=...
K(3*e19(n)-2:3*e19(n),3*e15(n)-2:3*e15(n))+k(55:57,43:45,n);
K(3*e19(n)-2:3*el9(n),3*e16(n)-2:3*e16(n))=...
K(3*e19(n)-2:3*e19(n),3*e16(n)-2:3*e16(n))+k(55:57,46:48,n);
K(3*e19(n)-2:3*el9(n),3*e17(n)-2:3*e17(n))=...
K(3*e19(n)-2:3*e19(n),3*e17(n)-2:3*e17(n))+k(55:57,49:51,n);
K(3*el9(n)-2:3*e19(n),3*el8(n)-2:3*e18(n))=...
K(3*e19(n)-2:3*e19(n),3*el8(n)-2:3*el8(n))+k(55:57,52:54,n);
K(3*el9(n)-2:3*e19(n),3*e19(n)-2:3*e19(n))=...
K(3*e19(n)-2:3*e19(n),3*el9(n)-2:3*e19(n))+k(55:57,55:57,n);
K(3*e19(n)-2:3*e19(n),3*e20(n)-2:3*e20(n))=...
K(3*e19(n)-2:3*e19(n),3*e20(n)-2:3*e20(n))+k(55:57,58:60,n);
K(3*e20(n)-2:3*e20(n),3*e1(n)-2:3*e1(n))=...
K(3*e20(n)-2:3*e20(n),3*el(n)-2:3*e1(n))+k(58:60,l:3,n);
K(3*e20(n)-2:3*e20(n),3*e2(n)-2:3*e2(n))=...
K(3*e20(n)-2:3*e20(n),3*e2(n)-2:3*e2(n))+k(58:60,4:6,n);
K(3*e20(n)-2:3*e20(n),3*e3(n)-2:3*e3(n))=...
K(3*e20(n)-2:3*e20(n),3*e3(n)-2:3*e3(n))+k(58:60,7:9,n);
K(3*e20(n)-2:3*e20(n),3*e4(n)-2:3*e4(n))=...
K(3*e20(n)-2:3*e20(n),3*e4(n)-2:3*e4(n))+k(58:60,10:12,n);
K(3*e20(n)-2:3*e20(n),3*e5(n)-2:3*e5(n))=...
K(3*e20(n)-2:3*e20(n),3*e5(n)-2:3*e5(n))+k(58:60,13:15,n);
K(3*e20(n)-2:3*e20(n),3*e6(n)-2:3*e6(n))=...
K(3*e20(n)-2:3*e20(n),3*e6(n)-2:3*e6(n))+k(58:60,l6:18,n);
K(3*e20(n)-2:3*e20(n),3*e7(n)-2:3*e7(n))=...
K(3*e20(n)-2:3*e20(n),3*e7(n)-2:3*e7(n))+k(58:60,19:21,n);
K(3*e20(n)-2:3*e20(n),3*e8(n)-2:3*e8(n))=...
K(3*e20(n)-2:3*e20(n),3*e8(n)-2:3*e8(n))+k(58:60,22:24,n);
K(3*e20(n)-2:3*e20(n),3*e9(n)-2:3*e9(n))=...
K(3*e20(n)-2:3*e20(n),3*e9(n)-2:3*e9(n))+k(58:60,25:27,n);
K(3*e20(n)-2:3*e20(n),3*e10(n)-2:3*e10(n))=...
K(3*e20(n)-2:3*e20(n),3*e10(n)-2:3*e10(n))+k(58:60,28:30,n);
K(3*e20(n)-2:3*e20(n),3*e11(n)-2:3*e11(n))=...
K(3*e20(n)-2:3*e20(n),3*ell(n)-2:3*ell(n))+k(58:60,31:33,n);
K(3*e20(n)-2:3*e20(n),3*e12(n)-2:3*e12(n))=...
K(3*e20(n)-2:3*e20(n),3*el2(n)-2:3*e12(n))+k(58:60,34:36,n);
K(3*e20(n)-2:3*e20(n),3*e13(n)-2:3*e13(n))=...
K(3*e20(n)-2:3*e20(n),3*el3(n)-2:3*e13(n))+k(58:60,37:39,n);
K(3*e20(n)-2:3*e20(n),3*e14(n)-2:3*e14(n))=...
K(3*e20(n)-2:3*e20(n),3*e14(n)-2:3*e14(n))+k(58:60,40:42,n);
K(3*e20(n)-2:3*e20(n),3*e15(n)-2:3*el5(n))=...
K(3*e20(n)-2:3*e20(n),3*e15(n)-2:3*e15(n))+k(58:60,43:45,n);
K(3*e20(n)-2:3*e20(n),3*el6(n)-2:3*e16(n))=...
K(3*e20(n)-2:3*e20(n),3*e16(n)-2:3*e16(n))+k(58:60,46:48,n);
K(3*e20(n)-2:3*e20(n),3*e17(n)-2:3*e17(n))=...
K(3*e20(n)-2:3*e20(n),3*e17(n)-2:3*e17(n))+k(58:60,49:51,n);
K(3*e20(n)-2:3*e20(n),3*e18(n)-2:3*e18(n))=...
K(3*e20(n)-2:3*e20(n),3*e18(n)-2:3*e18(n))+k(58:60,52:54,n);
K(3*e20(n)-2:3*e20(n),3*e19(n)-2:3*e19(n))=...
K(3*e20(n)-2:3*e20(n),3*e19(n)-2:3*e19(n))+k(58:60,55:57,n);
K(3*e20(n)-2:3*e20(n),3*e20(n)-2:3*e20(n))=...
K(3*e20(n)-2:3*e20(n),3*e20(n)-2:3*e20(n))+k(58:60,58:60,n);
end
%Definition of primary external nodal forces vector:
F=zeros(3*NN,1);
for i=1:NF
r=3*Forc(i,1);
F(r-2)=Forc(i,2);
F(r-1)=Forc(i,3);
F(r)=Forc(i,4);
end
%Elmination of.rows and columns of Kematrix concern to supports:
S=K;
for i=1:NC
r=3*Cons(i,1);
for j=2:4
if Cons(i,j)==l
S(r-4+j,:)=0; S(:,r-4+j)=0; S(r-4+j,r-4+j)=1;
end
end
end
%Solution of equation "{E7=[S](di}" by Gauss elimination method:
n=3*NN;
%Creation of upper triangular matrix
Sum=0;
for j=1:n-1
if S(j,j)==0
k=j;
for k=k+1:n
if S(k,j)==0
continue
end
break
end
B=S(j,:); C=F(j);
S(j,:)=S(k,:); F(j)=F(k);
S(k,:)=B; F(k)=C;
end
for i=1+Sum:n-1
L=S(i+1,j)/S(j:j);
S(i+1,:)=S(i+1,:)-L*S(j,:);
F(i+1)=F(i+1)-L*F(j);
end
Sum=Sum+1;
end
%solution of equations
 
di=zeros(3*NN,1);
 
di(n)=F(n)/S(n,n);
 
for i=n-1:-1:1
 
sum=0;
 
for j=i+1:n
sum=sum+S(i,j)*di(j);
 
end
 
di(i)=(1/S(i,i))*(F(i)-sum);
 
end
 

 
%Creation of external nodal fbrces vector:
 
W=K*di;
 
%Calculation of stress in Gaussian points far each elanent:
 
for n=l:NE
 
dis(:,l,n)=[di(3*e1(n)-2);di(3*e1(n)-1);di(3*el(n));...
di(3*e2(n)-2);di(3*e2(n)-1);di(3*e2(n));...
di(3*e3(n)-2);di(3*e3(n)-1);di(3*e3(n));...
di(3*e4(n)-2);di(3*e4(n)-1);di(3*e4(n));...
di(3*e5(n)-2);di(3*e5(n)-1);di(3*e5(n));...
di(3*e6(n)-2);di(3*e6(n)-1);di(3*e6(n));...
di(3*e7(n)-2);di(3*e7(n)-1);di(3*e7(n));...
di(3*e8(n)-2);di(3*e8(n)-1);di(3*e8(n));...
di(3*e9(n)-2);di(3*e9(n)-1);di(3*e9(n));...
di(3*el0(n)-2);di(3*e10(n)-1);di(3*e10(n));...
di(3*ell(n)-2);di(3*el1(n)-1);di(3*e11(n));...
di(3*el2(n)-2);di(3*e12(n)-1);di(3*e12(n));...
di(3*e13(n)-2);di(3*el3(n)-1);di(3*e13(n));...
di(3*e14(n)-2);di(3*e14(n)-1);di(3*e14(n));...
di(3*e15(n)-2);di(3*e15(n)-1);di(3*el5(n));...
di(3*e16(n)-2);di(3*e16(n)-1);di(3*e16(n));...
di(3*e17(n)-2);di(3*e17(n)-1);di(3*e17(n));...
di(3*el8(n)-2);di(3*el8(n)-1);di(3*e18(n));...
di(3*e19(n)-2);di(3*e19(n)-1);di(3*e19(n));...
di(3*e20(n)-2);di(3*e20(n)-1);di(3*e20(n))];
Sig1(:,1,n)=D*B1(:,:,n)*dis(:,1,n);
Sig2(:,1,n)=D*B2(:,:,n)*dis(:,1,n);
Sig3(:,1,n)=D*B3(:,:,n)*dis(:,l,n);
Sig4(:,1,n)=D*B4(:,:,n)*dis(:,1,n);
Sig5(:,1,n)=D*B5(:,:,n)*dis(:,1,n);
Sig6(:,1,n)=D*B6(:,:,n)*dis(:,1,n);
Sig7(:,1,n)=D*B7(:,:,n)*dis(:,1,n);
Sig8(:,1,n)=D*B8(:,:,n)*dis(:,l,n);
 
end
 
%Transportation of stress from Gaussian points to element nodes (inteppolation):
 
a=(71*sqrt(627)/2166);
 
b=(7*sqrt(627)/1083);
 
c=(11*sqrt(627)/722);
 
d=(4*sqrt(627)/361);
 
e=(15*sqrt(627)/361);
 
T=[(46/57+a) (-7/57-b) (-7/114+b) (-7/57-b) (-7/57-b) (-7/114+b) (85/114-a) (-7/114+b)...
(-7/114-b) (85/l14+a) (-7/114-b) (-7/57+b) (-7/57+b) (-7/114-b) (-7/57+b) (46/57-a)...
(-20/57+c) (1/6-c) (59/114+e) (1/6-c) (59/57-e) (-20/57+c) (1/6-c) (-20/57+c)...
(-7/114-b) (-7/57+b) (-7/114-b) (85/l14+a) (-7/57+b) (46/57-a) (-7/57+b) (-7/114-b)...
(-20/57+c) (1/6-C) (59/114-d) (1/6-c) (59/57+d) (-20/57+c) (1/6-C) (-20/57+c)...
(-7/114+b) (-7/57-b) (-7/114+b) (85/114-a) (-7/57-b) (46/57+a) (-7/57-b) (-7/114+b)...
(46/57-a) (-7/57+b) (-7/114-b) (-7/57+b) (-7/57+b) (-7/114-b) (85/114+a) (-7/114-b)...
(-7/114+b) (85/114-a) (-7/114+b) (-7/57-b) (-7/57-b) (-7/114+b) (-7/57-b) (46/57+a)];
for n=1:NE
Sig_X(:,1,n)=T*[Sig1(1,1,n);Sig2(1,1,n);Sig3(1,1,n);Sig4(1,1,n);...
sig5(1,1,n);sig6(1,1,n);sig7(1,1,n);sig8(1,1,n)];
Sig_Y(:,1,n)=T*[Sig1(2,1,n);Sig2(2,1,n);Sig3(2,1,n);Sig4(2,1,n);...
Sig5(2,1,n);Sig6(2,1,n);Sig7(2,1,n);Sig8(2,1,n)];
Sig_Z(:,1,n)=T*[Sig1(3,1,n);Sig2(3,1,n);Sig3(3,1,n);Sig4(3,1,n);...
Sig5(3,1,n);Sig6(3,1,n);Sig7(3,1,n);Sig8(3,1,n)];
Tau_XY(:,1,n)=T*[Sig1(4,1,n);Sig2(4,1,n);Sig3(4,1,n);Sig4(4,1,n);...
Sig5(4,1,n);Sig6(4,1,n);Sig7(4,1,n);Sig8(4,1,n)];
Tau_YZ(:,1,n)=T*[Sig1(5,1,n);Sig2(5,1,n);Sig3(5,1,n);Sig4(5,1,n);...
Sig5(5,1,n);Sig6(5,1,n);Sig7(5,1,n);Sig8(5,1,n)];
Tau_XZ(:,1,n)=T*[sigl(6,1,n);Sig2(6,1,n);Sig3(6,1,n);Sig4(6,1,n);...
Sig5(6,l,n);Sig6(6,1,n);Sig7(6,1,n);Sig8(6,1,n)];
Sigma_1(:,1,n)=[Sig_X(1,1,n);Sig_Y(1,l,n);Sig_Z(1,l,n);...
Tau_XY(1,1,n);Tau_YZ(1,1,n);Tau_XZ(1,l,n)];
Sigma_2(:,1,n)=[Sig_X(2,1,n);Sig_Y(2,l,n);Sig_Z(2,l,n);...
Tau_XY(2,1,n);Tau_YZ(2,1,n);Tau_XZ(2,l,n)];
Sigma_3(:,1,n)=[Sig_X(3,l,n);Sig_Y(3,1,n);Sig_Z(3,l,n);...
Tau_XY(3,1,n);Tau_XZ(3,1,n);Tau_XZ(3,1,n)];
Sigma_4(:,1,n)=[Sig_X(4,l,n);Sig_Y(4,1,n);Sig_Z(4,1,n);...
Tau_XY(4,1,n);Tau_XZ(4,1,n);Tau_XZ(4,1,n)];
Sigma_5(:,1,n)=[Sig_X(5,1,n);Sig_Y(5,1,n);Sig_Z(5,1,n);...
Tau_XY(5,1,n);Tau_YZ(5,l,n);Tau_XZ(5,1,n)];
Sigma_6(:,1,n)=[Sig_X(6,1,n);Sig_Y(6,1,n);Sig_Z(6,l,n);...
Tau_XY(6,1,n);Tau_XZ(6,1,n);Tau_XZ(6,l,n)];
Sigma_7(:,1,n)=[Sig_X(7,1,n);Sig_Y(7,l,n);Sig_Z(7,1,n);...
Tau_XY(7,1,n);Tau_XZ(7,1,n);Tau_XZ(7,l,n)];
Sigma_8(:,1,n)=[Sig_X(8,1,n);Sig_Y(8,1,n);Sig_Z(8,1,n);...
Tau_XY(8,1,n);Tau_XZ(8,1,n);Tau_XZ(8,1,n)];
end
%Ana1ysis results:
disp('Plane stress & Plane strain analysis(OUTPUT)')
disp('********************************************')
disp(' ')
disp('Node displacement :')
disp(' ')
for j=l:NN
fprintf('U%g\t',j); fprintf('= %G\n',di(3*j-2))
fprintf('V%g\t',j); fprintf('= %G\n',di(3*j-1))
fprintf('W%g\t',j); fprintf('= %G\n',di(3*j))
disp(' ')
end
disp('~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~')
disp('Elements stress')
for n=l:NE
disp(' ')
fprintf(' Element-%g\n',n)
disp(' Sigma(x) Sigma(y) Sigma(z) Tau(xy) Tau(yz) Tau(xz)')
fprintf('Node-%g \t',e2(n)); fprintf('%8.2f\t',Sigma_2(:,:,n))
disp(' ')
fprintf('Node-%g \t',e6(n)); fprintf('%8.3f\t',Sigma_6(:,:,n))
disp(' ')
fprintf('Node-%g \t',e7(n)); fprintf('%8.3f\t',Sigma_7(:,:,n))
disp(' ')
fprintf('Node-%g \t',e3(n)); fprintf('%8.3f\t',Sigma_3(:,:,n))
disp(' ')
fprintf('Node-%g \t',e1(n)); fprintf('%8.3f\t',Sigma_1(:,:,n))
disp(' ')
fprintf('Node-%g \t',e5(n)); fprintf('%8.3f\t',Sigma_5(:,:,n))
disp(' ')
fprintf('Node-%g \t',e8(n)); fprintf('%8.3f\t',Sigma_8(:,:,n))
disp(' ')
fprintf('Node-%g \t',e4(n)); fprintf('%8.3f\t',Sigma_4(:,:,n))
disp(' ')
end
disp('~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~')
disp('Finished!')
disp('Thank you for using of this Program')