program rigidez
implicit none
integer(4) :: i,j,Ne1,Ne
real(8) ::  E, NIU
real(8) :: NODO_r(4), NODO_z(4), K(16,16), Kcopy(16,16), dvgr(4), dvgz(4)
Real::sec1,sec2,seconds

open(3, file='s_anal8_axi.dat')
open(14,file='s_anal8_axi.res')

read(3,*) E, NIU,Ne1  !módulo de Young, coeficiente de Poisson y número de elementos

READ(3,*)(dvgr(i),dvgz(i),i=1,4)

NODO_r(1)=dvgr(1)
NODO_r(2)=dvgr(4)
NODO_r(3)=dvgr(3)
NODO_r(4)=dvgr(2)
NODO_z(1)=dvgz(1)
NODO_z(2)=dvgz(4)
NODO_z(3)=dvgz(3)
NODO_z(4)=dvgz(2)

! do i= 1, 4
!	read(3,*) NODO_r(i), NODO_z(i) ! lectura de los nodos r, z
! end do

call Etime(sec1)
   
  do Ne=1, Ne1
	
	call AXISYMMETRIC(NODO_r, NODO_z,E,NIU,K)
	
  end do
  
  DO i=1,16
    DO j=i,16
      K(i,j)=K(j,i)
    END DO
  END DO
  
  Kcopy=K
  
  K( 3,:)=Kcopy(15,:)
  K( 4,:)=Kcopy(16,:)
  K( 5,:)=Kcopy( 7,:)
  K( 6,:)=Kcopy( 8,:)
  K( 7,:)=Kcopy(13,:)
  K( 8,:)=Kcopy(14,:)
  K( 9,:)=Kcopy( 5,:)
  K(10,:)=Kcopy( 6,:)
  K(13,:)=Kcopy( 3,:)
  K(14,:)=Kcopy( 4,:)
  K(15,:)=Kcopy( 9,:)
  K(16,:)=Kcopy(10,:)
  
  Kcopy=K
  
  K(:, 3)=Kcopy(:,15)
  K(:, 4)=Kcopy(:,16)
  K(:, 5)=Kcopy(:, 7)
  K(:, 6)=Kcopy(:, 8)
  K(:, 7)=Kcopy(:,13)
  K(:, 8)=Kcopy(:,14)
  K(:, 9)=Kcopy(:, 5)
  K(:,10)=Kcopy(:, 6)
  K(:,13)=Kcopy(:, 3)
  K(:,14)=Kcopy(:, 4)
  K(:,15)=Kcopy(:, 9)
  K(:,16)=Kcopy(:,10)
 
call Etime(sec2)

seconds=sec2-sec1

 write(14,*) ' matriz de rigidez '
 write(14,*) ' EL Tiempo con ',Ne1,'elementos es',seconds


 DO i=1,16
  DO j=1,i
    WRITE(14,'(2i3,f20.4)')i,j,K(i,j)
  END DO
END DO 
2000 format (16(' ',e12.4))

stop
CONTAINS 

subroutine AXISYMMETRIC(NODO_r, NODO_z,E,NIU, K)
implicit none
real(8), intent(IN):: E,NIU
real(8),intent(IN) :: NODO_r(4), NODO_z(4)
real(8) :: NOD_r(4), NOD_z(4)
real(8), intent(INOUT) :: K(16,16)
real::A,E1,E2,E3
Real:: r21,r31,r32,r41,r42,r43,Z21,Z31,Z32,Z41,Z42,Z43,rs13,rs24,A1,A2,A3,f1,f2
Real::L1,L2,L3,L4,P1,P2,P3,P5,P6,P7,P8,P10,P11,P12,P13,P14,P15,P16,P17,P18,P19,P20
Real::P21,P22,P23,P24,P25,P26,P27,P28,P29,P30,P32,P34,P35,P36,P37,P38,P40 
Real::P43,P45,P46,P47,P48,P49,P50,P51,P52,P53,P54,P55,P56,P57,P58,P59,P60,P61,P62,P63,P64,P65 
Real::Q1,Q2,Q3,Q4,Q5,Q6,Q7,Q8,Q9,Q11,Q12,Q13,Q14,Q15,Q16,Q17,Q18,Q19,Q20,Q21,Q22,Q23,Q24,Q25 
Real::Q27,Q29,Q30,Q31,Q32,Q33,Q34,Q35,Q36,Q37,Q38,Q39,Q40,Q41,Q42,Q43,Q44,Q45,Q46,Q47,Q48 
Real::Q49,Q50,Q51,Q52,Q53,Q54,Q55,Q56,Q57,Q58,Q59,Q60,Q61,Q62,Q63,Q64,Q65,Q66,Q67,Q68,Q69,Q70,Q71 
Real::Q72,M1,M2,M3,M4,M5,M6,M7,M8,M10,M11,M12,C1,C2,C3,C4,C6,C7,C8,Y1,Y2,Y3,Y4,Y5,Y6,Y7,Y8,Y9,Y10
Real::Y11,Y12,Y13,Y14,E9,E10,E11,E12,J9,J10,J11,J12,J13,J14,J15,J16,J17,J18,J19,J20,J21,J22,J23,J24
Real::U9,U10,U11,U12,U13,U14,U15,U16,U17,U18,U19,U20,S13,S14,S15,S16 
real::D1,D2,D3,D4,D5,D6,D7,D8,x10,x11,x12,E4,E5,E6,E7,E8,J1,J2,J3,J4,J5,J6,J7,J8,U1,U2,U3,U4,U5
real::U6,U7,U8,II1,II2,II3,II4,II5,II6,II7,II8,S1,S2,S3,S4,S5,S6,S7,S8,S9,S10,S11,S12,R1,R2,R3
Real::R4,R5,R6,R7,R8,R9,R10,RL11,RL12,NN1,NN2,NN3,NN4,NN5,NN6,NN7,NN8,PP1,PP2,PP3,PP4,PP5,PP6,PP7
Real::PP8,O1,O2,O3,O4,O5,O6,O7,O8,O9,O10,O11,O12,UU1,UU2,UU3,UU4,UU5,UU6,UU7,UU8,UU9,UU10,UU11
Real::X8,I1,I2,I3,I4,I5,I6,I7,I8,I9,I10,I11,I12,N1,N2,N3,N4,N5,N6,N7,N8,N9,N10,N11,N12
Real::UU12,SL111,SL112,SL113,SL114,SL115,SL116,SL117,SL118,V1,V2,V3,V4,V5,V6,V7,V8,V9,V10,V11,V12,X1,X2,X3,X4,X5,X6,X7
Real::S1A1,S2A1,S3A1,S4A1,S5A1,S6A1,T1A1,T2A1,T3A1,T4A1,T5A1,T6A1
Real::S1A11,S2A11,S3A11,S4A11,S5A11,S6A11,T1A11,T2A11,T3A11,T4A11,T5A11,T6A11
Real::S1A2,S2A2,S3A2,S4A2,S5A2,S6A2,T1A2,T2A2,T3A2,T4A2,T5A2,T6A2
Real::S1A22,S2A22,S3A22,S4A22,S5A22,S6A22,T1A22,T2A22,T3A22,T4A22,T5A22,T6A22
Real::S1B1,S2B1,S3B1,S4B1,S5B1,S6B1,T1B1,T2B1,T3B1,T4B1,T5B1,T6B1
Real::S1B2,S2B2,S3B2,S4B2,S5B2,S6B2,T1B2,T2B2,T3B2,T4B2,T5B2,T6B2
Real::S1C1,S2C1,S3C1,S4C1,S5C1,S6C1,T1C1,T2C1,T3C1,T4C1,T5C1,T6C1
Real::S1C11,S2C11,S3C11,S4C11,S5C11,S6C11,T1C11,T2C11,T3C11,T4C11,T5C11,T6C11
Real::S1C2,S2C2,S3C2,S4C2,S5C2,S6C2,T1C2,T2C2,T3C2,T4C2,T5C2,T6C2
Real::S1C22,S2C22,S3C22,S4C22,S5C22,S6C22,T1C22,T2C22,T3C22,T4C22,T5C22,T6C22
Real::S1D1,S2D1,S3D1,S4D1,S5D1,S6D1,T1D1,T2D1,T3D1,T4D1,T5D1,T6D1
Real::S1D11,S2D11,S3D11,S4D11,S5D11,S6D11,T1D11,T2D11,T3D11,T4D11,T5D11,T6D11
Real::S1D2,S2D2,S3D2,S4D2,S5D2,S6D2,T1D2,T2D2,T3D2,T4D2,T5D2,T6D2
Real::S1D22,S2D22,S3D22,S4D22,S5D22,S6D22,T1D22,T2D22,T3D22,T4D22,T5D22,T6D22
Real::S1E1,S2E1,S3E1,S4E1,S5E1,S6E1,T1E1,T2E1,T3E1,T4E1,T5E1,T6E1
Real::S1E11,S2E11,S3E11,S4E11,S5E11,S6E11,T1E11,T2E11,T3E11,T4E11,T5E11,T6E11
Real::S1E2,S2E2,S3E2,S4E2,S5E2,S6E2,T1E2,T2E2,T3E2,T4E2,T5E2,T6E2
Real::S1E22,S2E22,S3E22,S4E22,S5E22,S6E22,T1E22,T2E22,T3E22,T4E22,T5E22,T6E22
Real::S1F1,S2F1,S3F1,S4F1,S5F1,S6F1,T1F1,T2F1,T3F1,T4F1,T5F1,T6F1
Real::S1F2,S2F2,S3F2,S4F2,S5F2,S6F2,T1F2,T2F2,T3F2,T4F2,T5F2,T6F2
Real::S1G1,S2G1,S3G1,S4G1,S5G1,S6G1,T1G1,T2G1,T3G1,T4G1,T5G1,T6G1
Real::S1G11,S2G11,S3G11,S4G11,S5G11,S6G11,T1G11,T2G11,T3G11,T4G11,T5G11,T6G11
Real::S1G2,S2G2,S3G2,S4G2,S5G2,S6G2,T1G2,T2G2,T3G2,T4G2,T5G2,T6G2
Real::S1G22,S2G22,S3G22,S4G22,S5G22,S6G22,T1G22,T2G22,T3G22,T4G22,T5G22,T6G22
Real::S1H1,S2H1,S3H1,S4H1,S5H1,S6H1,T1H1,T2H1,T3H1,T4H1,T5H1,T6H1
Real::S1H11,S2H11,S3H11,S4H11,S5H11,S6H11,T1H11,T2H11,T3H11,T4H11,T5H11,T6H11
Real::S1H2,S2H2,S3H2,S4H2,S5H2,S6H2,T1H2,T2H2,T3H2,T4H2,T5H2,T6H2
Real::S1H22,S2H22,S3H22,S4H22,S5H22,S6H22,T1H22,T2H22,T3H22,T4H22,T5H22,T6H22
Real::S1I1,S2I1,S3I1,S4I1,S5I1,S6I1,T1I1,T2I1,T3I1,T4I1,T5I1,T6I1
Real::S1I11,S2I11,S3I11,S4I11,S5I11,S6I11,T1I11,T2I11,T3I11,T4I11,T5I11,T6I11
Real::S1I2,S2I2,S3I2,S4I2,S5I2,S6I2,T1I2,T2I2,T3I2,T4I2,T5I2,T6I2
Real::S1I22,S2I22,S3I22,S4I22,S5I22,S6I22,T1I22,T2I22,T3I22,T4I22,T5I22,T6I22 
Real::S1J1,S2J1,S3J1,S4J1,S5J1,S6J1,T1J1,T2J1,T3J1,T4J1,T5J1,T6J1
Real::S1J11,S2J11,S3J11,S4J11,S5J11,S6J11,T1J11,T2J11,T3J11,T4J11,T5J11,T6J11
Real::S1J2,S2J2,S3J2,S4J2,S5J2,S6J2,T1J2,T2J2,T3J2,T4J2,T5J2,T6J2
Real::S1J22,S2J22,S3J22,S4J22,S5J22,S6J22,T1J22,T2J22,T3J22,T4J22,T5J22,T6J22
Real:: H11,H22,f1a,f2a,L6,L7,L8,L9,C9,C10,RL17,RL18,RL19,RL20,RL22,RL23,RL24,RL25,RL26
Real:: RL27,RL28,RL29,RL30,RL31,RL32,RL13,RL14,RL15,RL16,RL21,S17,S18,S19,S20,S21,S22
Real::S23,S24,S25,S26,S27,S28,S29,S30,S31,S32,S33,S34,S35,S36,S37,S38,S39,S40
Real::S41,S42,S43,S44,S45,S46,S47,S48,S49,S50,S51,S52,S53,S54,S55,S56,S57,S58,S59,S60
Real::S61,S62,S63,S64,S65,S66,S67,S68,S69,S70,S71,S72,S73,S74,S75,S76,S77,S78,S79,S80
Real::S81,S82,S83,S84,S85,S86,S87,S88,S89,S90,S91,S92,S93,S94,S95,S96,S97,S98,S99,S100
Real::S101,S102,S103,S104,S105,S106,S107,S108,S109,S110,S111,S112,S113,S114,S115,S116,S117
Real::S118,S119,S120,S121,S122,S123,S124,S125,S126,S127,S128,S129,S130,S131,S132,S133,S134
Real::S135,S136,S137,S138,S139,S140,S141,S142,S143,S144,S145,S146,S147,S148,S149,S150
Real::S151,S152,S153,S154,S155,S156,S157,S158,S159,S160,S161,S162,S163,S164


Real::X,Y,Z,M01,M02,M03,M04,M05,M06,M07,M08,M09,M010,M011,M012,M013,M014,M015
Real:: M016,M017,M018,M019,M020,M021,M022,M023,M024,M025,M026,M027,M028,M029
Real::M030,M031,M032,M033,M034,M035,M036,M037,M038,M039
integer (4)::i



!COEFICIENTES DE LAS EXPRESIONES COMUNES

K=0.0
X=sqrt(3.)
Y=1./3*X
Z=2*Y
M01=1./3+1./2*Y
M02=1./3-1./2*Y
M03=1./6
M04=5./36*X+1./4
M05=1./12-1./36*X
M06=1./36*X+1./12
M07=1./12+5./36*X
M08=-5./36*X+1./4
M09=-1./36*X+1./12
M010=-1./12+5./36*X
M011=4./27-2./9*Y
M012=4./27+2./27*X
M013=2./3-2./3*Y 
M014=2./3+2./9*X
M015=1./9*X
M016=-1./8+1./8*X
M017=1./8+1./8*X
M018=1./8*X-1./8
M019=1./36*X
M020=1./18*X
M021=-1./72*X+1./24
M022=1./72*X+1./24
M023=1./8*(-1+X)
M024=1./8*(X+1)
M025=1./3+1./9*X
M026=1./3-1./9*X
M027=1./108*X
M028=1./24-1./72*X
M029=1./24+1./72*X
M030=2./27
M031=5./36*X+1./4
M032=5./36*X+1./12
M033=1./108
M034=1./36
M035=1./18+1./54*X
M036=1./18-1./54*X
M037=1./54*(1+X)
M038=1./54*(-1+X)
M039=1./9*X


! VALORES EN LOS CUALES QUEDAN EXPRESADAS LAS COMPONENTES DE LA MATRIZ
! DE LAS CONSTANTES ELÁSTICAS   

 A=E/((1+niu)*(1-2*niu))
 E1=A*(1-niu)
 E2=A*niu
 E3=A*((1-2*niu)/2)

!FUNCIONES GENERADORAS DE LOS GRUPOS

r21 = NODO_r(2)-NODO_r(1)
r31 = NODO_r(3)-NODO_r(1)
r32 = NODO_r(3)-NODO_r(2)
r41 = NODO_r(4)-NODO_r(1)
r42 = NODO_r(4)-NODO_r(2)
r43 = NODO_r(4)-NODO_r(3)
Z21 = NODO_z(2)-NODO_z(1)
Z31 = NODO_z(3)-NODO_z(1)
Z32 = NODO_z(3)-NODO_z(2)
Z41 = NODO_z(4)-NODO_z(1)
Z42 = NODO_z(4)-NODO_z(2)
Z43 = NODO_z(4)-NODO_z(3)
rs13= NODO_r(1)+NODO_r(3)
rs24= NODO_r(2)+NODO_r(4)
A1=1./8*(Z41*r32-Z32*r41)
A2=1./8*(Z21*r43-Z43*r21)
A3=1./8*(Z42*r31-Z31*r42)
f1=A1+A2
f2=A1-A2
h11=3*A3*A3-f1*f1
h22=3*A3*A3-f2*f2
 
DO I=1, 4

NOD_r = NODO_r
NOD_z = NODO_z	

	select case(i)
	
	 	case(1) 	 !TRANS P0
			NOD_r = NODO_r
			NOD_z = NODO_z				
		case(2)
							
			nod_r(1) = nodo_r(4)
			nod_r(2) = nodo_r(1)
			nod_r(3) = nodo_r(2)
			nod_r(4) = nodo_r(3)
			nod_z(1) = nodo_z(4)
			nod_z(2) = nodo_z(1)
			nod_z(3) = nodo_z(2)
			nod_z(4) = nodo_z(3)
		case(3)
			nod_r(1) = nodo_r(3)
			nod_r(2) = nodo_r(4)
			nod_r(3) = nodo_r(1)
			nod_r(4) = nodo_r(2)
			nod_z(1) = nodo_z(3)
			nod_z(2) = nodo_z(4)
			nod_z(3) = nodo_z(1)
			nod_z(4) = nodo_z(2)
		case(4)
		
			nod_r(1) = nodo_r(2)
			nod_r(2) = nodo_r(3)
			nod_r(3) = nodo_r(4)
			nod_r(4) = nodo_r(1)
			nod_z(1) = nodo_z(2)
			nod_z(2) = nodo_z(3)
			nod_z(3) = nodo_z(4)
			nod_z(4) = nodo_z(1)
					
		endselect

r21 = nod_r(2)-nod_r(1)
r31 = nod_r(3)-nod_r(1)
r32 = nod_r(3)-nod_r(2)
r41 = nod_r(4)-nod_r(1)
r42 = nod_r(4)-nod_r(2)
r43 = nod_r(4)-nod_r(3)
Z21 = nod_z(2)-nod_z(1)
Z31 = nod_z(3)-nod_z(1)
Z32 = nod_z(3)-nod_z(2)
Z41 = nod_z(4)-nod_z(1)
Z42 = nod_z(4)-nod_z(2)
Z43 = nod_z(4)-nod_z(3)
rs13= nod_r(1)+nod_r(3)
rs24= nod_r(2)+nod_r(4)
A1=1./8*(Z41*r32-Z32*r41)
A2=1./8*(Z21*r43-Z43*r21)
f1=A1+A2
f2=A1-A2 

!ECUACIONES COMUNES

L1=Y*f1+A3 
L2=-Y*f1+A3
L3=Y*f2+A3
L4=-Y*f2+A3 
L6=L1*L1
L7=L2*L2
L8=L3*L3
L9=L4*L4
P1=M019*nod_z(1)-M020*nod_z(3)
P2=M019*nod_r(1)-M020*nod_r(3) 
P3=M023*z42
P5=M024*z42
P6=M03*rs24
P7=M03*rs13
P8=M023*r42
P10=M024*r42
P11=M08*nod_z(1)
P12=M06*nod_z(2)
P13=M010*nod_z(4)
P14=M06*nod_z(3)
P15=M05*nod_z(2)
P16=M010*nod_z(3)
P17=M06*nod_z(1)
P18=M08*nod_z(2)
P19=M05*nod_z(1)
P20=M06*nod_z(4)
P21=M010*nod_r(3)
P22=M08*nod_r(2)
P23=M06*nod_r(1)
P24=M06*nod_r(4)
P25=M05*nod_r(1)
P26=M06*nod_r(2)
P27=M08*nod_r(1)
P28=M06*nod_r(3)
P29=M05*nod_r(2)
P30=M023*z31
P32=M024*r31
P34=M023*r31
P35=M024*z31
P36=M019*nod_r(2)-M020*nod_r(4)
P37=M019*nod_z(2)-M020*nod_z(4)
P38=M019*nod_z(3)-M020*nod_z(1)
P40=M019*nod_r(3)-M020*nod_r(1)
P43=M019*nod_z(4)-M020*nod_z(2)
P45=M019*nod_r(4)-M020*nod_r(2)
P46=M05*nod_r(3)
P47=M05*nod_r(4)
P48=M010*nod_z(1)
P49=M010*nod_r(1)
P50=M05*nod_z(4)
P51=M05*nod_z(3)
P52=M08*nod_r(3)
P53=M010*nod_r(4)
P54=M08*nod_z(3)
P55=M08*nod_r(4)
P56=M010*nod_r(2)
P57=M010*nod_z(2)
Q1=P6+M02*nod_r(1)+M01*nod_r(3)
Q2=P3
Q4=-P5
Q3=P6+M01*nod_r(1)+M02*nod_r(3)
Q5=-P1-M022*nod_z(2)+M021*nod_z(4)
Q6=P1-M021*nod_z(2)+M022*nod_z(4)
Q7=P7+M02*nod_r(2)+M01*nod_r(4)
Q8=P7+M01*nod_r(2)+M02*nod_r(4)
Q9=-P2+M021*nod_r(2)-M022*nod_r(4)
Q11=P2+M022*nod_r(2)-M021*nod_r(4)
Q12=-P8
Q13=P10
Q14=M09*nod_z(3)-P11+P12-P13
Q15=P14+P15-M04*nod_z(1)+M07*nod_z(4) 
Q16=P16-P17+P18-M09*nod_z(4) 
Q17=-M07*nod_z(3)+M04*nod_z(2)-P19-P20 
Q18=-P21-P22+P23+M09*nod_r(4)
Q19=P24+P25+M07*nod_r(3)-M04*nod_r(2) 
Q20=-P26+P27-M09*nod_r(3)+P53
Q21=-P28-P29-M07*nod_r(4)+M04*nod_r(1)
Q22=-P30
Q23=-P32 
Q24=P34 
Q25=P35  
Q27=-P36-M022*nod_r(1)+M021*nod_r(3)
Q29=P36-M028*nod_r(1)+M029*nod_r(3) 
Q30=P37+M022*nod_z(1)-M021*nod_z(3)
Q31=-P37+M028*nod_z(1)-M029*nod_z(3) 
Q32=P38+M029*nod_z(2)-M028*nod_z(4)  
Q33=P30  
Q34=-P38+M021*nod_z(2)-M022*nod_z(4)
Q35=-P35 
Q36=-P40-M029*nod_r(2)+M028*nod_r(4) 
Q37=-P34
Q38=P40-M021*nod_r(2)+M022*nod_r(4)
Q39=P32 
Q40=-P3
Q41=-P43-M029*nod_z(1)+M028*nod_z(3)
Q42=P5 
Q43=P43-M021*nod_z(1)+M022*nod_z(3)
Q44=P8 
Q45=P45+M029*nod_r(1)-M028*nod_r(3)
Q46=-P10 
Q47=-P45+M021*nod_r(1)-M022*nod_r(3)  
Q48=-P24-P46+M031*nod_r(2)-M032*nod_r(1)
Q49=-P28+P22-P47+P49
Q50=P14-P18-P48+P50
Q51=P51+M032*nod_z(1)-M031*nod_z(2)+P20
Q52=M031*nod_z(3)-P17-P15-M032*nod_z(4) 
Q53=-P52-P53+P26+P25
Q54=-M031*nod_r(3)+M032*nod_r(4)+P29+P23
Q55=P54-P19-P12+P13 
Q56=-P14+P48-M09*nod_z(2)+M08*nod_z(4)
Q57=-P16+P17+P15-M08*nod_z(4)
Q58=P21+P55-P29-P23
Q59=-P51-M07*nod_z(1)-P12+M04*nod_z(4)
Q60=-M032*nod_r(3)+M031*nod_r(4)-P26-P25
Q62=M032*nod_z(3)+P19+P12-M031*nod_z(4)
Q61=P28-P55+M09*nod_r(2)-P49
Q63=P46-M04*nod_r(4)+P26+M07*nod_r(1)
Q64=-P54+M09*nod_z(1)-P57+P20
Q65=P28+P47+M032*nod_r(2)-M031*nod_r(1) 
Q66=-M04*nod_z(3)+P17+M07*nod_z(2)+P50
Q67=P46+P24-P56-P27 
Q68=-P14+M031*nod_z(1)-M032*nod_z(2)-P50
Q69=P52-P24+P56-M09*nod_r(1)
Q70=-P51+P11+P57-P20 
Q71=M04*nod_r(3)-P47-M07*nod_r(2)-P23
Q72=Q4*Q40 
P58=L6/Q1
P59=L7/Q3
P60=L8/Q7
P61=Q2*L1
P62=Q4*L2
P63=Q5*L4
P64=Q6*L3
P65=Q7*Q9*Q9
M1=M033*P58+Q1*Q2*Q2
M2=M033*P59+Q3*Q4*Q4 
M3=M015*P61
M4=M039*P62
M5=Q1*Q12*Q12
M6=Q3*Q13*Q13
M7=(1./3)*P63
M8=(1./3)*P64
M10=Q8*Q11**2
M11=(1./36)*P60+Q7*Q6*Q6
M12=(1./36)*L9/Q8+Q8*Q5*Q5
c1=Q3*M017*M017 
c2=Q1*M016*M016 
C3=Q1*M018*M018 
C4=C1
C6=Q8*Q11*Q11
C7=Q7*Q6*Q6
C8=Q8*Q5*Q5
C9=r42*r42
C10=z42*z42
Y1=M013*Q14*L1
Y2=M014*Q15*L2 
Y3=Q1*Q20*Q20 
Y4=Q3*Q21*Q21
Y5=M013*Q16*L3 
Y6=M014*Q17*L4
Y7=Q8*Q19*Q19
Y8=Q7*Q18*Q18
Y9=Q1*Q14*Q14 
Y10=Q3*Q15*Q15
Y11=M012*P59+Y10
Y12=M011*P58+Y9
Y13=M011*P60+Q7*Q16*Q16
Y14=M012*L9/Q8+Q8*Q17*Q17
D1=Q1*Q20*Q20
D2=Q3*Q21*Q21 
D3=Q1*Q14*Q14 
D4=Q3*Q15*Q15
D5=Q7*Q18*Q18
D6=Q8*Q19*Q19
D7=Q7*Q16*Q16
D8=Q8*Q17*Q17 
x10=Q1*Q12*Q2
x11=Q3*Q13*Q4 
x12=M020*Q13*L2
E4=M020*Q12*L1
E5=Q8*Q11*Q5
E6=Q7*Q9*Q6
E7=M03*Q11*L4
E8=M03*Q9*L3
E9=x11+x12
E10=x10-E4
E11=E5-E7
E12=E6-E8
J1=Q1*Q20*Q14
J2=M026*Q20*L1
J3=Q3*Q21*Q15 
J4=M025*Q21*L2 
J5=Q7*Q18*Q16
J6=Q8*Q19*Q17 
J7=M026*Q18*L3
J8=M025*Q19*L4
J9=J1+J2
J10=J3+J4
J11=J6+J8
J12=J5+J7
J13=M027*P58+Q1*Q2* Q31
J14=-M027*P59+Q3*Q4*Q30
J15=Q1*Q29*Q12
J16=Q3*Q27*Q13
J17=M03*P61+M020*L1*Q31
J18=-M03*P62+M020*L2*Q30
J19=M020*P63-M03*L4*Q25
J20=M020*P64+M03*L3*Q22
J21=Q8*Q23*Q11
J22=Q7*Q24*Q9
J23=-M027*L9/Q8+Q8*Q5*Q25
J24=M027*P60+Q7*Q6*Q22
U1=Q1*Q47*Q46
U2=Q3*Q45*Q44
U3=Q1*Q43*Q42
U4=Q3*Q41*Q40
U5=Q7*Q39*Q38
U6=Q8*Q37*Q36
U7=Q7*Q35*Q34
U8=Q8*Q33*Q32
U9=M030*P58+Q1*Q14 *Q52
U10=M030*P59+Q3*Q15*Q55
U11=M026*Q15*L2+M025*L2*Q55
U12=M025*Q14*L1+M026*L1*Q52
U13=Q1*Q54*Q20
U14=Q3*Q53*Q21
U15=M011*P60+Q7*Q16 *Q50
U16=M012*L9/Q8+Q8*Q17*Q51
U17=Q7*Q49*Q18
U18=Q8*Q48*Q19
U19=M026*Q16*L3+M026*L3*Q50
U20=M025*Q17*L4+M025*L4*Q51
II1=Q1*Q54*Q20
II2=Q3*Q53*Q21 
II3=Q1*Q52*Q14
II4=Q3*Q55*Q15 
II5=Q7*Q49*Q18 
II6=Q8*Q48*Q19
II7=Q7*Q50*Q16
II8=Q8*Q51*Q17 
S1=Q1*Q29*Q2
S2=Q1*Q31*Q12   
S3=Q3*Q30*Q13 
S4=M020*Q29*L1
S5=Q3*Q27*Q4  
S6=M020*Q27*L2
S7=Q7*Q24*Q6  
S8=M03*Q24*L3
S9=Q8*Q23*Q5 
S10=M03*Q23*L4 
S11=Q8*Q25*Q11  
S12=Q7*Q22*Q9
S13=S5+S6
S14=S1-S4
S15=S7-S8
S16=S9-S10
R1=Q1*Q52*Q20
R2=Q3*Q55*Q21
R3=Q1*Q54*Q14 
R4=M026*Q54*L1  
R5=Q3*Q53*Q15
R6=M025*Q53*L2
R7=Q7*Q49*Q16
R8=M026*Q49*L3
R9=Q8*Q48*Q17 
R10=M025*Q48*L4 
RL11=Q7*Q50*Q18
RL12=Q8*Q51*Q19
RL13=R5+R6
RL14=R3+R4
RL15=R9+R10
RL16=R7+R8
RL17=Q1*Q47*Q2-M020*Q47 *L1
RL18=Q3*Q45*Q4+M020*Q45*L2
RL19=Q1*Q43*Q12
RL20=Q3*Q41*Q13
RL21=Q7*Q39*Q6-M03*Q39*L3
RL22=Q8*Q37*Q5-M03*Q37*L4
RL23=Q7*Q35*Q9
RL24=Q8*Q33*Q11 
RL25=Q1*Q71*Q70+M026*Q71*L1
RL26=Q3*Q69*Q68+M025*Q69*L2
RL27=Q3*Q65*Q64
RL28=Q1*Q67*Q66
RL29=Q7*Q63*Q62+M025*Q63*L3
RL30=Q8*Q61*Q57+M026*Q61*L4
RL31=Q8*Q58*Q56
RL32=Q7*Q60*Q59
S17=M020*P61-M020*L1*Q42
S18=-M020*P62+M020*L2*Q40
S19=Q1*Q46*Q12
S20=Q3*Q44*Q13
S21=M033*P58-1.*Q1*Q2*Q42
S22=-M033*P59+Q3 *Q72
S23=Q8*Q36*Q11
S24=Q7*Q38*Q9
S25=M034*L9/Q8+Q8*Q5*Q32
S26=M034*P60+Q7*Q6*Q34
S27=-M03*P63-M03*L4*Q32
S28=M03*P64+M03*L3 *Q34
nn1=Q1*Q47*Q29 
nn2=Q3*Q45*Q27
nn3=Q1*Q43*Q31
nn4=Q3*Q41*Q30
nn5=Q7*Q39*Q24
nn6=Q8*Q37*Q23
nn7=Q7*Q35*Q22 
nn8=Q8*Q33*Q25 
S29=M030*P58+Q1*Q14*Q66
S30=M030*P59+Q3 *Q15*Q64
S31=Q1*Q71*Q20
S32=Q3*Q69*Q21
S33=M025*Q14*L1+M026*L1*Q66
S34=M026*Q15*L2+M025*L2*Q64
S35=M030*P60+Q7*Q16*Q59
S36=M030*L9/Q8 +Q8*Q17*Q56
S37=Q7*Q63*Q18
S38=Q8*Q61*Q19
S39=M025*Q16*L3+M026*L3*Q59
S40=M026*Q17*L4+M025*L4 *Q56
PP1=Q1*Q67*Q54
PP2=Q3*Q65*Q53 
PP3=Q1*Q70*Q52
PP4=Q3*Q68*Q55 
PP5=Q7*Q60*Q49
PP6=Q8*Q58*Q48
PP7=Q7*Q62*Q50
PP8=Q8*Q57*Q51
O1=Q1*Q46*Q2
O2=M020*Q46*L1
O3=Q3*Q44*Q4
O4=M020*Q44*L2
O5=-Q1*Q42*Q44 
O6=Q3*Q40*Q13 
O7=Q7*Q38*Q6
O8=M03*Q38*L3
O9=Q8*Q36*Q5
O10=M03*Q36* L4 
O11=Q7*Q34*Q9
O12=Q8*Q32*Q11
S41=O1-O2
S42=O3+O4
S43=O9-O10
S44=O7-O8
uu1=Q1*Q71*Q14
uu2=M026*Q71 *L1
uu3=Q3*Q69*Q15
uu4=M025*Q69*L2
uu5=Q1*Q66*Q20 
uu6=Q3*Q64*Q21
uu7=Q7*Q63*Q16
uu8=M026*Q63*L3
uu9=Q8*Q61*Q17 
uu10=M025*Q61*L4
uu11=Q7*Q59*Q18
uu12=Q8*Q56*Q19 
S45=uu1+uu2
S46=uu3+uu4
S47=uu7+uu8
S48=uu9+uu10
S49=M035*P58-1.*Q1*Q43*Q66
S50=-M036*P59+Q3*Q41 *Q64
S51=M026*Q41*L2-M03*L2*Q64
S52=M025*Q43*L1-M03*L1*Q66
S53=Q1*Q71*Q47
S54=Q3*Q69*Q45
S55=M037*P60+Q7*Q35*Q59
S56=-M038*L9/Q8+Q8*Q33 *Q56
S57=M026*Q33*L4-M020*L4*Q56
S58=M025*Q35*L3+M020*L3*Q59
S59=Q7*Q63*Q39
S60=Q8*Q61*Q37
S61=-M035*P59+Q3*Q41*Q68
S62=M036*P58-1.*Q1*Q43*Q70
S63=Q3*Q65*Q45
S64=Q1*Q67*Q47
S65=M026*Q43*L1-M03*L1*Q70
S66=M025*Q41*L2-M03*L2*Q68
S67=M037*P60+Q7*Q35*Q62
S68=-M038*L9/Q8+Q8*Q33*Q57
S69=M025*Q35*L3+M020*L3*Q62
S70=M026*Q33 *L4-M020*L4*Q57
S71=Q7*Q60*Q39
S72=Q8*Q58*Q37
SL111=Q1*Q20*Q12
SL112=Q3*Q21*Q13
SL113=Q1*Q14*Q2
SL114=Q3*Q15*Q4 
SL115=Q7*Q18*Q9 
SL116=Q8*Q19*Q11
SL117=Q7*Q16*Q6
SL118=Q8*Q17*Q5
S73=Q1*Q54*Q46
S74=Q3*Q53*Q44
S75=Q1*Q52*Q42
S76=Q3*Q55*Q40
S77=Q7*Q49*Q38
S78=Q8*Q48*Q36
S79=Q8*Q51*Q32
S80=Q7*Q50*Q34
S81=Q1*Q46*Q52+M025*Q46*L1
S82=Q3*Q44*Q55+M026*Q44*L2
S83=Q1*Q54*Q42
S84=Q3*Q53*Q40
S85=Q7*Q38*Q50+M026*Q38 *L3
S86=Q8*Q36*Q51+M025*Q36 *L4
S87=Q7*Q49*Q34
S88=Q8*Q48*Q32 
V1=Q1*Q12*Q14
V2=M026*Q12*L1 
V3=Q3*Q13*Q15
V4=M025*Q13*L2
V5=Q1*Q20*Q2 
V6=Q3*Q21*Q4
V7=Q7*Q9*Q16
V8=M026*Q9*L3
V9=Q8*Q11*Q17
V10=M025*Q11*L4 
V11=Q7*Q18*Q6
V12=Q8*Q19*Q5
S89=V1+V2
S90=V3+V4
S91=V7+V8
S92=V9+V10
S93=Q1*Q20*Q2-M020*Q20 *L1
S94=Q3*Q21*Q4+M020*Q21*L2
S95=Q1*Q14*Q12
S96=Q3*Q15*Q13
S97=Q7*Q18*Q6-M03*Q18*L3
S98=Q8*Q19 *Q5-M03*Q19 *L4
S99=Q7*Q16*Q9
S100=Q8*Q17*Q11
S101=Q1*Q54*Q42+M020*Q54*L1
S102=Q3*Q53*Q40-M020*Q53*L2
S103=Q3*Q55*Q44
S104=Q1*Q52*Q46
S105=Q7*Q49*Q34-M03*Q49 *L3
S106=Q8*Q48*Q32-M03*Q48 *L4
S107=Q8*Q51*Q36
S108=Q7*Q50*Q38
S109=M038*P59+Q3*Q4*Q55
S110=M037*P58-1.*Q1*Q2*Q52
S111=Q3*Q53*Q13
S112=Q1*Q54 *Q12
S113=M025*P61-M020*L1*Q52
S114=M026*P62+M020*L2*Q55
S115=Q7*Q49*Q9
S116=Q8*Q48 *Q11
S117=M025*P63-M03*L4*Q51
S118=M026*P64-M03*L3*Q50
S119=-M035*L9/Q8+Q8*Q5*Q51
S120=M036*P60-1.*Q7*Q6*Q50
S121=Q8*Q19*Q36
S122=Q7*Q18*Q38
S123=M025*Q32*L4-M03*L4*Q17
S124=M026*Q34*L3-M03*L3*Q16
S125=-M035*L9/Q8+Q8 *Q32*Q17
S126=M036*P60-1.*Q7*Q34*Q16
S127=Q3*Q21*Q44
S128=Q1*Q20*Q46 
S129=M025*Q40*L2-M020*L2*Q15
S130=M026*Q42*L1+M020*L1*Q14
S131=M038*P58+Q1*Q42*Q14
S132=-M037*P59+Q3*Q40 *Q15
X1=Q1*Q54*Q12
X2=Q3*Q53*Q13
X3=Q1*Q52*Q2 
X4=Q3*Q55*Q4
X5=Q7*Q49*Q9 
X6=Q8*Q48*Q11
X7=Q7*Q50*Q6
X8=Q8*Q51*Q5   
S133=Q1*Q20*Q46
S134=Q3*Q21*Q44
S135=Q1*Q14*Q42
S136=Q3*Q15*Q40
S137=Q7*Q18*Q38
S138=Q8*Q19*Q36
S139=Q7*Q16*Q34
S140=Q8*Q17*Q32
I1=Q1*Q12*Q52
I2=M025*Q12*L1 
I3=Q3*Q13*Q55 
I4=M026*Q13*L2
I5=Q1*Q54*Q2
I6=Q3*Q53*Q4 
I7=Q7*Q9*Q50
I8=M026*Q9*L3
I9=Q8*Q11*Q51
I10=M025*Q11*L4 
I11=Q7*Q49*Q6 
I12=Q8*Q48*Q5
S141=I1+I2
S142=I3+I4
S143=I7+I8
S144=I9+I10
N1=Q1*Q54*Q2
N2=M020*Q54*L1
N3=Q3*Q53*Q4 
N4=M020*Q53*L2
N5=Q1*Q52*Q12
N6=Q3*Q55*Q13
N7=Q7*Q49*Q6
N8=M03*Q49*L3
N9=Q8*Q48*Q5
N10=M03*Q48*L4
N11=Q7*Q50*Q9 
N12=Q8*Q51*Q11
S145=N1-N2
S146=N3+N4
S147=N7-N8
S148=N9-N10
S149=Q1*Q12*Q66+M025*Q12*L1
S150=Q3*Q13*Q64+M026*Q13*L2
S151=Q1*Q71*Q2
S152=Q3*Q69*Q4
S153=Q7*Q9*Q59+M025*Q9 *L3
S154=Q8*Q11*Q56+M026*Q11*L4
S155=Q7*Q63*Q6
S156=Q8*Q61*Q5
S157=Q1*Q71*Q2-M020*Q71*L1
S158=Q3*Q69*Q4+M020*Q69*L2
S159=Q1*Q66*Q12
S160=Q3*Q64*Q13
S161=Q7*Q63*Q6-M03*Q63 *L3
S162=Q8*Q61*Q5-M03*Q61*L4
S163=Q7*Q59*Q9
S164=Q8*Q56*Q11 

		 
!Grupo A (FUNCIONES GENERADORAS DE  K(1,1))  

s1A1=M1+M2
s2A1=-M3+M4
s3A1=M5+M6
s4A1=M2-M1 
s5A1=M4+M3 
s6A1=M6-M5 
t1A1=M11+M12
t2A1=-M8-M7  
t3A1=P65+M10 
t4A1=M12-M11 
t5A1=-M7+M8 
t6A1=M10-P65     

!Grupo A	(FUNCIONES GENERADORAS DE  K(2,2))


s1A11=(c1+c2)*C9
s2A11= 0
s3A11=(C3+C4)*C10
s4A11=(c1-c2)*C9
s5A11=0  
s6A11=(C4-C3)*C10 
t1A11=P65+C6
t2A11 =0 
t3A11=C7+C8
t4A11=C6-P65 
t5A11 =0 
t6A11=C8-C7           

!Grupo A (FUNCIONES GENERADORAS DE  K(9,9))


s1A2=Y12 +Y11 
s2A2=Y1+Y2
s3A2=Y3+Y4 
s4A2=Y11-Y12 
s5A2=Y2-Y1 
s6A2=Y4-Y3 
t1A2=Y13+Y14 
t2A2=Y5+Y6 
t3A2=Y8+Y7
t4A2=Y14-Y13 
t5A2=Y6-Y5
t6A2=Y7-Y8    

!Grupo A2 (FUNCIONES GENERADORAS DE  K(10,10))

s1A22=D1+D2
s2A22= 0
s3A22=D3+D4
s4A22=D2-D1
s5A22=0 
s6A22=D4-D3
t1A22=D5+D6
t2A22 =0
t3A22=D7+D8 
t4A22=D6-D5
t5A22 = 0
t6A22=D8-D7       

!Grupo B (FUNCIONES GENERADORAS DE  K(1,2))
 


s1B1=0
s2B1=E9+E10 
s3B1=x10+x11 
s4B1=0 
s5B1=E9-E10 
s6B1=x11-x10   
t1B1 =0
t2B1=E11+E12
t3B1=E6+E5  
t4B1 =0 
t5B1=E11-E12
t6B1=E5-E6 
  
!Grupo B (FUNCIONES GENERADORAS DE  K(9,10))


s1B2=0
s2B2=J10+J9
s3B2=J1+J3 
s4B2=0 
s5B2=J10-J9
s6B2=J3-J1 
t1B2=0 
t2B2=J11+J12
t3B2=J5+J6 
t4B2=0
t5B2 =J11-J12 
t6B2=J6-J5           

!Grupo C (FUNCIONES GENERADORAS DE  K(3,1))


s1C1=J13+J14
s2C1=-J17+J18 
s3C1=J15+J16 
s4C1=J14-J13
s5C1=J18+J17 
s6C1=J16-J15  
t1C1=J24+J23
t2C1=-J20+J19  
t3C1=J22+J21 
t4C1=J23-J24 
t5C1=J19+J20 
t6C1=J21-J22
      

!Grupo C (FUNCIONES GENERADORAS DE  K(6,8))



s1C11=U1+U2
s2C11=0 
s3C11=U3+U4
s4C11=U2-U1
s5C11=0 
s6C11=U4-U3
t1C11=U5+U6 
t2C11 =0
t3C11=U7+U8 
t4C11=U6-U5 
t5C11=0 
t6C11=U8-U7    


!Grupo C (FUNCIONES GENERADORAS DE  K(9,11))




s1C2=U9+U10
s2C2=U12+U11 
s3C2=U13+U14
s4C2=U10-U9
s5C2=U11-U12
s6C2=U14-U13  
t1C2=U15+U16
t2C2=U19+U20 
t3C2=U17+U18 
t4C2=U16-U15
t5C2=U20-U19 
t6C2=U18-U17

  

!Grupo C (FUNCIONES GENERADORAS DE  K(10,12))


s1C22=II1+II2 
s2C22=0 
s3C22=II3+II4 
s4C22=II2-II1
s5C22=0 
s6C22=II4-II3
t1C22=II5+II6 
t2C22=0 
t3C22=II7+II8
t4C22=II6-II5
t5C22=0
t6C22=II8-II7    

!grupo D


       
s1D1=0 
s2D1=S14+S13
s3D1=S2+S3
s4D1=0
s5D1=S13-S14
s6D1=S3-S2  
t1D1=0  
t2D1=S16+S15
t3D1=S12+S11 
t4D1=0   
t5D1=S16-S15
t6D1=S11-S12        

!grupo D

s1D2=0
s2D2=RL14+RL13
s3D2=R1+R2
s4D2=0 
s5D2=RL13-RL14
s6D2=R2-R1 
t1D2 =0   
t2D2=RL16+RL15 
t3D2=RL11+RL12
t4D2=0   
t5D2=RL15-RL16
t6D2=RL12-RL11        
 
!grupo D


s1D11=0
s2D11=RL17+RL18
s3D11=RL19+RL20
s4D11=0
s5D11=RL18-RL17 
s6D11=RL20-RL19
t1D11=0 
t2D11=RL21+RL22
t3D11=RL23+RL24
t4D11=0 
t5D11=RL22-RL21
t6D11=RL24-RL23 

!grupo D


s1D22=0
s2D22=RL25+RL26
s3D22=RL28+RL27
s4D22=0
s5D22=RL26-RL25 
s6D22=RL27-RL28 
t1D22=0 
t2D22=RL29+RL30 
t3D22=RL32+RL31 
t4D22=0 
t5D22=RL30-RL29
t6D22=RL31-RL32 

!Grupo E1


s1E1=-S21+S22
s2E1=S17+S18
s3E1=S19+S20 
s4E1=S22+S21
s5E1=S18-S17
s6E1=S20-S19  
t1E1=S26+S25
t2E1=-S28+S27
t3E1=S24+S23
t4E1=S25-S26 
t5E1=S27+S28
t6E1=S23-S24  

!Grupo E11

s1E11=nn1+nn2
s2E11=0 
s3E11=nn3+nn4
s4E11=nn2-nn1
s5E11=0
s6E11=nn4-nn3
t1E11=nn5+nn6
t2E11=0 
t3E11=nn7+nn8
t4E11=nn6-nn5 
t5E11=0 
t6E11=nn8-nn7  

!Grupo E2
 
 s1E2=S29+S30 
s2E2=S33+S34 
s3E2=S31+S32
s4E2=S30-S29
s5E2=S34-S33 
s6E2=S32-S31
t1E2 =S35+S36
t2E2 =S39+S40 
t3E2 =S37+S38 
t4E2 =S36-S35
t5E2=S40-S39
t6E2=S38-S37 
  
!Grupo E22


s1E22=PP1+PP2
s2E22=0
s3E22=PP3+PP4
s4E22=PP2-PP1
s5E22=0
s6E22=PP4-PP3
t1E22=PP5+PP6
t2E22=0
t3E22=PP7+PP8
t4E22=PP6-PP5
t5E22=0
t6E22=PP8-PP7    

!Grupo F1

 
s1F1=0 
s2F1=S41+S42
s3F1=O5+O6
s4F1=0
s5F1=S42-S41
s6F1=O6-O5
t1F1=0 
t2F1=S44+S43
t3F1=O11+O12 
t4F1=0
t5F1=S43-S44
t6F1=O12-O11  

!Grupo F2


 s1F2=0
 s2F2=S45+S46 
 s3F2=uu5+uu6
 s4F2=0 
 s5F2=S46-S45
 s6F2=uu6-uu5
 t1F2=0
 t2F2=S47+S48
 t3F2=uu11+uu12
 t4F2=0
 t5F2=S48-S47
 t6F2=uu12-uu11 
 
 !GRUPOG

s1G2=SL111+SL112
s2G2=0 
s3G2=SL113+SL114 
s4G2=SL112-SL111
s5G2=0
s6G2=SL114-SL113 
t1G2=SL115+SL116
t2G2=0 
t3G2=SL117+SL118
t4G2=SL116-SL115
t5G2=0 
t6G2=SL118-SL117      

!grupo G11


 
s1G11=-S49+S50
s2G11=S52+S51 
s3G11=S53+S54 
s4G11=S50+S49
s5G11=S51-S52
s6G11=S54-S53 
t1G11=S55+S56
t2G11=S58+S57
t3G11=S59+S60
t4G11=S56-S55 
t5G11=S57-S58 
t6G11=S60-S59  

!grupo G1

s1G1=-S62+S61
s2G1=S65+S66 
s3G1=S64+S63
s4G1=S61+S62
s5G1=S66-S65 
s6G1=S63-S64
t1G1=S67+S68
t2G1=S69+S70
t3G1=S71+S72 
t4G1=S68-S67
t5G1=S70-S69 
t6G1=S72-S71 
 
!grupo G2
  
s1G22=S73+S74
s2G22=0
s3G22=S75+S76
s4G22=S74-S73
s5G22=0
s6G22=S76-S75
t1G22=S77+S78
t2G22=0
t3G22=S80+S79
t4G22=S78-S77
t5G22=0 
t6G22=S79-S80 

!grupo H11
    
s1H11=0 
s2H11=S81+S82
s3H11=S83+S84
s4H11=0
s5H11=S82-S81
s6H11=S84-S83 
t1H11=0
t2H11=S85+S86
t3H11=S87+S88
t4H11=0
t5H11=S86-S85
t6H11=S88-S87  

!grupo H1

s1H1=0
s2H1=S89+S90
s3H1=V5+V6
s4H1=0
s5H1=S90-S89
s6H1=V6-V5
t1H1=0
t2H1=S91+S92
t3H1=V11+V12
t4H1=0
t5H1=S92-S91
t6H1=V12-V11        
 
!grupo H2


s1H2=0
s2H2=S93+S94
s3H2=S95+S96
s4H2=0
s5H2=S94-S93
s6H2=S96-S95
t1H2 =0 
t2H2=S97+S98
t3H2=S99+S100
t4H2=0 
t5H2=S98-S97
t6H2=S100-S99

!grupo H22


s1H22=0 
s2H22=S101+S102
s3H22=S104+S103 
s4H22=0 
s5H22=S102-S101 
s6H22=S103-S104
t1H22=0 
t2H22=S105+S106
t3H22=S108+S107
t4H22=0
t5H22=S106-S105
t6H22=S107-S108

!grupo I1


s1I1=-S110+S109
s2I1=S113+S114
s3I1=S112+S111
s4I1=S109+S110
s5I1=S114-S113
s6I1=S111-S112
t1I1=-S120+S119
t2I1=S118+S117
t3I1=S115+S116
t4I1=S119+S120
t5I1=S117-S118
t6I1=S116-S115  

!grupo I11


s1I11=S131+S132
s2I11=S130+S129
s3I11=S128+S127 
s4I11=S132-S131
s5I11=S129-S130
s6I11=S127-S128
t1I11=-S126+S125
t2I11=S124+S123
t3I11=S122+S121
t4I11=S125+S126
t5I11=S123-S124
t6I11=S121-S122
 
!grupo I2


s1I2=X1+X2
s2I2=0 
s3I2=X3+X4
s4I2=X2-X1 
s5I2=0
s6I2=X4-X3
t1I2=X5+X6 
t2I2=0 
t3I2=X7+X8
t4I2=X6-X5
t5I2=0
t6I2=X8-X7  

!grupo I22

s1I22=S133+S134
s2I22=0 
s3I22=S135+S136
s4I22=S134-S133 
s5I22=0 
s6I22=S136-S135
t1I22=S137+S138
t2I22=0 
t3I22=S139+S140 
t4I22=S138-S137 
t5I22=0 
t6I22=S140-S139

!grupo J1
 

s1J1=0 
s2J1=S141+S142
s3J1=I5+I6 
s4J1=0
s5J1=S142-S141
s6J1=I6-I5
t1J1=0
t2J1=S143+S144
t3J1=I11+I12
t4J1=0 
t5J1=S144-S143
t6J1=I12-I11       

!grupo J2

 
s1J2=0 
s2J2=S145+S146
s3J2=N5+N6 
s4J2=0 
s5J2=S146-S145
s6J2=N6-N5  
t1J2=0
t2J2=S147+S148
t3J2=N11+N12
t4J2=0 
t5J2=S148-S147
t6J2=N12-N11      
 
!grupo J11

s1J11=0 
s2J11=S149+S150 
s3J11=S151+S152
s4J11=0
s5J11=S150-S149 
s6J11=S152-S151
t1J11=0
t2J11=S153+S154
t3J11=S155+S156 
t4J11=0
t5J11=S154-S153
t6J11=S156-S155

!grupo J22

s1J22=0
s2J22=S157+S158
s3J22=S159+S160
s4J22=0
s5J22=S158-S157
s6J22=S160-S159
t1J22=0
t2J22=S161+S162
t3J22=S163+S164
t4J22=0
t5J22=S162-S161
t6J22=S164-S163

    
	select case(i)
		case(1)
		F1A= F1
		F2A= F2
  			
 K(1,1)=((3*A3*(E1*s1A1+E2*s2A1+E3*s3A1)+X*f1a*(E1*s4A1+E2*s5A1+E3*s6A1))/(h11))&
        &+((3*A3*(E1*t1A1+E2*t2A1+E3*t3A1)+X*f2a*(E1*t4A1+E2*t5A1+E3*t6A1))/(h22))
				
 K(2,2)=((3*A3*(E1*s1A11+E2*s2A11+E3*s3A11)+X*f1a*(E1*s4A11+E2*s5A11+E3*s6A11))/(h11))&
    	  &+((3*A3*(E1*t1A11+E2*t2A11+E3*t3A11)+X*f2a*(E1*t4A11+E2*t5A11+E3*t6A11))/(h22))
  
 K(9,9)= ((3*A3*(E1*s1A2+E2*s2A2+E3*s3A2)+X*f1a*(E1*s4A2+E2*s5A2+E3*s6A2))/(h11))&
         &+((3*A3*(E1*t1A2+E2*t2A2+E3*t3A2)+X*f2a*(E1*t4A2+E2*t5A2+E3*t6A2))/(h22))
				
 K(10,10)=((3*A3*(E1*s1A22+E2*s2A22+E3*s3A22)+X*f1a*(E1*s4A22+E2*s5A22+E3*s6A22))/(h11))&
    	    &+((3*A3*(E1*t1A22+E2*t2A22+E3*t3A22)+X*f2a*(E1*t4A22+E2*t5A22+E3*t6A22))/(h22))	  		
	   
 K(2,1)=((3*A3*(E2*s2B1+E3*s3B1)+X*f1a*(E2*s5B1+E3*s6B1))/(h11))&
         &+((3*A3*(E2*t2B1+E3*t3B1)+X*f2a*(E2*t5B1+E3*t6B1))/(h22))
  
 K(10,9)=((3*A3*(E2*s2B2+E3*s3B2)+X*f1a*(E2*s5B2+E3*s6B2))/(h11))&
         &+((3*A3*(E2*t2B2+E3*t3B2)+X*f2a*(E2*t5B2+E3*t6B2))/(h22)) 
	     	
 K(3,1) =((3*A3*(E1*s1C1+E2*s2C1+E3*s3C1)+X*f1a*(E1*s4C1+E2*s5C1+E3*s6C1))/(h11))&
        &+((3*A3*(E1*t1C1+E2*t2C1+E3*t3C1)+X*f2a*(E1*t4C1+E2*t5C1+E3*t6C1))/(h22))

 K(8,6) =((3*A3*(E1*s1C11+E2*s2C11+E3*s3C11)+X*f1a*(E1*s4C11+E2*s5C11+E3*s6C11))/(h11))&
        &+((3*A3*(E1*t1C11+E2*t2C11+E3*t3C11)+X*f2a*(E1*t4C11+E2*t5C11+E3*t6C11))/(h22))
  
 K(11,9) =((3*A3*(E1*s1C2+E2*s2C2+E3*s3C2)+X*f1a*(E1*s4C2+E2*s5C2+E3*s6C2))/(h11))&
         &+((3*A3*(E1*t1C2+E2*t2C2+E3*t3C2)+X*f2a*(E1*t4C2+E2*t5C2+E3*t6C2))/(h22))

 K(12,10) =((3*A3*(E1*s1C22+E2*s2C22+E3*s3C22)+X*f1a*(E1*s4C22+E2*s5C22+E3*s6C22))/(h11))&
          &+((3*A3*(E1*t1C22+E2*t2C22+E3*t3C22)+X*f2a*(E1*t4C22+E2*t5C22+E3*t6C22))/(h22))

 K(4,1) =((3*A3*(E1*s1D1+E2*s2D1+E3*s3D1)+X*f1a*(E1*s4D1+E2*s5D1+E3*s6D1))/(h11))&
        &+((3*A3*(E1*t1D1+E2*t2D1+E3*t3D1)+X*f2a*(E1*t4D1+E2*t5D1+E3*t6D1))/(h22))

 K(8,1) =((3*A3*(E1*s1D11+E2*s2D11+E3*s3D11)+X*f1a*(E1*s4D11+E2*s5D11+E3*s6D11))/(h11))&
        &+((3*A3*(E1*t1D11+E2*t2D11+E3*t3D11)+X*f2a*(E1*t4D11+E2*t5D11+E3*t6D11))/(h22))
  
 K(12,9) =((3*A3*(E1*s1D2+E2*s2D2+E3*s3D2)+X*f1a*(E1*s4D2+E2*s5D2+E3*s6D2))/(h11))&
         &+((3*A3*(E1*t1D2+E2*t2D2+E3*t3D2)+X*f2a*(E1*t4D2+E2*t5D2+E3*t6D2))/(h22))

 K(15,14) =((3*A3*(E1*s1D22+E2*s2D22+E3*s3D22)+X*f1a*(E1*s4D22+E2*s5D22+E3*s6D22))/(h11))&
          &+((3*A3*(E1*t1D22+E2*t2D22+E3*t3D22)+X*f2a*(E1*t4D22+E2*t5D22+E3*t6D22))/(h22))

 K(5,1) =((3*A3*(E1*s1E1+E2*s2E1+E3*s3E1)+X*f1a*(E1*s4E1+E2*s5E1+E3*s6E1))/(h11))&
        &+((3*A3*(E1*t1E1+E2*t2E1+E3*t3E1)+X*f2a*(E1*t4E1+E2*t5E1+E3*t6E1))/(h22))

 K(8,4) =((3*A3*(E1*s1E11+E2*s2E11+E3*s3E11)+X*f1a*(E1*s4E11+E2*s5E11+E3*s6E11))/(h11))&
        &+((3*A3*(E1*t1E11+E2*t2E11+E3*t3E11)+X*f2a*(E1*t4E11+E2*t5E11+E3*t6E11))/(h22))
  
 K(13,9) =((3*A3*(E1*s1E2+E2*s2E2+E3*s3E2)+X*f1a*(E1*s4E2+E2*s5E2+E3*s6E2))/(h11))&
         &+((3*A3*(E1*t1E2+E2*t2E2+E3*t3E2)+X*f2a*(E1*t4E2+E2*t5E2+E3*t6E2))/(h22))

 K(16,12) =((3*A3*(E1*s1E22+E2*s2E22+E3*s3E22)+X*f1a*(E1*s4E22+E2*s5E22+E3*s6E22))/(h11))&
          &+((3*A3*(E1*t1E22+E2*t2E22+E3*t3E22)+X*f2a*(E1*t4E22+E2*t5E22+E3*t6E22))/(h22))

 K(6,1) =((3*A3*(E1*s1F1+E2*s2F1+E3*s3F1)+X*f1a*(E1*s4F1+E2*s5F1+E3*s6F1))/(h11))&
        &+((3*A3*(E1*t1F1+E2*t2F1+E3*t3F1)+X*f2a*(E1*t4F1+E2*t5F1+E3*t6F1))/(h22))  
 
 K(14,9) =((3*A3*(E1*s1F2+E2*s2F2+E3*s3F2)+X*f1a*(E1*s4F2+E2*s5F2+E3*s6F2))/(h11))&
         &+((3*A3*(E1*t1F2+E2*t2F2+E3*t3F2)+X*f2a*(E1*t4F2+E2*t5F2+E3*t6F2))/(h22))

 K(15,7) =((3*A3*(E1*s1G1+E2*s2G1+E3*s3G1)+X*f1a*(E1*s4G1+E2*s5G1+E3*s6G1))/(h11))&
        &+((3*A3*(E1*t1G1+E2*t2G1+E3*t3G1)+X*f2a*(E1*t4G1+E2*t5G1+E3*t6G1))/(h22))

 K(13,7) =((3*A3*(E1*s1G11+E2*s2G11+E3*s3G11)+X*f1a*(E1*s4G11+E2*s5G11+E3*s6G11))/(h11))&
        &+((3*A3*(E1*t1G11+E2*t2G11+E3*t3G11)+X*f2a*(E1*t4G11+E2*t5G11+E3*t6G11))/(h22))
  
 K(10,2) =((3*A3*(E1*s1G2+E2*s2G2+E3*s3G2)+X*f1a*(E1*s4G2+E2*s5G2+E3*s6G2))/(h11))&
         &+((3*A3*(E1*t1G2+E2*t2G2+E3*t3G2)+X*f2a*(E1*t4G2+E2*t5G2+E3*t6G2))/(h22))

 K(12,6) =((3*A3*(E1*s1G22+E2*s2G22+E3*s3G22)+X*f1a*(E1*s4G22+E2*s5G22+E3*s6G22))/(h11))&
          &+((3*A3*(E1*t1G22+E2*t2G22+E3*t3G22)+X*f2a*(E1*t4G22+E2*t5G22+E3*t6G22))/(h22))
  
 K(9,2) =((3*A3*(E1*s1H1+E2*s2H1+E3*s3H1)+X*f1a*(E1*s4H1+E2*s5H1+E3*s6H1))/(h11))&
        &+((3*A3*(E1*t1H1+E2*t2H1+E3*t3H1)+X*f2a*(E1*t4H1+E2*t5H1+E3*t6H1))/(h22))

 K(11,6) =((3*A3*(E1*s1H11+E2*s2H11+E3*s3H11)+X*f1a*(E1*s4H11+E2*s5H11+E3*s6H11))/(h11))&
        &+((3*A3*(E1*t1H11+E2*t2H11+E3*t3H11)+X*f2a*(E1*t4H11+E2*t5H11+E3*t6H11))/(h22))
  
 K(10,1) =((3*A3*(E1*s1H2+E2*s2H2+E3*s3H2)+X*f1a*(E1*s4H2+E2*s5H2+E3*s6H2))/(h11))&
         &+((3*A3*(E1*t1H2+E2*t2H2+E3*t3H2)+X*f2a*(E1*t4H2+E2*t5H2+E3*t6H2))/(h22))

 K(12,5) =((3*A3*(E1*s1H22+E2*s2H22+E3*s3H22)+X*f1a*(E1*s4H22+E2*s5H22+E3*s6H22))/(h11))&
          &+((3*A3*(E1*t1H22+E2*t2H22+E3*t3H22)+X*f2a*(E1*t4H22+E2*t5H22+E3*t6H22))/(h22))
  
 K(11,1) =((3*A3*(E1*s1I1+E2*s2I1+E3*s3I1)+X*f1a*(E1*s4I1+E2*s5I1+E3*s6I1))/(h11))&
        &+((3*A3*(E1*t1I1+E2*t2I1+E3*t3I1)+X*f2a*(E1*t4I1+E2*t5I1+E3*t6I1))/(h22))

 K(9,5) =((3*A3*(E1*s1I11+E2*s2I11+E3*s3I11)+X*f1a*(E1*s4I11+E2*s5I11+E3*s6I11))/(h11))&
        &+((3*A3*(E1*t1I11+E2*t2I11+E3*t3I11)+X*f2a*(E1*t4I11+E2*t5I11+E3*t6I11))/(h22))
  
 K(12,2) =((3*A3*(E1*s1I2+E2*s2I2+E3*s3I2)+X*f1a*(E1*s4I2+E2*s5I2+E3*s6I2))/(h11))&
         &+((3*A3*(E1*t1I2+E2*t2I2+E3*t3I2)+X*f2a*(E1*t4I2+E2*t5I2+E3*t6I2))/(h22))

 K(10,6) =((3*A3*(E1*s1I22+E2*s2I22+E3*s3I22)+X*f1a*(E1*s4I22+E2*s5I22+E3*s6I22))/(h11))&
          &+((3*A3*(E1*t1I22+E2*t2I22+E3*t3I22)+X*f2a*(E1*t4I22+E2*t5I22+E3*t6I22))/(h22))
    
 K(11,2) =((3*A3*(E1*s1J1+E2*s2J1+E3*s3J1)+X*f1a*(E1*s4J1+E2*s5J1+E3*s6J1))/(h11))&
        &+((3*A3*(E1*t1J1+E2*t2J1+E3*t3J1)+X*f2a*(E1*t4J1+E2*t5J1+E3*t6J1))/(h22))

 K(13,2) =((3*A3*(E1*s1J11+E2*s2J11+E3*s3J11)+X*f1a*(E1*s4J11+E2*s5J11+E3*s6J11))/(h11))&
        &+((3*A3*(E1*t1J11+E2*t2J11+E3*t3J11)+X*f2a*(E1*t4J11+E2*t5J11+E3*t6J11))/(h22))
  
 K(12,1) =((3*A3*(E1*s1J2+E2*s2J2+E3*s3J2)+X*f1a*(E1*s4J2+E2*s5J2+E3*s6J2))/(h11))&
         &+((3*A3*(E1*t1J2+E2*t2J2+E3*t3J2)+X*f2a*(E1*t4J2+E2*t5J2+E3*t6J2))/(h22))

 K(14,1) =((3*A3*(E1*s1J22+E2*s2J22+E3*s3J22)+X*f1a*(E1*s4J22+E2*s5J22+E3*s6J22))/(h11))&
          &+((3*A3*(E1*t1J22+E2*t2J22+E3*t3J22)+X*f2a*(E1*t4J22+E2*t5J22+E3*t6J22))/(h22))

 

		case(2)		


  	K(7,7) = ((3*A3*(E1*s1A1+E2*s2A1+E3*s3A1)+X*f1*(E1*s4A1+E2*s5A1+E3*s6A1))/(h22))&
       &+((3*A3*(E1*t1A1+E2*t2A1+E3*t3A1)+X*f2*(E1*t4A1+E2*t5A1+E3*t6A1))/(h11))
  			
	K(8,8)=((3*A3*(E1*s1A11+E2*s2A11+E3*s3A11)+X*f1*(E1*s4A11+E2*s5A11+E3*s6A11))/(h22))&
	  &+((3*A3*(E1*t1A11+E2*t2A11+E3*t3A11)+X*f2*(E1*t4A11+E2*t5A11+E3*t6A11))/(h11))
	   
   	K(15,15)=((3*A3*(E1*s1A2+E2*s2A2+E3*s3A2)+X*f1*(E1*s4A2+E2*s5A2+E3*s6A2))/(h22))&
       &+((3*A3*(E1*t1A2+E2*t2A2+E3*t3A2)+X*f2*(E1*t4A2+E2*t5A2+E3*t6A2))/(h11))
  			
	K(16,16)=((3*A3*(E1*s1A22+E2*s2A22+E3*s3A22)+X*f1*(E1*s4A22+E2*s5A22+E3*s6A22))/(h22))&
	  &+((3*A3*(E1*t1A22+E2*t2A22+E3*t3A22)+X*f2*(E1*t4A22+E2*t5A22+E3*t6A22))/(h11))
	   
    K(8,7) = ((3*A3*(E1*s1B1+E2*s2B1+E3*s3B1)+X*f1*(E1*s4B1+E2*s5B1+E3*s6B1))/(h22))&
       &+((3*A3*(E1*t1B1+E2*t2B1+E3*t3B1)+X*f2*(E1*t4B1+E2*t5B1+E3*t6B1))/(h11))
	
	K(16,15)=((3*A3*(E1*s1B2+E2*s2B2+E3*s3B2)+X*f1*(E1*s4B2+E2*s5B2+E3*s6B2))/(h22))&
       &+((3*A3*(E1*t1B2+E2*t2B2+E3*t3B2)+X*f2*(E1*t4B2+E2*t5B2+E3*t6B2))/(h11))

	K(7,1) = ((3*A3*(E1*s1C1+E2*s2C1+E3*s3C1)+X*f1*(E1*s4C1+E2*s5C1+E3*s6C1))/(h22))&
       &+((3*A3*(E1*t1C1+E2*t2C1+E3*t3C1)+X*f2*(E1*t4C1+E2*t5C1+E3*t6C1))/(h11))
  			
	K(6,4)=((3*A3*(E1*s1C11+E2*s2C11+E3*s3C11)+X*f1*(E1*s4C11+E2*s5C11+E3*s6C11))/(h22))&
	  &+((3*A3*(E1*t1C11+E2*t2C11+E3*t3C11)+X*f2*(E1*t4C11+E2*t5C11+E3*t6C11))/(h11))
	   
   	K(15,9)=((3*A3*(E1*s1C2+E2*s2C2+E3*s3C2)+X*f1*(E1*s4C2+E2*s5C2+E3*s6C2))/(h22))&
       &+((3*A3*(E1*t1C2+E2*t2C2+E3*t3C2)+X*f2*(E1*t4C2+E2*t5C2+E3*t6C2))/(h11))
  			
	K(16,10)=((3*A3*(E1*s1C22+E2*s2C22+E3*s3C22)+X*f1*(E1*s4C22+E2*s5C22+E3*s6C22))/(h22))&
	  &+((3*A3*(E1*t1C22+E2*t2C22+E3*t3C22)+X*f2*(E1*t4C22+E2*t5C22+E3*t6C22))/(h11))

   K(7,2) = ((3*A3*(E1*s1D1+E2*s2D1+E3*s3D1)+X*f1*(E1*s4D1+E2*s5D1+E3*s6D1))/(h22))&
       &+((3*A3*(E1*t1D1+E2*t2D1+E3*t3D1)+X*f2*(E1*t4D1+E2*t5D1+E3*t6D1))/(h11))

  K(7,3) =((3*A3*(E1*s1E1+E2*s2E1+E3*s3E1)+X*f1*(E1*s4E1+E2*s5E1+E3*s6E1))/(h22))&
        &+((3*A3*(E1*t1E1+E2*t2E1+E3*t3E1)+X*f2*(E1*t4E1+E2*t5E1+E3*t6E1))/(h11))


  K(6,2) =((3*A3*(E1*s1E11+E2*s2E11+E3*s3E11)+X*f1*(E1*s4E11+E2*s5E11+E3*s6E11))/(h22))&
        &+((3*A3*(E1*t1E11+E2*t2E11+E3*t3E11)+X*f2*(E1*t4E11+E2*t5E11+E3*t6E11))/(h11))

  K(15,11) =((3*A3*(E1*s1E2+E2*s2E2+E3*s3E2)+X*f1*(E1*s4E2+E2*s5E2+E3*s6E2))/(h22))&
         &+((3*A3*(E1*t1E2+E2*t2E2+E3*t3E2)+X*f2*(E1*t4E2+E2*t5E2+E3*t6E2))/(h11))

  K(14,10) =((3*A3*(E1*s1E22+E2*s2E22+E3*s3E22)+X*f1*(E1*s4E22+E2*s5E22+E3*s6E22))/(h22))&
          &+((3*A3*(E1*t1E22+E2*t2E22+E3*t3E22)+X*f2*(E1*t4E22+E2*t5E22+E3*t6E22))/(h11))


  			
   K(7,6)=((3*A3*(E1*s1D11+E2*s2D11+E3*s3D11)+X*f1*(E1*s4D11+E2*s5D11+E3*s6D11))/(h22))&
	  &+((3*A3*(E1*t1D11+E2*t2D11+E3*t3D11)+X*f2*(E1*t4D11+E2*t5D11+E3*t6D11))/(h11))
	   
   K(15,10)=((3*A3*(E1*s1D2+E2*s2D2+E3*s3D2)+X*f1*(E1*s4D2+E2*s5D2+E3*s6D2))/(h22))&
       &+((3*A3*(E1*t1D2+E2*t2D2+E3*t3D2)+X*f2*(E1*t4D2+E2*t5D2+E3*t6D2))/(h11))
  			
	K(13,12)=((3*A3*(E1*s1D22+E2*s2D22+E3*s3D22)+X*f1*(E1*s4D22+E2*s5D22+E3*s6D22))/(h22))&
	  &+((3*A3*(E1*t1D22+E2*t2D22+E3*t3D22)+X*f2*(E1*t4D22+E2*t5D22+E3*t6D22))/(h11))
	    
    K(7,4) = ((3*A3*(E1*s1F1+E2*s2F1+E3*s3F1)+X*f1*(E1*s4F1+E2*s5F1+E3*s6F1))/(h22))&
       &+((3*A3*(E1*t1F1+E2*t2F1+E3*t3F1)+X*f2*(E1*t4F1+E2*t5F1+E3*t6F1))/(h11))
	
	K(15,12)=((3*A3*(E1*s1F2+E2*s2F2+E3*s3F2)+X*f1*(E1*s4F2+E2*s5F2+E3*s6F2))/(h22))&
       &+((3*A3*(E1*t1F2+E2*t2F2+E3*t3F2)+X*f2*(E1*t4F2+E2*t5F2+E3*t6F2))/(h11))
   
  K(13,5) = ((3*A3*(E1*s1G1+E2*s2G1+E3*s3G1)+X*f1*(E1*s4G1+E2*s5G1+E3*s6G1))/(h22))&
       &+((3*A3*(E1*t1G1+E2*t2G1+E3*t3G1)+X*f2*(E1*t4G1+E2*t5G1+E3*t6G1))/(h11))
  			
  K(11,5)=((3*A3*(E1*s1G11+E2*s2G11+E3*s3G11)+X*f1*(E1*s4G11+E2*s5G11+E3*s6G11))/(h22))&
	  &+((3*A3*(E1*t1G11+E2*t2G11+E3*t3G11)+X*f2*(E1*t4G11+E2*t5G11+E3*t6G11))/(h11))
	   
  K(16,8)=((3*A3*(E1*s1G2+E2*s2G2+E3*s3G2)+X*f1*(E1*s4G2+E2*s5G2+E3*s6G2))/(h22))&
       &+((3*A3*(E1*t1G2+E2*t2G2+E3*t3G2)+X*f2*(E1*t4G2+E2*t5G2+E3*t6G2))/(h11))
  			
  K(10,4)=((3*A3*(E1*s1G22+E2*s2G22+E3*s3G22)+X*f1*(E1*s4G22+E2*s5G22+E3*s6G22))/(h22))&
	  &+((3*A3*(E1*t1G22+E2*t2G22+E3*t3G22)+X*f2*(E1*t4G22+E2*t5G22+E3*t6G22))/(h11))
		
   K(15,8) = ((3*A3*(E1*s1H1+E2*s2H1+E3*s3H1)+X*f1*(E1*s4H1+E2*s5H1+E3*s6H1))/(h22))&
          &+((3*A3*(E1*t1H1+E2*t2H1+E3*t3H1)+X*f2*(E1*t4H1+E2*t5H1+E3*t6H1))/(h11))
  			
   K(9,4)=((3*A3*(E1*s1H11+E2*s2H11+E3*s3H11)+X*f1*(E1*s4H11+E2*s5H11+E3*s6H11))/(h22))&
	     &+((3*A3*(E1*t1H11+E2*t2H11+E3*t3H11)+X*f2*(E1*t4H11+E2*t5H11+E3*t6H11))/(h11))
	   
   K(16,7)=((3*A3*(E1*s1H2+E2*s2H2+E3*s3H2)+X*f1*(E1*s4H2+E2*s5H2+E3*s6H2))/(h22))&
          &+((3*A3*(E1*t1H2+E2*t2H2+E3*t3H2)+X*f2*(E1*t4H2+E2*t5H2+E3*t6H2))/(h11))
  			
   K(10,3)=((3*A3*(E1*s1H22+E2*s2H22+E3*s3H22)+X*f1*(E1*s4H22+E2*s5H22+E3*s6H22))/(h22))&
	      &+((3*A3*(E1*t1H22+E2*t2H22+E3*t3H22)+X*f2*(E1*t4H22+E2*t5H22+E3*t6H22))/(h11))
   
  K(9,7) = ((3*A3*(E1*s1I1+E2*s2I1+E3*s3I1)+X*f1*(E1*s4I1+E2*s5I1+E3*s6I1))/(h22))&
       &+((3*A3*(E1*t1I1+E2*t2I1+E3*t3I1)+X*f2*(E1*t4I1+E2*t5I1+E3*t6I1))/(h11))
  			
  K(15,3)=((3*A3*(E1*s1I11+E2*s2I11+E3*s3I11)+X*f1*(E1*s4I11+E2*s5I11+E3*s6I11))/(h22))&
	  &+((3*A3*(E1*t1I11+E2*t2I11+E3*t3I11)+X*f2*(E1*t4I11+E2*t5I11+E3*t6I11))/(h11))
	   
  K(10,8)=((3*A3*(E1*s1I2+E2*s2I2+E3*s3I2)+X*f1*(E1*s4I2+E2*s5I2+E3*s6I2))/(h22))&
       &+((3*A3*(E1*t1I2+E2*t2I2+E3*t3I2)+X*f2*(E1*t4I2+E2*t5I2+E3*t6I2))/(h11))
  			
  K(16,4)=((3*A3*(E1*s1I22+E2*s2I22+E3*s3I22)+X*f1*(E1*s4I22+E2*s5I22+E3*s6I22))/(h22))&
	  &+((3*A3*(E1*t1I22+E2*t2I22+E3*t3I22)+X*f2*(E1*t4I22+E2*t5I22+E3*t6I22))/(h11))


  K(9,8) = ((3*A3*(E1*s1J1+E2*s2J1+E3*s3J1)+X*f1*(E1*s4J1+E2*s5J1+E3*s6J1))/(h22))&
        &+((3*A3*(E1*t1J1+E2*t2J1+E3*t3J1)+X*f2*(E1*t4J1+E2*t5J1+E3*t6J1))/(h11))
  	 		
  K(11,8)=((3*A3*(E1*s1J11+E2*s2J11+E3*s3J11)+X*f1*(E1*s4J11+E2*s5J11+E3*s6J11))/(h22))&
	   &+((3*A3*(E1*t1J11+E2*t2J11+E3*t3J11)+X*f2*(E1*t4J11+E2*t5J11+E3*t6J11))/(h11))
	   
  K(10,7)=((3*A3*(E1*s1J2+E2*s2J2+E3*s3J2)+X*f1*(E1*s4J2+E2*s5J2+E3*s6J2))/(h22))&
       &+((3*A3*(E1*t1J2+E2*t2J2+E3*t3J2)+X*f2*(E1*t4J2+E2*t5J2+E3*t6J2))/(h11))
  			
  K(12,7)=((3*A3*(E1*s1J22+E2*s2J22+E3*s3J22)+X*f1*(E1*s4J22+E2*s5J22+E3*s6J22))/(h22))&
	    &+((3*A3*(E1*t1J22+E2*t2J22+E3*t3J22)+X*f2*(E1*t4J22+E2*t5J22+E3*t6J22))/(h11))

    
		case(3)
		
  
 K(5,5)=((3*A3*(E1*s1A1+E2*s2A1+E3*s3A1)+X*f1*(E1*s4A1+E2*s5A1+E3*s6A1))/(h11))&
        &+((3*A3*(E1*t1A1+E2*t2A1+E3*t3A1)+X*f2*(E1*t4A1+E2*t5A1+E3*t6A1))/(h22))
				
 K(6,6)=((3*A3*(E1*s1A11+E2*s2A11+E3*s3A11)+X*f1*(E1*s4A11+E2*s5A11+E3*s6A11))/(h11))&
    	  &+((3*A3*(E1*t1A11+E2*t2A11+E3*t3A11)+X*f2*(E1*t4A11+E2*t5A11+E3*t6A11))/(h22))
  
 K(13,13)= ((3*A3*(E1*s1A2+E2*s2A2+E3*s3A2)+X*f1*(E1*s4A2+E2*s5A2+E3*s6A2))/(h11))&
         &+((3*A3*(E1*t1A2+E2*t2A2+E3*t3A2)+X*f2*(E1*t4A2+E2*t5A2+E3*t6A2))/(h22))
				
 K(14,14)=((3*A3*(E1*s1A22+E2*s2A22+E3*s3A22)+X*f1*(E1*s4A22+E2*s5A22+E3*s6A22))/(h11))&
    	    &+((3*A3*(E1*t1A22+E2*t2A22+E3*t3A22)+X*f2*(E1*t4A22+E2*t5A22+E3*t6A22))/(h22))	  		
	   
 K(6,5)=((3*A3*(E2*s2B1+E3*s3B1)+X*f1*(E2*s5B1+E3*s6B1))/(h11))&
         &+((3*A3*(E2*t2B1+E3*t3B1)+X*f2*(E2*t5B1+E3*t6B1))/(h22))
  
! K(14,13)=((3*A3*(E2*s2B1+E3*s3B1)+X*f1*(E2*s5B1+E3*s6B1))/(h11))&
!         &+((3*A3*(E2*t2B1+E3*t3B1)+X*f2*(E2*t5B1+E3*t6B1))/(h22)) 

 K(14,13)=((3*A3*(E2*s2B2+E3*s3B2)+X*f1*(E2*s5B2+E3*s6B2))/(h11))&
         &+((3*A3*(E2*t2B2+E3*t3B2)+X*f2*(E2*t5B2+E3*t6B2))/(h22))
	     	
 K(7,5) =((3*A3*(E1*s1C1+E2*s2C1+E3*s3C1)+X*f1*(E1*s4C1+E2*s5C1+E3*s6C1))/(h11))&
        &+((3*A3*(E1*t1C1+E2*t2C1+E3*t3C1)+X*f2*(E1*t4C1+E2*t5C1+E3*t6C1))/(h22))

 K(4,2) =((3*A3*(E1*s1C11+E2*s2C11+E3*s3C11)+X*f1*(E1*s4C11+E2*s5C11+E3*s6C11))/(h11))&
        &+((3*A3*(E1*t1C11+E2*t2C11+E3*t3C11)+X*f2*(E1*t4C11+E2*t5C11+E3*t6C11))/(h22))
  
 K(15,13) =((3*A3*(E1*s1C2+E2*s2C2+E3*s3C2)+X*f1*(E1*s4C2+E2*s5C2+E3*s6C2))/(h11))&
         &+((3*A3*(E1*t1C2+E2*t2C2+E3*t3C2)+X*f2*(E1*t4C2+E2*t5C2+E3*t6C2))/(h22))

 K(16,14) =((3*A3*(E1*s1C22+E2*s2C22+E3*s3C22)+X*f1*(E1*s4C22+E2*s5C22+E3*s6C22))/(h11))&
          &+((3*A3*(E1*t1C22+E2*t2C22+E3*t3C22)+X*f2*(E1*t4C22+E2*t5C22+E3*t6C22))/(h22))

 K(8,5) =((3*A3*(E1*s1D1+E2*s2D1+E3*s3D1)+X*f1*(E1*s4D1+E2*s5D1+E3*s6D1))/(h11))&
        &+((3*A3*(E1*t1D1+E2*t2D1+E3*t3D1)+X*f2*(E1*t4D1+E2*t5D1+E3*t6D1))/(h22))

 K(5,4) =((3*A3*(E1*s1D11+E2*s2D11+E3*s3D11)+X*f1*(E1*s4D11+E2*s5D11+E3*s6D11))/(h11))&
        &+((3*A3*(E1*t1D11+E2*t2D11+E3*t3D11)+X*f2*(E1*t4D11+E2*t5D11+E3*t6D11))/(h22))
  
 K(16,13) =((3*A3*(E1*s1D2+E2*s2D2+E3*s3D2)+X*f1*(E1*s4D2+E2*s5D2+E3*s6D2))/(h11))&
         &+((3*A3*(E1*t1D2+E2*t2D2+E3*t3D2)+X*f2*(E1*t4D2+E2*t5D2+E3*t6D2))/(h22))

 K(11,10) =((3*A3*(E1*s1D22+E2*s2D22+E3*s3D22)+X*f1*(E1*s4D22+E2*s5D22+E3*s6D22))/(h11))&
          &+((3*A3*(E1*t1D22+E2*t2D22+E3*t3D22)+X*f2*(E1*t4D22+E2*t5D22+E3*t6D22))/(h22))
  

 K(5,2) =((3*A3*(E1*s1F1+E2*s2F1+E3*s3F1)+X*f1*(E1*s4F1+E2*s5F1+E3*s6F1))/(h11))&
        &+((3*A3*(E1*t1F1+E2*t2F1+E3*t3F1)+X*f2*(E1*t4F1+E2*t5F1+E3*t6F1))/(h22))  
 
 K(13,10) =((3*A3*(E1*s1F2+E2*s2F2+E3*s3F2)+X*f1*(E1*s4F2+E2*s5F2+E3*s6F2))/(h11))&
         &+((3*A3*(E1*t1F2+E2*t2F2+E3*t3F2)+X*f2*(E1*t4F2+E2*t5F2+E3*t6F2))/(h22))

 K(11,3) =((3*A3*(E1*s1G1+E2*s2G1+E3*s3G1)+X*f1*(E1*s4G1+E2*s5G1+E3*s6G1))/(h11))&
        &+((3*A3*(E1*t1G1+E2*t2G1+E3*t3G1)+X*f2*(E1*t4G1+E2*t5G1+E3*t6G1))/(h22))

 K(9,3) =((3*A3*(E1*s1G11+E2*s2G11+E3*s3G11)+X*f1*(E1*s4G11+E2*s5G11+E3*s6G11))/(h11))&
        &+((3*A3*(E1*t1G11+E2*t2G11+E3*t3G11)+X*f2*(E1*t4G11+E2*t5G11+E3*t6G11))/(h22))
  
 K(14,6) =((3*A3*(E1*s1G2+E2*s2G2+E3*s3G2)+X*f1*(E1*s4G2+E2*s5G2+E3*s6G2))/(h11))&
         &+((3*A3*(E1*t1G2+E2*t2G2+E3*t3G2)+X*f2*(E1*t4G2+E2*t5G2+E3*t6G2))/(h22))

 K(16,2) =((3*A3*(E1*s1G22+E2*s2G22+E3*s3G22)+X*f1*(E1*s4G22+E2*s5G22+E3*s6G22))/(h11))&
          &+((3*A3*(E1*t1G22+E2*t2G22+E3*t3G22)+X*f2*(E1*t4G22+E2*t5G22+E3*t6G22))/(h22))
  
 K(13,6) =((3*A3*(E1*s1H1+E2*s2H1+E3*s3H1)+X*f1*(E1*s4H1+E2*s5H1+E3*s6H1))/(h11))&
        &+((3*A3*(E1*t1H1+E2*t2H1+E3*t3H1)+X*f2*(E1*t4H1+E2*t5H1+E3*t6H1))/(h22))

 K(15,2) =((3*A3*(E1*s1H11+E2*s2H11+E3*s3H11)+X*f1*(E1*s4H11+E2*s5H11+E3*s6H11))/(h11))&
        &+((3*A3*(E1*t1H11+E2*t2H11+E3*t3H11)+X*f2*(E1*t4H11+E2*t5H11+E3*t6H11))/(h22))
  
 K(14,5) =((3*A3*(E1*s1H2+E2*s2H2+E3*s3H2)+X*f1*(E1*s4H2+E2*s5H2+E3*s6H2))/(h11))&
         &+((3*A3*(E1*t1H2+E2*t2H2+E3*t3H2)+X*f2*(E1*t4H2+E2*t5H2+E3*t6H2))/(h22))

 K(16,1) =((3*A3*(E1*s1H22+E2*s2H22+E3*s3H22)+X*f1*(E1*s4H22+E2*s5H22+E3*s6H22))/(h11))&
          &+((3*A3*(E1*t1H22+E2*t2H22+E3*t3H22)+X*f2*(E1*t4H22+E2*t5H22+E3*t6H22))/(h22))
  
 K(15,5) =((3*A3*(E1*s1I1+E2*s2I1+E3*s3I1)+X*f1*(E1*s4I1+E2*s5I1+E3*s6I1))/(h11))&
        &+((3*A3*(E1*t1I1+E2*t2I1+E3*t3I1)+X*f2*(E1*t4I1+E2*t5I1+E3*t6I1))/(h22))

 K(13,1) =((3*A3*(E1*s1I11+E2*s2I11+E3*s3I11)+X*f1*(E1*s4I11+E2*s5I11+E3*s6I11))/(h11))&
        &+((3*A3*(E1*t1I11+E2*t2I11+E3*t3I11)+X*f2*(E1*t4I11+E2*t5I11+E3*t6I11))/(h22))
  
 K(16,6) =((3*A3*(E1*s1I2+E2*s2I2+E3*s3I2)+X*f1*(E1*s4I2+E2*s5I2+E3*s6I2))/(h11))&
         &+((3*A3*(E1*t1I2+E2*t2I2+E3*t3I2)+X*f2*(E1*t4I2+E2*t5I2+E3*t6I2))/(h22))

 K(14,2) =((3*A3*(E1*s1I22+E2*s2I22+E3*s3I22)+X*f1*(E1*s4I22+E2*s5I22+E3*s6I22))/(h11))&
          &+((3*A3*(E1*t1I22+E2*t2I22+E3*t3I22)+X*f2*(E1*t4I22+E2*t5I22+E3*t6I22))/(h22))
    
 K(15,6) =((3*A3*(E1*s1J1+E2*s2J1+E3*s3J1)+X*f1*(E1*s4J1+E2*s5J1+E3*s6J1))/(h11))&
        &+((3*A3*(E1*t1J1+E2*t2J1+E3*t3J1)+X*f2*(E1*t4J1+E2*t5J1+E3*t6J1))/(h22))

 K(9,6) =((3*A3*(E1*s1J11+E2*s2J11+E3*s3J11)+X*f1*(E1*s4J11+E2*s5J11+E3*s6J11))/(h11))&
        &+((3*A3*(E1*t1J11+E2*t2J11+E3*t3J11)+X*f2*(E1*t4J11+E2*t5J11+E3*t6J11))/(h22))
  
 K(16,5) =((3*A3*(E1*s1J2+E2*s2J2+E3*s3J2)+X*f1*(E1*s4J2+E2*s5J2+E3*s6J2))/(h11))&
         &+((3*A3*(E1*t1J2+E2*t2J2+E3*t3J2)+X*f2*(E1*t4J2+E2*t5J2+E3*t6J2))/(h22))

 K(10,5) =((3*A3*(E1*s1J22+E2*s2J22+E3*s3J22)+X*f1*(E1*s4J22+E2*s5J22+E3*s6J22))/(h11))&
          &+((3*A3*(E1*t1J22+E2*t2J22+E3*t3J22)+X*f2*(E1*t4J22+E2*t5J22+E3*t6J22))/(h22))

 

	
		case(4)
		
	   	
   K(3,3) = ((3*A3*(E1*s1A1+E2*s2A1+E3*s3A1)+X*f1*(E1*s4A1+E2*s5A1+E3*s6A1))/(h22))&
         &+((3*A3*(E1*t1A1+E2*t2A1+E3*t3A1)+X*f2*(E1*t4A1+E2*t5A1+E3*t6A1))/(h11))
  	 		
	K(4,4)=((3*A3*(E1*s1A11+E2*s2A11+E3*s3A11)+X*f1*(E1*s4A11+E2*s5A11+E3*s6A11))/(h22))&
	  &+((3*A3*(E1*t1A11+E2*t2A11+E3*t3A11)+X*f2*(E1*t4A11+E2*t5A11+E3*t6A11))/(h11))
	   
   	K(11,11)=((3*A3*(E1*s1A2+E2*s2A2+E3*s3A2)+X*f1*(E1*s4A2+E2*s5A2+E3*s6A2))/(h22))&
       &+((3*A3*(E1*t1A2+E2*t2A2+E3*t3A2)+X*f2*(E1*t4A2+E2*t5A2+E3*t6A2))/(h11))
  			
	K(12,12)=((3*A3*(E1*s1A22+E2*s2A22+E3*s3A22)+X*f1*(E1*s4A22+E2*s5A22+E3*s6A22))/(h22))&
	  &+((3*A3*(E1*t1A22+E2*t2A22+E3*t3A22)+X*f2*(E1*t4A22+E2*t5A22+E3*t6A22))/(h11))
	   
    K(4,3) = ((3*A3*(E1*s1B1+E2*s2B1+E3*s3B1)+X*f1*(E1*s4B1+E2*s5B1+E3*s6B1))/(h22))&
       &+((3*A3*(E1*t1B1+E2*t2B1+E3*t3B1)+X*f2*(E1*t4B1+E2*t5B1+E3*t6B1))/(h11))
	
	K(12,11)=((3*A3*(E1*s1B2+E2*s2B2+E3*s3B2)+X*f1*(E1*s4B2+E2*s5B2+E3*s6B2))/(h22))&
       &+((3*A3*(E1*t1B2+E2*t2B2+E3*t3B2)+X*f2*(E1*t4B2+E2*t5B2+E3*t6B2))/(h11))

	K(5,3) = ((3*A3*(E1*s1C1+E2*s2C1+E3*s3C1)+X*f1*(E1*s4C1+E2*s5C1+E3*s6C1))/(h22))&
       &+((3*A3*(E1*t1C1+E2*t2C1+E3*t3C1)+X*f2*(E1*t4C1+E2*t5C1+E3*t6C1))/(h11))
  			
	K(8,2)=((3*A3*(E1*s1C11+E2*s2C11+E3*s3C11)+X*f1*(E1*s4C11+E2*s5C11+E3*s6C11))/(h22))&
	  &+((3*A3*(E1*t1C11+E2*t2C11+E3*t3C11)+X*f2*(E1*t4C11+E2*t5C11+E3*t6C11))/(h11))
	   
   	K(13,11)=((3*A3*(E1*s1C2+E2*s2C2+E3*s3C2)+X*f1*(E1*s4C2+E2*s5C2+E3*s6C2))/(h22))&
       &+((3*A3*(E1*t1C2+E2*t2C2+E3*t3C2)+X*f2*(E1*t4C2+E2*t5C2+E3*t6C2))/(h11))
  			
	K(14,12)=((3*A3*(E1*s1C22+E2*s2C22+E3*s3C22)+X*f1*(E1*s4C22+E2*s5C22+E3*s6C22))/(h22))&
	  &+((3*A3*(E1*t1C22+E2*t2C22+E3*t3C22)+X*f2*(E1*t4C22+E2*t5C22+E3*t6C22))/(h11))

    K(6,3) = ((3*A3*(E1*s1D1+E2*s2D1+E3*s3D1)+X*f1*(E1*s4D1+E2*s5D1+E3*s6D1))/(h22))&
       &+((3*A3*(E1*t1D1+E2*t2D1+E3*t3D1)+X*f2*(E1*t4D1+E2*t5D1+E3*t6D1))/(h11))
  			
   K(3,2)=((3*A3*(E1*s1D11+E2*s2D11+E3*s3D11)+X*f1*(E1*s4D11+E2*s5D11+E3*s6D11))/(h22))&
	  &+((3*A3*(E1*t1D11+E2*t2D11+E3*t3D11)+X*f2*(E1*t4D11+E2*t5D11+E3*t6D11))/(h11))
	   
   K(14,11)=((3*A3*(E1*s1D2+E2*s2D2+E3*s3D2)+X*f1*(E1*s4D2+E2*s5D2+E3*s6D2))/(h22))&
       &+((3*A3*(E1*t1D2+E2*t2D2+E3*t3D2)+X*f2*(E1*t4D2+E2*t5D2+E3*t6D2))/(h11))
  			
	K(16,9)=((3*A3*(E1*s1D22+E2*s2D22+E3*s3D22)+X*f1*(E1*s4D22+E2*s5D22+E3*s6D22))/(h22))&
	  &+((3*A3*(E1*t1D22+E2*t2D22+E3*t3D22)+X*f2*(E1*t4D22+E2*t5D22+E3*t6D22))/(h11))

    
    K(8,3) = ((3*A3*(E1*s1F1+E2*s2F1+E3*s3F1)+X*f1*(E1*s4F1+E2*s5F1+E3*s6F1))/(h22))&
       &+((3*A3*(E1*t1F1+E2*t2F1+E3*t3F1)+X*f2*(E1*t4F1+E2*t5F1+E3*t6F1))/(h11))
	
	K(16,11)=((3*A3*(E1*s1F2+E2*s2F2+E3*s3F2)+X*f1*(E1*s4F2+E2*s5F2+E3*s6F2))/(h22))&
       &+((3*A3*(E1*t1F2+E2*t2F2+E3*t3F2)+X*f2*(E1*t4F2+E2*t5F2+E3*t6F2))/(h11))
   
  K(9,1) = ((3*A3*(E1*s1G1+E2*s2G1+E3*s3G1)+X*f1*(E1*s4G1+E2*s5G1+E3*s6G1))/(h22))&
       &+((3*A3*(E1*t1G1+E2*t2G1+E3*t3G1)+X*f2*(E1*t4G1+E2*t5G1+E3*t6G1))/(h11))
  			
  K(15,1)=((3*A3*(E1*s1G11+E2*s2G11+E3*s3G11)+X*f1*(E1*s4G11+E2*s5G11+E3*s6G11))/(h22))&
	  &+((3*A3*(E1*t1G11+E2*t2G11+E3*t3G11)+X*f2*(E1*t4G11+E2*t5G11+E3*t6G11))/(h11))
	   
  K(12,4)=((3*A3*(E1*s1G2+E2*s2G2+E3*s3G2)+X*f1*(E1*s4G2+E2*s5G2+E3*s6G2))/(h22))&
       &+((3*A3*(E1*t1G2+E2*t2G2+E3*t3G2)+X*f2*(E1*t4G2+E2*t5G2+E3*t6G2))/(h11))
  			
  K(14,8)=((3*A3*(E1*s1G22+E2*s2G22+E3*s3G22)+X*f1*(E1*s4G22+E2*s5G22+E3*s6G22))/(h22))&
	  &+((3*A3*(E1*t1G22+E2*t2G22+E3*t3G22)+X*f2*(E1*t4G22+E2*t5G22+E3*t6G22))/(h11))
		
   K(11,4) = ((3*A3*(E1*s1H1+E2*s2H1+E3*s3H1)+X*f1*(E1*s4H1+E2*s5H1+E3*s6H1))/(h22))&
          &+((3*A3*(E1*t1H1+E2*t2H1+E3*t3H1)+X*f2*(E1*t4H1+E2*t5H1+E3*t6H1))/(h11))
  			
   K(13,8)=((3*A3*(E1*s1H11+E2*s2H11+E3*s3H11)+X*f1*(E1*s4H11+E2*s5H11+E3*s6H11))/(h22))&
	     &+((3*A3*(E1*t1H11+E2*t2H11+E3*t3H11)+X*f2*(E1*t4H11+E2*t5H11+E3*t6H11))/(h11))
	   
   K(12,3)=((3*A3*(E1*s1H2+E2*s2H2+E3*s3H2)+X*f1*(E1*s4H2+E2*s5H2+E3*s6H2))/(h22))&
          &+((3*A3*(E1*t1H2+E2*t2H2+E3*t3H2)+X*f2*(E1*t4H2+E2*t5H2+E3*t6H2))/(h11))
  			
   K(14,7)=((3*A3*(E1*s1H22+E2*s2H22+E3*s3H22)+X*f1*(E1*s4H22+E2*s5H22+E3*s6H22))/(h22))&
	      &+((3*A3*(E1*t1H22+E2*t2H22+E3*t3H22)+X*f2*(E1*t4H22+E2*t5H22+E3*t6H22))/(h11))
   
  K(13,3) = ((3*A3*(E1*s1I1+E2*s2I1+E3*s3I1)+X*f1*(E1*s4I1+E2*s5I1+E3*s6I1))/(h22))&
       &+((3*A3*(E1*t1I1+E2*t2I1+E3*t3I1)+X*f2*(E1*t4I1+E2*t5I1+E3*t6I1))/(h11))
  			
  K(11,7)=((3*A3*(E1*s1I11+E2*s2I11+E3*s3I11)+X*f1*(E1*s4I11+E2*s5I11+E3*s6I11))/(h22))&
	  &+((3*A3*(E1*t1I11+E2*t2I11+E3*t3I11)+X*f2*(E1*t4I11+E2*t5I11+E3*t6I11))/(h11))
	   
  K(14,4)=((3*A3*(E1*s1I2+E2*s2I2+E3*s3I2)+X*f1*(E1*s4I2+E2*s5I2+E3*s6I2))/(h22))&
       &+((3*A3*(E1*t1I2+E2*t2I2+E3*t3I2)+X*f2*(E1*t4I2+E2*t5I2+E3*t6I2))/(h11))
  			
  K(12,8)=((3*A3*(E1*s1I22+E2*s2I22+E3*s3I22)+X*f1*(E1*s4I22+E2*s5I22+E3*s6I22))/(h22))&
	  &+((3*A3*(E1*t1I22+E2*t2I22+E3*t3I22)+X*f2*(E1*t4I22+E2*t5I22+E3*t6I22))/(h11))


  K(13,4) = ((3*A3*(E1*s1J1+E2*s2J1+E3*s3J1)+X*f1*(E1*s4J1+E2*s5J1+E3*s6J1))/(h22))&
        &+((3*A3*(E1*t1J1+E2*t2J1+E3*t3J1)+X*f2*(E1*t4J1+E2*t5J1+E3*t6J1))/(h11))
  	 		
  K(15,4)=((3*A3*(E1*s1J11+E2*s2J11+E3*s3J11)+X*f1*(E1*s4J11+E2*s5J11+E3*s6J11))/(h22))&
	   &+((3*A3*(E1*t1J11+E2*t2J11+E3*t3J11)+X*f2*(E1*t4J11+E2*t5J11+E3*t6J11))/(h11))
	   
  K(14,3)=((3*A3*(E1*s1J2+E2*s2J2+E3*s3J2)+X*f1*(E1*s4J2+E2*s5J2+E3*s6J2))/(h22))&
       &+((3*A3*(E1*t1J2+E2*t2J2+E3*t3J2)+X*f2*(E1*t4J2+E2*t5J2+E3*t6J2))/(h11))
  			
  K(16,3)=((3*A3*(E1*s1J22+E2*s2J22+E3*s3J22)+X*f1*(E1*s4J22+E2*s5J22+E3*s6J22))/(h22))&
	    &+((3*A3*(E1*t1J22+E2*t2J22+E3*t3J22)+X*f2*(E1*t4J22+E2*t5J22+E3*t6J22))/(h11))

 
  	
end select
END DO
	
end subroutine

end program 
