/PhysicalParameters/matlab/whip_eig.m

Large files files are truncated, but you can click here to view the full file

% File: whip_eig.m
% Date: February 20, 2008
% Author: Jason Moore
% Description: A function that takes the 25 benchmark parameters and
% returns the eigenvalues and eigenvectors of the system
function [S,D]=whip_eig(bench_par)
% Declare the benchmark parameters
%-------------------------------+--------------------------+---------------
% Quantity                      | Value                    | Units         
%-------------------------------|--------------------------|---------------
C                               =  bench_par(1);            % M             
G                               =  bench_par(2);            % N/KG          
IBXX                            =  bench_par(3);            % KG*M^2        
IBXZ                            =  bench_par(4);            % KG*M^2        
IBZZ                            =  bench_par(5);            % KG*M^2        
IFXX                            =  bench_par(6);            % KG*M^2        
IFYY                            =  bench_par(7);            % KG*M^2        
IHXX                            =  bench_par(8);            % KG*M^2       
IHXZ                            =  bench_par(9);            % KG*M^2       
IHZZ                            =  bench_par(10);           % KG*M^2       
IRXX                            =  bench_par(11);           % KG*M^2       
IRYY                            =  bench_par(12);           % KG*M^2       
LAMBDA                          =  bench_par(13);           % RAD          
M_B                             =  bench_par(14);           % KG           
M_F                             =  bench_par(15);           % KG           
M_H                             =  bench_par(16);           % KG           
M_R                             =  bench_par(17);           % KG           
RF                              =  bench_par(18);           % M            
RR                              =  bench_par(19);           % M            
V                               =  bench_par(20);           % M/S          
W                               =  bench_par(21);           % M            
XB                              =  bench_par(22);           % M            
XH                              =  bench_par(23);           % M            
ZB                              =  bench_par(24);           % M            
ZH                              =  bench_par(25);           % M            
%-------------------------------+--------------------------+---------------
% Reserve space and initialize matrices
z = zeros(519,1);

% Evaluate constants
D1 = cos(LAMBDA)*(C+W-RR*tan(LAMBDA));
D3 = -cos(LAMBDA)*(C-RF*tan(LAMBDA));
IC11 = IRXX;
IC22 = IRYY;
IC33 = IRXX;
ID11 = IBXX;
ID31 = IBXZ;
ID33 = IBZZ;
IF11 = IHZZ*sin(LAMBDA)^2 + cos(LAMBDA)*(IHXX*cos(LAMBDA)-2*IHXZ*sin(LAMBDA));
IF31 = sin(LAMBDA)*(IHXX*cos(LAMBDA)-IHXZ*sin(LAMBDA)) + cos(LAMBDA)*(IHXZ*cos(LAMBDA)-IHZZ*sin(LAMBDA));
IF33 = IHZZ*cos(LAMBDA)^2 + sin(LAMBDA)*(IHXX*sin(LAMBDA)+2*IHXZ*cos(LAMBDA));
IG11 = IFXX;
IG22 = IFYY;
IG33 = IFXX;
L1 = XB;
L2 = RR + ZB;
L3 = XH*cos(LAMBDA) - C*cos(LAMBDA) - W*cos(LAMBDA) - ZH*sin(LAMBDA);
L4 = RR*cos(LAMBDA) + XH*sin(LAMBDA) + ZH*cos(LAMBDA);
MC = M_R;
MD = M_B;
MF = M_H;
MG = M_F;
z(9) = cos(LAMBDA);
z(10) = sin(LAMBDA);
z(15) = (RF-RR-D1*z(10)-D3*z(10))/z(9);
z(516) = G*MC;
z(517) = G*MD;
z(518) = G*MF;
z(519) = G*MG;
% Calculate the A matrix in x'' = Ax + Bu, the state vector is
% [q4 q7 u4 u7]' and the input vector is u=[T4 T7]'
A(1,1) = 0;
A(1,2) = 0;
A(1,3) = 1;
A(1,4) = 0;
A(2,1) = 0;
A(2,2) = 0;
A(2,3) = 0;
A(2,4) = 1;
A(3,1) = -((L1*z(517)*(D3-RF*z(10))/(D1*z(9)+D3*z(9)-z(10)*z(15))-RF*z(10)*z(519)-z(518)*(L3-(D3-RF*z(10))*(D1*z(9)+L3*z(9)+L4*z(10))/(D1*z(9)+D3*z(9)-z(10)*z(15))))*(D3*MG*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))-IF31*z(9)-IF33*z(10)-IG33*z(10)-D3*MG*(RR+D1*z(10)+D3*z(10)+z(9)*z(15))-L3*MF*(RR+D1*z(10)+L3*z(10)-L4*z(9))-((RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(IF31*z(10)-IF33*z(9)-IG33*z(9)-L3*MF*(D1*z(9)+L3*z(9)+L4*z(10)))+(D3-RF*z(10))*(IF11*z(9)*z(10)+IF31*z(10)^2+IG11*z(9)*z(10)+MF*(D1*L4*z(9)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(L3*z(9)+L4*z(10))-D1*z(10)*(L4*z(10)+2*L3*z(9))-(L3*z(9)+L4*z(10))*(L3*z(10)-L4*z(9)))+MG*(D1*z(15)*z(10)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(D3*z(9)-z(10)*z(15))-D1*z(9)*(z(9)*z(15)+2*D3*z(10))-(D3*z(10)+z(9)*z(15))*(D3*z(9)-z(10)*z(15)))-ID31-IF31*z(9)^2-IF33*z(9)*z(10)-IG33*z(9)*z(10)-L1*MD*(RR-L2)-(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(2*IF31*z(9)*z(10)-IC33-ID33-IF11*z(10)^2-IF33*z(9)^2-IG11*z(10)^2-IG33*z(9)^2-MD*L1^2-MF*(D1^2*z(9)^2+(L3*z(9)+L4*z(10))^2+2*D1*z(9)*(L3*z(9)+L4*z(10)))-MG*(D1^2*z(9)^2+(D3*z(9)-z(10)*z(15))^2+2*D1*z(9)*(D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15)))+(2*D3*MG*(D3-RF*z(10))-IF33-IG33-MF*L3^2-MG*D3^2-(D3-RF*z(10))*(2*IF31*z(10)-2*IF33*z(9)-2*IG33*z(9)-2*L3*MF*(D1*z(9)+L3*z(9)+L4*z(10))-(D3-RF*z(10))*(2*IF31*z(9)*z(10)-IC33-ID33-IF11*z(10)^2-IF33*z(9)^2-IG11*z(10)^2-IG33*z(9)^2-MD*L1^2-MF*(D1^2*z(9)^2+(L3*z(9)+L4*z(10))^2+2*D1*z(9)*(L3*z(9)+L4*z(10)))-MG*(D1^2*z(9)^2+(D3*z(9)-z(10)*z(15))^2+2*D1*z(9)*(D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15)))/(D1*z(9)+D3*z(9)-z(10)*z(15)))*(RR*z(516)+z(519)*(RR+D3*z(10)+z(9)*z(15)+D1*(z(10)-z(9)*(RR+D1*z(10)+D3*z(10)+z(9)*z(15))/(D1*z(9)+D3*z(9)-z(10)*z(15)))+(RF*D1*z(9)+RF*D3*z(9)*(z(9)^2+z(10)^2)*(1-z(9)^2-z(10)^2)-RF*z(10)*z(15)*(z(9)^2+z(10)^2)*(1-z(9)^2-z(10)^2)-(D3*z(9)-z(10)*z(15))*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2)))/(D1*z(9)+D3*z(9)-z(10)*z(15)))-z(517)*(L2-RR-L1*(RF-RR-D1*z(10)-D3*z(10)-z(9)*z(15))/(D1*z(9)+D3*z(9)-z(10)*z(15)))-z(518)*(L4*z(9)-RR-L3*z(10)-D1*(z(10)-z(9)*(RR+D1*z(10)+D3*z(10)+z(9)*z(15))/(D1*z(9)+D3*z(9)-z(10)*z(15)))-(RF*D1*z(9)+RF*L3*z(9)*(z(9)^2+z(10)^2)*(1-z(9)^2-z(10)^2)-RF*L4*z(10)*(z(9)^2+z(10)^2)*(-1+z(9)^2+z(10)^2)-(L3*z(9)+L4*z(10))*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2)))/(D1*z(9)+D3*z(9)-z(10)*z(15)))))/((2*D3*MG*(D3-RF*z(10))-IF33-IG33-MF*L3^2-MG*D3^2-(D3-RF*z(10))*(2*IF31*z(10)-2*IF33*z(9)-2*IG33*z(9)-2*L3*MF*(D1*z(9)+L3*z(9)+L4*z(10))-(D3-RF*z(10))*(2*IF31*z(9)*z(10)-IC33-ID33-IF11*z(10)^2-IF33*z(9)^2-IG11*z(10)^2-IG33*z(9)^2-MD*L1^2-MF*(D1^2*z(9)^2+(L3*z(9)+L4*z(10))^2+2*D1*z(9)*(L3*z(9)+L4*z(10)))-MG*(D1^2*z(9)^2+(D3*z(9)-z(10)*z(15))^2+2*D1*z(9)*(D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15)))/(D1*z(9)+D3*z(9)-z(10)*z(15)))*(MD*(2*RR*L2-RR^2-L2^2)-IC11-ID11-2*IF31*z(9)*z(10)-RR^2*MC-IF11*z(9)^2-IF33*z(10)^2-IG11*z(9)^2-IG33*z(10)^2-MG*(RR^2+2*RR*D1*z(10)+D1^2*z(10)^2+(D3*z(10)+z(9)*z(15))^2+2*RR*(D3*z(10)+z(9)*z(15))+2*D1*z(10)*(D3*z(10)+z(9)*z(15)))-MF*(RR^2+2*RR*D1*z(10)+D1^2*z(10)^2+(L3*z(10)-L4*z(9))^2+2*RR*(L3*z(10)-L4*z(9))+2*D1*z(10)*(L3*z(10)-L4*z(9)))-(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(2*IF11*z(9)*z(10)+2*IF31*z(10)^2+2*IG11*z(9)*z(10)+2*MF*(D1*L4*z(9)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(L3*z(9)+L4*z(10))-D1*z(10)*(L4*z(10)+2*L3*z(9))-(L3*z(9)+L4*z(10))*(L3*z(10)-L4*z(9)))+2*MG*(D1*z(15)*z(10)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(D3*z(9)-z(10)*z(15))-D1*z(9)*(z(9)*z(15)+2*D3*z(10))-(D3*z(10)+z(9)*z(15))*(D3*z(9)-z(10)*z(15)))-2*ID31-2*IF31*z(9)^2-2*IF33*z(9)*z(10)-2*IG33*z(9)*z(10)-2*L1*MD*(RR-L2)-(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(2*IF31*z(9)*z(10)-IC33-ID33-IF11*z(10)^2-IF33*z(9)^2-IG11*z(10)^2-IG33*z(9)^2-MD*L1^2-MF*(D1^2*z(9)^2+(L3*z(9)+L4*z(10))^2+2*D1*z(9)*(L3*z(9)+L4*z(10)))-MG*(D1^2*z(9)^2+(D3*z(9)-z(10)*z(15))^2+2*D1*z(9)*(D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15)))/(D1*z(9)+D3*z(9)-z(10)*z(15)))-(D3*MG*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))-IF31*z(9)-IF33*z(10)-IG33*z(10)-D3*MG*(RR+D1*z(10)+D3*z(10)+z(9)*z(15))-L3*MF*(RR+D1*z(10)+L3*z(10)-L4*z(9))-((RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(IF31*z(10)-IF33*z(9)-IG33*z(9)-L3*MF*(D1*z(9)+L3*z(9)+L4*z(10))-(D3-RF*z(10))*(2*IF31*z(9)*z(10)-IC33-ID33-IF11*z(10)^2-IF33*z(9)^2-IG11*z(10)^2-IG33*z(9)^2-MD*L1^2-MF*(D1^2*z(9)^2+(L3*z(9)+L4*z(10))^2+2*D1*z(9)*(L3*z(9)+L4*z(10)))-MG*(D1^2*z(9)^2+(D3*z(9)-z(10)*z(15))^2+2*D1*z(9)*(D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15)))-(D3-RF*z(10))*(ID31+IF31*z(9)^2+IF33*z(9)*z(10)+IG33*z(9)*z(10)+L1*MD*(RR-L2)-IF11*z(9)*z(10)-IF31*z(10)^2-IG11*z(9)*z(10)-MF*(D1*L4*z(9)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(L3*z(9)+L4*z(10))-D1*z(10)*(L4*z(10)+2*L3*z(9))-(L3*z(9)+L4*z(10))*(L3*z(10)-L4*z(9)))-MG*(D1*z(15)*z(10)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(D3*z(9)-z(10)*z(15))-D1*z(9)*(z(9)*z(15)+2*D3*z(10))-(D3*z(10)+z(9)*z(15))*(D3*z(9)-z(10)*z(15)))))/(D1*z(9)+D3*z(9)-z(10)*z(15)))*(D3*MG*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))-IF31*z(9)-IF33*z(10)-IG33*z(10)-D3*MG*(RR+D1*z(10)+D3*z(10)+z(9)*z(15))-L3*MF*(RR+D1*z(10)+L3*z(10)-L4*z(9))-((RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(IF31*z(10)-IF33*z(9)-IG33*z(9)-L3*MF*(D1*z(9)+L3*z(9)+L4*z(10)))+(D3-RF*z(10))*(IF11*z(9)*z(10)+IF31*z(10)^2+IG11*z(9)*z(10)+MF*(D1*L4*z(9)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(L3*z(9)+L4*z(10))-D1*z(10)*(L4*z(10)+2*L3*z(9))-(L3*z(9)+L4*z(10))*(L3*z(10)-L4*z(9)))+MG*(D1*z(15)*z(10)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(D3*z(9)-z(10)*z(15))-D1*z(9)*(z(9)*z(15)+2*D3*z(10))-(D3*z(10)+z(9)*z(15))*(D3*z(9)-z(10)*z(15)))-ID31-IF31*z(9)^2-IF33*z(9)*z(10)-IG33*z(9)*z(10)-L1*MD*(RR-L2)-(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(2*IF31*z(9)*z(10)-IC33-ID33-IF11*z(10)^2-IF33*z(9)^2-IG11*z(10)^2-IG33*z(9)^2-MD*L1^2-MF*(D1^2*z(9)^2+(L3*z(9)+L4*z(10))^2+2*D1*z(9)*(L3*z(9)+L4*z(10)))-MG*(D1^2*z(9)^2+(D3*z(9)-z(10)*z(15))^2+2*D1*z(9)*(D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15))));
A(3,2) = ((2*D3*MG*(D3-RF*z(10))-IF33-IG33-MF*L3^2-MG*D3^2-(D3-RF*z(10))*(2*IF31*z(10)-2*IF33*z(9)-2*IG33*z(9)-2*L3*MF*(D1*z(9)+L3*z(9)+L4*z(10))-(D3-RF*z(10))*(2*IF31*z(9)*z(10)-IC33-ID33-IF11*z(10)^2-IF33*z(9)^2-IG11*z(10)^2-IG33*z(9)^2-MD*L1^2-MF*(D1^2*z(9)^2+(L3*z(9)+L4*z(10))^2+2*D1*z(9)*(L3*z(9)+L4*z(10)))-MG*(D1^2*z(9)^2+(D3*z(9)-z(10)*z(15))^2+2*D1*z(9)*(D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15)))/(D1*z(9)+D3*z(9)-z(10)*z(15)))*(L1*z(517)*(D3-RF*z(10))/(D1*z(9)+D3*z(9)-z(10)*z(15))+z(519)*((D1*z(9)*(D3-RF*z(10))+z(15)*z(10)^2*(RR+D1*z(10)+z(9)*z(15)))/(D1*z(9)+D3*z(9)-z(10)*z(15))-D3*z(9)*(z(9)-(D3*z(9)^2-z(10)*(RR+D1*z(10)+2*z(9)*z(15)))/(D1*z(9)+D3*z(9)-z(10)*z(15)))-z(10)*(D3*z(10)-(D3*z(9)-z(10)*z(15))*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))/(D1*z(9)+D3*z(9)-z(10)*z(15))))+z(518)*((D1*z(9)*(D3-RF*z(10))+L4*z(10)*(D3*z(9)^2-z(10)*(RR+D1*z(10)+z(9)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15))-L3*z(9)*(z(9)-(D3*z(9)^2-z(10)*(RR+D1*z(10)+z(9)*z(15)))/(D1*z(9)+D3*z(9)-z(10)*z(15)))-z(10)*(L3*z(10)-(L3*z(9)+L4*z(10))*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))/(D1*z(9)+D3*z(9)-z(10)*z(15))))+V^2*z(9)*((IC11+IC22-IC33)/(RR*(D1*z(9)+D3*z(9)-z(10)*z(15)))+(IG11*z(9)^2+IG33*z(10)^2+z(9)^2*(IG22-IG33)-z(10)^2*(IG11-IG22))/(RF*(D1*z(9)+D3*z(9)-z(10)*z(15)))+2*z(9)*z(10)*(IG11-IG33)*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))/(RF*(D1*z(9)+D3*z(9)-z(10)*z(15))^2)+(RR*MC+MD*(RR-L2)+MG*(RR+D1*z(10)+D3*z(10)+z(9)*z(15))+MF*(RR+D1*z(10)+L3*z(10)-L4*z(9))-MG*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))-(L1*MD+MF*(D1*z(9)+L3*z(9)+L4*z(10)))*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))/(D1*z(9)+D3*z(9)-z(10)*z(15)))/(D1*z(9)+D3*z(9)-z(10)*z(15))))+(V^2*z(9)*(z(10)*((IG11-IG22-IG33)/(RF*(D1*z(9)+D3*z(9)-z(10)*z(15)))-2*z(9)*(IG11-IG33)*(D3-RF*z(10))/(RF*(D1*z(9)+D3*z(9)-z(10)*z(15))^2))-(D3*MG+L3*MF-MG*(D3-RF*z(10))-(D3-RF*z(10))*(L1*MD+MF*(D1*z(9)+L3*z(9)+L4*z(10)))/(D1*z(9)+D3*z(9)-z(10)*z(15)))/(D1*z(9)+D3*z(9)-z(10)*z(15)))-z(10)*(L1*z(517)*(D3-RF*z(10))/(D1*z(9)+D3*z(9)-z(10)*z(15))-z(518)*(L3-(D3-RF*z(10))*(D1*z(9)+L3*z(9)+L4*z(10)+L4*z(10)*(-1+z(9)^2+z(10)^2)-L3*z(9)*(1-z(9)^2-z(10)^2))/(D1*z(9)+D3*z(9)-z(10)*z(15)))-z(519)*(D3+(D3-RF*z(10))*(z(10)*z(15)+D3*z(9)*(1-z(9)^2-z(10)^2)-D1*z(9)-D3*z(9)-z(10)*z(15)*(1-z(9)^2-z(10)^2))/(D1*z(9)+D3*z(9)-z(10)*z(15)))))*(D3*MG*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))-IF31*z(9)-IF33*z(10)-IG33*z(10)-D3*MG*(RR+D1*z(10)+D3*z(10)+z(9)*z(15))-L3*MF*(RR+D1*z(10)+L3*z(10)-L4*z(9))-((RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(IF31*z(10)-IF33*z(9)-IG33*z(9)-L3*MF*(D1*z(9)+L3*z(9)+L4*z(10)))+(D3-RF*z(10))*(IF11*z(9)*z(10)+IF31*z(10)^2+IG11*z(9)*z(10)+MF*(D1*L4*z(9)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(L3*z(9)+L4*z(10))-D1*z(10)*(L4*z(10)+2*L3*z(9))-(L3*z(9)+L4*z(10))*(L3*z(10)-L4*z(9)))+MG*(D1*z(15)*z(10)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(D3*z(9)-z(10)*z(15))-D1*z(9)*(z(9)*z(15)+2*D3*z(10))-(D3*z(10)+z(9)*z(15))*(D3*z(9)-z(10)*z(15)))-ID31-IF31*z(9)^2-IF33*z(9)*z(10)-IG33*z(9)*z(10)-L1*MD*(RR-L2)-(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(2*IF31*z(9)*z(10)-IC33-ID33-IF11*z(10)^2-IF33*z(9)^2-IG11*z(10)^2-IG33*z(9)^2-MD*L1^2-MF*(D1^2*z(9)^2+(L3*z(9)+L4*z(10))^2+2*D1*z(9)*(L3*z(9)+L4*z(10)))-MG*(D1^2*z(9)^2+(D3*z(9)-z(10)*z(15))^2+2*D1*z(9)*(D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15))))/((2*D3*MG*(D3-RF*z(10))-IF33-IG33-MF*L3^2-MG*D3^2-(D3-RF*z(10))*(2*IF31*z(10)-2*IF33*z(9)-2*IG33*z(9)-2*L3*MF*(D1*z(9)+L3*z(9)+L4*z(10))-(D3-RF*z(10))*(2*IF31*z(9)*z(10)-IC33-ID33-IF11*z(10)^2-IF33*z(9)^2-IG11*z(10)^2-IG33*z(9)^2-MD*L1^2-MF*(D1^2*z(9)^2+(L3*z(9)+L4*z(10))^2+2*D1*z(9)*(L3*z(9)+L4*z(10)))-MG*(D1^2*z(9)^2+(D3*z(9)-z(10)*z(15))^2+2*D1*z(9)*(D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15)))/(D1*z(9)+D3*z(9)-z(10)*z(15)))*(MD*(2*RR*L2-RR^2-L2^2)-IC11-ID11-2*IF31*z(9)*z(10)-RR^2*MC-IF11*z(9)^2-IF33*z(10)^2-IG11*z(9)^2-IG33*z(10)^2-MG*(RR^2+2*RR*D1*z(10)+D1^2*z(10)^2+(D3*z(10)+z(9)*z(15))^2+2*RR*(D3*z(10)+z(9)*z(15))+2*D1*z(10)*(D3*z(10)+z(9)*z(15)))-MF*(RR^2+2*RR*D1*z(10)+D1^2*z(10)^2+(L3*z(10)-L4*z(9))^2+2*RR*(L3*z(10)-L4*z(9))+2*D1*z(10)*(L3*z(10)-L4*z(9)))-(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(2*IF11*z(9)*z(10)+2*IF31*z(10)^2+2*IG11*z(9)*z(10)+2*MF*(D1*L4*z(9)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(L3*z(9)+L4*z(10))-D1*z(10)*(L4*z(10)+2*L3*z(9))-(L3*z(9)+L4*z(10))*(L3*z(10)-L4*z(9)))+2*MG*(D1*z(15)*z(10)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(D3*z(9)-z(10)*z(15))-D1*z(9)*(z(9)*z(15)+2*D3*z(10))-(D3*z(10)+z(9)*z(15))*(D3*z(9)-z(10)*z(15)))-2*ID31-2*IF31*z(9)^2-2*IF33*z(9)*z(10)-2*IG33*z(9)*z(10)-2*L1*MD*(RR-L2)-(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(2*IF31*z(9)*z(10)-IC33-ID33-IF11*z(10)^2-IF33*z(9)^2-IG11*z(10)^2-IG33*z(9)^2-MD*L1^2-MF*(D1^2*z(9)^2+(L3*z(9)+L4*z(10))^2+2*D1*z(9)*(L3*z(9)+L4*z(10)))-MG*(D1^2*z(9)^2+(D3*z(9)-z(10)*z(15))^2+2*D1*z(9)*(D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15)))/(D1*z(9)+D3*z(9)-z(10)*z(15)))-(D3*MG*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))-IF31*z(9)-IF33*z(10)-IG33*z(10)-D3*MG*(RR+D1*z(10)+D3*z(10)+z(9)*z(15))-L3*MF*(RR+D1*z(10)+L3*z(10)-L4*z(9))-((RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(IF31*z(10)-IF33*z(9)-IG33*z(9)-L3*MF*(D1*z(9)+L3*z(9)+L4*z(10))-(D3-RF*z(10))*(2*IF31*z(9)*z(10)-IC33-ID33-IF11*z(10)^2-IF33*z(9)^2-IG11*z(10)^2-IG33*z(9)^2-MD*L1^2-MF*(D1^2*z(9)^2+(L3*z(9)+L4*z(10))^2+2*D1*z(9)*(L3*z(9)+L4*z(10)))-MG*(D1^2*z(9)^2+(D3*z(9)-z(10)*z(15))^2+2*D1*z(9)*(D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15)))-(D3-RF*z(10))*(ID31+IF31*z(9)^2+IF33*z(9)*z(10)+IG33*z(9)*z(10)+L1*MD*(RR-L2)-IF11*z(9)*z(10)-IF31*z(10)^2-IG11*z(9)*z(10)-MF*(D1*L4*z(9)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(L3*z(9)+L4*z(10))-D1*z(10)*(L4*z(10)+2*L3*z(9))-(L3*z(9)+L4*z(10))*(L3*z(10)-L4*z(9)))-MG*(D1*z(15)*z(10)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(D3*z(9)-z(10)*z(15))-D1*z(9)*(z(9)*z(15)+2*D3*z(10))-(D3*z(10)+z(9)*z(15))*(D3*z(9)-z(10)*z(15)))))/(D1*z(9)+D3*z(9)-z(10)*z(15)))*(D3*MG*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))-IF31*z(9)-IF33*z(10)-IG33*z(10)-D3*MG*(RR+D1*z(10)+D3*z(10)+z(9)*z(15))-L3*MF*(RR+D1*z(10)+L3*z(10)-L4*z(9))-((RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(IF31*z(10)-IF33*z(9)-IG33*z(9)-L3*MF*(D1*z(9)+L3*z(9)+L4*z(10)))+(D3-RF*z(10))*(IF11*z(9)*z(10)+IF31*z(10)^2+IG11*z(9)*z(10)+MF*(D1*L4*z(9)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(L3*z(9)+L4*z(10))-D1*z(10)*(L4*z(10)+2*L3*z(9))-(L3*z(9)+L4*z(10))*(L3*z(10)-L4*z(9)))+MG*(D1*z(15)*z(10)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(D3*z(9)-z(10)*z(15))-D1*z(9)*(z(9)*z(15)+2*D3*z(10))-(D3*z(10)+z(9)*z(15))*(D3*z(9)-z(10)*z(15)))-ID31-IF31*z(9)^2-IF33*z(9)*z(10)-IG33*z(9)*z(10)-L1*MD*(RR-L2)-(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(2*IF31*z(9)*z(10)-IC33-ID33-IF11*z(10)^2-IF33*z(9)^2-IG11*z(10)^2-IG33*z(9)^2-MD*L1^2-MF*(D1^2*z(9)^2+(L3*z(9)+L4*z(10))^2+2*D1*z(9)*(L3*z(9)+L4*z(10)))-MG*(D1^2*z(9)^2+(D3*z(9)-z(10)*z(15))^2+2*D1*z(9)*(D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15))));
A(3,3) = -V*((2*D3*MG*(D3-RF*z(10))-IF33-IG33-MF*L3^2-MG*D3^2-(D3-RF*z(10))*(2*IF31*z(10)-2*IF33*z(9)-2*IG33*z(9)-2*L3*MF*(D1*z(9)+L3*z(9)+L4*z(10))-(D3-RF*z(10))*(2*IF31*z(9)*z(10)-IC33-ID33-IF11*z(10)^2-IF33*z(9)^2-IG11*z(10)^2-IG33*z(9)^2-MD*L1^2-MF*(D1^2*z(9)^2+(L3*z(9)+L4*z(10))^2+2*D1*z(9)*(L3*z(9)+L4*z(10)))-MG*(D1^2*z(9)^2+(D3*z(9)-z(10)*z(15))^2+2*D1*z(9)*(D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15)))/(D1*z(9)+D3*z(9)-z(10)*z(15)))*(2*IC11*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))/(RR*(D1*z(9)+D3*z(9)-z(10)*z(15)))+(IG33*z(10)*(z(9)+z(10)*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))/(D1*z(9)+D3*z(9)-z(10)*z(15)))-IG11*z(9)*(z(10)-z(9)*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))/(D1*z(9)+D3*z(9)-z(10)*z(15)))-z(10)*(IG11-IG22)*(z(9)+z(10)*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))/(D1*z(9)+D3*z(9)-z(10)*z(15)))-z(9)*(IG22-IG33)*(z(10)-z(9)*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))/(D1*z(9)+D3*z(9)-z(10)*z(15))))/RF-(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(2*IC33/RR+MG*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))+(L1*MD+MF*(D1*z(9)+L3*z(9)+L4*z(10)))*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))/(D1*z(9)+D3*z(9)-z(10)*z(15))-RR*MC-MD*(RR-L2)-MG*(RR+D1*z(10)+D3*z(10)+z(9)*z(15))-MF*(RR+D1*z(10)+L3*z(10)-L4*z(9))-(z(9)*(IG11-IG22)*(z(9)+z(10)*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))/(D1*z(9)+D3*z(9)-z(10)*z(15)))-IG33*z(9)*(z(9)+z(10)*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))/(D1*z(9)+D3*z(9)-z(10)*z(15)))-IG11*z(10)*(z(10)-z(9)*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))/(D1*z(9)+D3*z(9)-z(10)*z(15)))-z(10)*(IG22-IG33)*(z(10)-z(9)*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))/(D1*z(9)+D3*z(9)-z(10)*z(15))))/RF)/(D1*z(9)+D3*z(9)-z(10)*z(15)))+((IG11-IG22-IG33)*(z(9)+z(10)*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))/(D1*z(9)+D3*z(9)-z(10)*z(15)))/RF-((RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(D3*MG+L3*MF-MG*(D3-RF*z(10))-(D3-RF*z(10))*(L1*MD+MF*(D1*z(9)+L3*z(9)+L4*z(10)))/(D1*z(9)+D3*z(9)-z(10)*z(15)))+(D3-RF*z(10))*((IC11-IC22-IC33)/RR+(z(9)*(IG11-IG22)*(z(9)+z(10)*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))/(D1*z(9)+D3*z(9)-z(10)*z(15)))-IG33*z(9)*(z(9)+z(10)*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))/(D1*z(9)+D3*z(9)-z(10)*z(15)))-IG11*z(10)*(z(10)-z(9)*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))/(D1*z(9)+D3*z(9)-z(10)*z(15)))-z(10)*(IG22-IG33)*(z(10)-z(9)*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))/(D1*z(9)+D3*z(9)-z(10)*z(15))))/RF))/(D1*z(9)+D3*z(9)-z(10)*z(15)))*(D3*MG*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))-IF31*z(9)-IF33*z(10)-IG33*z(10)-D3*MG*(RR+D1*z(10)+D3*z(10)+z(9)*z(15))-L3*MF*(RR+D1*z(10)+L3*z(10)-L4*z(9))-((RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(IF31*z(10)-IF33*z(9)-IG33*z(9)-L3*MF*(D1*z(9)+L3*z(9)+L4*z(10)))+(D3-RF*z(10))*(IF11*z(9)*z(10)+IF31*z(10)^2+IG11*z(9)*z(10)+MF*(D1*L4*z(9)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(L3*z(9)+L4*z(10))-D1*z(10)*(L4*z(10)+2*L3*z(9))-(L3*z(9)+L4*z(10))*(L3*z(10)-L4*z(9)))+MG*(D1*z(15)*z(10)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(D3*z(9)-z(10)*z(15))-D1*z(9)*(z(9)*z(15)+2*D3*z(10))-(D3*z(10)+z(9)*z(15))*(D3*z(9)-z(10)*z(15)))-ID31-IF31*z(9)^2-IF33*z(9)*z(10)-IG33*z(9)*z(10)-L1*MD*(RR-L2)-(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(2*IF31*z(9)*z(10)-IC33-ID33-IF11*z(10)^2-IF33*z(9)^2-IG11*z(10)^2-IG33*z(9)^2-MD*L1^2-MF*(D1^2*z(9)^2+(L3*z(9)+L4*z(10))^2+2*D1*z(9)*(L3*z(9)+L4*z(10)))-MG*(D1^2*z(9)^2+(D3*z(9)-z(10)*z(15))^2+2*D1*z(9)*(D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15))))/((2*D3*MG*(D3-RF*z(10))-IF33-IG33-MF*L3^2-MG*D3^2-(D3-RF*z(10))*(2*IF31*z(10)-2*IF33*z(9)-2*IG33*z(9)-2*L3*MF*(D1*z(9)+L3*z(9)+L4*z(10))-(D3-RF*z(10))*(2*IF31*z(9)*z(10)-IC33-ID33-IF11*z(10)^2-IF33*z(9)^2-IG11*z(10)^2-IG33*z(9)^2-MD*L1^2-MF*(D1^2*z(9)^2+(L3*z(9)+L4*z(10))^2+2*D1*z(9)*(L3*z(9)+L4*z(10)))-MG*(D1^2*z(9)^2+(D3*z(9)-z(10)*z(15))^2+2*D1*z(9)*(D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15)))/(D1*z(9)+D3*z(9)-z(10)*z(15)))*(MD*(2*RR*L2-RR^2-L2^2)-IC11-ID11-2*IF31*z(9)*z(10)-RR^2*MC-IF11*z(9)^2-IF33*z(10)^2-IG11*z(9)^2-IG33*z(10)^2-MG*(RR^2+2*RR*D1*z(10)+D1^2*z(10)^2+(D3*z(10)+z(9)*z(15))^2+2*RR*(D3*z(10)+z(9)*z(15))+2*D1*z(10)*(D3*z(10)+z(9)*z(15)))-MF*(RR^2+2*RR*D1*z(10)+D1^2*z(10)^2+(L3*z(10)-L4*z(9))^2+2*RR*(L3*z(10)-L4*z(9))+2*D1*z(10)*(L3*z(10)-L4*z(9)))-(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(2*IF11*z(9)*z(10)+2*IF31*z(10)^2+2*IG11*z(9)*z(10)+2*MF*(D1*L4*z(9)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(L3*z(9)+L4*z(10))-D1*z(10)*(L4*z(10)+2*L3*z(9))-(L3*z(9)+L4*z(10))*(L3*z(10)-L4*z(9)))+2*MG*(D1*z(15)*z(10)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(D3*z(9)-z(10)*z(15))-D1*z(9)*(z(9)*z(15)+2*D3*z(10))-(D3*z(10)+z(9)*z(15))*(D3*z(9)-z(10)*z(15)))-2*ID31-2*IF31*z(9)^2-2*IF33*z(9)*z(10)-2*IG33*z(9)*z(10)-2*L1*MD*(RR-L2)-(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(2*IF31*z(9)*z(10)-IC33-ID33-IF11*z(10)^2-IF33*z(9)^2-IG11*z(10)^2-IG33*z(9)^2-MD*L1^2-MF*(D1^2*z(9)^2+(L3*z(9)+L4*z(10))^2+2*D1*z(9)*(L3*z(9)+L4*z(10)))-MG*(D1^2*z(9)^2+(D3*z(9)-z(10)*z(15))^2+2*D1*z(9)*(D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15)))/(D1*z(9)+D3*z(9)-z(10)*z(15)))-(D3*MG*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))-IF31*z(9)-IF33*z(10)-IG33*z(10)-D3*MG*(RR+D1*z(10)+D3*z(10)+z(9)*z(15))-L3*MF*(RR+D1*z(10)+L3*z(10)-L4*z(9))-((RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(IF31*z(10)-IF33*z(9)-IG33*z(9)-L3*MF*(D1*z(9)+L3*z(9)+L4*z(10))-(D3-RF*z(10))*(2*IF31*z(9)*z(10)-IC33-ID33-IF11*z(10)^2-IF33*z(9)^2-IG11*z(10)^2-IG33*z(9)^2-MD*L1^2-MF*(D1^2*z(9)^2+(L3*z(9)+L4*z(10))^2+2*D1*z(9)*(L3*z(9)+L4*z(10)))-MG*(D1^2*z(9)^2+(D3*z(9)-z(10)*z(15))^2+2*D1*z(9)*(D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15)))-(D3-RF*z(10))*(ID31+IF31*z(9)^2+IF33*z(9)*z(10)+IG33*z(9)*z(10)+L1*MD*(RR-L2)-IF11*z(9)*z(10)-IF31*z(10)^2-IG11*z(9)*z(10)-MF*(D1*L4*z(9)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(L3*z(9)+L4*z(10))-D1*z(10)*(L4*z(10)+2*L3*z(9))-(L3*z(9)+L4*z(10))*(L3*z(10)-L4*z(9)))-MG*(D1*z(15)*z(10)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(D3*z(9)-z(10)*z(15))-D1*z(9)*(z(9)*z(15)+2*D3*z(10))-(D3*z(10)+z(9)*z(15))*(D3*z(9)-z(10)*z(15)))))/(D1*z(9)+D3*z(9)-z(10)*z(15)))*(D3*MG*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))-IF31*z(9)-IF33*z(10)-IG33*z(10)-D3*MG*(RR+D1*z(10)+D3*z(10)+z(9)*z(15))-L3*MF*(RR+D1*z(10)+L3*z(10)-L4*z(9))-((RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(IF31*z(10)-IF33*z(9)-IG33*z(9)-L3*MF*(D1*z(9)+L3*z(9)+L4*z(10)))+(D3-RF*z(10))*(IF11*z(9)*z(10)+IF31*z(10)^2+IG11*z(9)*z(10)+MF*(D1*L4*z(9)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(L3*z(9)+L4*z(10))-D1*z(10)*(L4*z(10)+2*L3*z(9))-(L3*z(9)+L4*z(10))*(L3*z(10)-L4*z(9)))+MG*(D1*z(15)*z(10)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(D3*z(9)-z(10)*z(15))-D1*z(9)*(z(9)*z(15)+2*D3*z(10))-(D3*z(10)+z(9)*z(15))*(D3*z(9)-z(10)*z(15)))-ID31-IF31*z(9)^2-IF33*z(9)*z(10)-IG33*z(9)*z(10)-L1*MD*(RR-L2)-(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(2*IF31*z(9)*z(10)-IC33-ID33-IF11*z(10)^2-IF33*z(9)^2-IG11*z(10)^2-IG33*z(9)^2-MD*L1^2-MF*(D1^2*z(9)^2+(L3*z(9)+L4*z(10))^2+2*D1*z(9)*(L3*z(9)+L4*z(10)))-MG*(D1^2*z(9)^2+(D3*z(9)-z(10)*z(15))^2+2*D1*z(9)*(D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15))));
A(3,4) = -V*((D3*MG*z(9)+z(10)*(D3-RF*z(10))*((IG11-IG22-IG33)/(RF*(D1*z(9)+D3*z(9)-z(10)*z(15)))-(z(9)*(D3-RF*z(10))*(IG11-IG22-IG33)/(RF*(D1*z(9)+D3*z(9)-z(10)*z(15)))-(IG11+IG22-IG33)*(1-z(9)*(D3-RF*z(10))/(D1*z(9)+D3*z(9)-z(10)*z(15)))/RF)/(D1*z(9)+D3*z(9)-z(10)*z(15)))-((D3-RF*z(10))*(D3*MG+L3*MF-MG*(D3-RF*z(10))-(D3-RF*z(10))*(L1*MD+MF*(D1*z(9)+L3*z(9)+L4*z(10)))/(D1*z(9)+D3*z(9)-z(10)*z(15)))+z(9)*(IF31*z(10)-IF33*z(9)-IG33*z(9)-L3*MF*(D1*z(9)+L3*z(9)+L4*z(10))-(D3-RF*z(10))*(2*IF31*z(9)*z(10)-IC33-ID33-IF11*z(10)^2-IF33*z(9)^2-IG11*z(10)^2-IG33*z(9)^2-MD*L1^2-MF*(D1^2*z(9)^2+(L3*z(9)+L4*z(10))^2+2*D1*z(9)*(L3*z(9)+L4*z(10)))-MG*(D1^2*z(9)^2+(D3*z(9)-z(10)*z(15))^2+2*D1*z(9)*(D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15)))*(D3*MG*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))-IF31*z(9)-IF33*z(10)-IG33*z(10)-D3*MG*(RR+D1*z(10)+D3*z(10)+z(9)*z(15))-L3*MF*(RR+D1*z(10)+L3*z(10)-L4*z(9))-((RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(IF31*z(10)-IF33*z(9)-IG33*z(9)-L3*MF*(D1*z(9)+L3*z(9)+L4*z(10)))+(D3-RF*z(10))*(IF11*z(9)*z(10)+IF31*z(10)^2+IG11*z(9)*z(10)+MF*(D1*L4*z(9)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(L3*z(9)+L4*z(10))-D1*z(10)*(L4*z(10)+2*L3*z(9))-(L3*z(9)+L4*z(10))*(L3*z(10)-L4*z(9)))+MG*(D1*z(15)*z(10)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(D3*z(9)-z(10)*z(15))-D1*z(9)*(z(9)*z(15)+2*D3*z(10))-(D3*z(10)+z(9)*z(15))*(D3*z(9)-z(10)*z(15)))-ID31-IF31*z(9)^2-IF33*z(9)*z(10)-IG33*z(9)*z(10)-L1*MD*(RR-L2)-(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(2*IF31*z(9)*z(10)-IC33-ID33-IF11*z(10)^2-IF33*z(9)^2-IG11*z(10)^2-IG33*z(9)^2-MD*L1^2-MF*(D1^2*z(9)^2+(L3*z(9)+L4*z(10))^2+2*D1*z(9)*(L3*z(9)+L4*z(10)))-MG*(D1^2*z(9)^2+(D3*z(9)-z(10)*z(15))^2+2*D1*z(9)*(D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15)))-(2*D3*MG*(D3-RF*z(10))-IF33-IG33-MF*L3^2-MG*D3^2-(D3-RF*z(10))*(2*IF31*z(10)-2*IF33*z(9)-2*IG33*z(9)-2*L3*MF*(D1*z(9)+L3*z(9)+L4*z(10))-(D3-RF*z(10))*(2*IF31*z(9)*z(10)-IC33-ID33-IF11*z(10)^2-IF33*z(9)^2-IG11*z(10)^2-IG33*z(9)^2-MD*L1^2-MF*(D1^2*z(9)^2+(L3*z(9)+L4*z(10))^2+2*D1*z(9)*(L3*z(9)+L4*z(10)))-MG*(D1^2*z(9)^2+(D3*z(9)-z(10)*z(15))^2+2*D1*z(9)*(D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15)))/(D1*z(9)+D3*z(9)-z(10)*z(15)))*(z(10)^2*(D3-RF*z(10))*(IG11-IG22-IG33)/(RF*(D1*z(9)+D3*z(9)-z(10)*z(15)))+z(9)*(IG11+IG22-IG33)*(1-z(9)*(D3-RF*z(10))/(D1*z(9)+D3*z(9)-z(10)*z(15)))/RF-(D3-RF*z(10))*(IC11+IC22-IC33)/(RR*(D1*z(9)+D3*z(9)-z(10)*z(15)))-(z(10)*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(z(9)*(D3-RF*z(10))*(IG11-IG22-IG33)/(RF*(D1*z(9)+D3*z(9)-z(10)*z(15)))-(IG11+IG22-IG33)*(1-z(9)*(D3-RF*z(10))/(D1*z(9)+D3*z(9)-z(10)*z(15)))/RF)+(D3-RF*z(10))*(RR*MC+MD*(RR-L2)+MG*(RR+D1*z(10)+D3*z(10)+z(9)*z(15))+MF*(RR+D1*z(10)+L3*z(10)-L4*z(9))-MG*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))-(L1*MD+MF*(D1*z(9)+L3*z(9)+L4*z(10)))*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))/(D1*z(9)+D3*z(9)-z(10)*z(15)))+z(9)*(IF11*z(9)*z(10)+IF31*z(10)^2+IG11*z(9)*z(10)+MF*(D1*L4*z(9)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(L3*z(9)+L4*z(10))-D1*z(10)*(L4*z(10)+2*L3*z(9))-(L3*z(9)+L4*z(10))*(L3*z(10)-L4*z(9)))+MG*(D1*z(15)*z(10)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(D3*z(9)-z(10)*z(15))-D1*z(9)*(z(9)*z(15)+2*D3*z(10))-(D3*z(10)+z(9)*z(15))*(D3*z(9)-z(10)*z(15)))-ID31-IF31*z(9)^2-IF33*z(9)*z(10)-IG33*z(9)*z(10)-L1*MD*(RR-L2)-(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(2*IF31*z(9)*z(10)-IC33-ID33-IF11*z(10)^2-IF33*z(9)^2-IG11*z(10)^2-IG33*z(9)^2-MD*L1^2-MF*(D1^2*z(9)^2+(L3*z(9)+L4*z(10))^2+2*D1*z(9)*(L3*z(9)+L4*z(10)))-MG*(D1^2*z(9)^2+(D3*z(9)-z(10)*z(15))^2+2*D1*z(9)*(D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15))))/((2*D3*MG*(D3-RF*z(10))-IF33-IG33-MF*L3^2-MG*D3^2-(D3-RF*z(10))*(2*IF31*z(10)-2*IF33*z(9)-2*IG33*z(9)-2*L3*MF*(D1*z(9)+L3*z(9)+L4*z(10))-(D3-RF*z(10))*(2*IF31*z(9)*z(10)-IC33-ID33-IF11*z(10)^2-IF33*z(9)^2-IG11*z(10)^2-IG33*z(9)^2-MD*L1^2-MF*(D1^2*z(9)^2+(L3*z(9)+L4*z(10))^2+2*D1*z(9)*(L3*z(9)+L4*z(10)))-MG*(D1^2*z(9)^2+(D3*z(9)-z(10)*z(15))^2+2*D1*z(9)*(D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15)))/(D1*z(9)+D3*z(9)-z(10)*z(15)))*(MD*(2*RR*L2-RR^2-L2^2)-IC11-ID11-2*IF31*z(9)*z(10)-RR^2*MC-IF11*z(9)^2-IF33*z(10)^2-IG11*z(9)^2-IG33*z(10)^2-MG*(RR^2+2*RR*D1*z(10)+D1^2*z(10)^2+(D3*z(10)+z(9)*z(15))^2+2*RR*(D3*z(10)+z(9)*z(15))+2*D1*z(10)*(D3*z(10)+z(9)*z(15)))-MF*(RR^2+2*RR*D1*z(10)+D1^2*z(10)^2+(L3*z(10)-L4*z(9))^2+2*RR*(L3*z(10)-L4*z(9))+2*D1*z(10)*(L3*z(10)-L4*z(9)))-(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(2*IF11*z(9)*z(10)+2*IF31*z(10)^2+2*IG11*z(9)*z(10)+2*MF*(D1*L4*z(9)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(L3*z(9)+L4*z(10))-D1*z(10)*(L4*z(10)+2*L3*z(9))-(L3*z(9)+L4*z(10))*(L3*z(10)-L4*z(9)))+2*MG*(D1*z(15)*z(10)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(D3*z(9)-z(10)*z(15))-D1*z(9)*(z(9)*z(15)+2*D3*z(10))-(D3*z(10)+z(9)*z(15))*(D3*z(9)-z(10)*z(15)))-2*ID31-2*IF31*z(9)^2-2*IF33*z(9)*z(10)-2*IG33*z(9)*z(10)-2*L1*MD*(RR-L2)-(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(2*IF31*z(9)*z(10)-IC33-ID33-IF11*z(10)^2-IF33*z(9)^2-IG11*z(10)^2-IG33*z(9)^2-MD*L1^2-MF*(D1^2*z(9)^2+(L3*z(9)+L4*z(10))^2+2*D1*z(9)*(L3*z(9)+L4*z(10)))-MG*(D1^2*z(9)^2+(D3*z(9)-z(10)*z(15))^2+2*D1*z(9)*(D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15)))/(D1*z(9)+D3*z(9)-z(10)*z(15)))-(D3*MG*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))-IF31*z(9)-IF33*z(10)-IG33*z(10)-D3*MG*(RR+D1*z(10)+D3*z(10)+z(9)*z(15))-L3*MF*(RR+D1*z(10)+L3*z(10)-L4*z(9))-((RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(IF31*z(10)-IF33*z(9)-IG33*z(9)-L3*MF*(D1*z(9)+L3*z(9)+L4*z(10))-(D3-RF*z(10))*(2*IF31*z(9)*z(10)-IC33-ID33-IF11*z(10)^2-IF33*z(9)^2-IG11*z(10)^2-IG33*z(9)^2-MD*L1^2-MF*(D1^2*z(9)^2+(L3*z(9)+L4*z(10))^2+2*D1*z(9)*(L3*z(9)+L4*z(10)))-MG*(D1^2*z(9)^2+(D3*z(9)-z(10)*z(15))^2+2*D1*z(9)*(D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15)))-(D3-RF*z(10))*(ID31+IF31*z(9)^2+IF33*z(9)*z(10)+IG33*z(9)*z(10)+L1*MD*(RR-L2)-IF11*z(9)*z(10)-IF31*z(10)^2-IG11*z(9)*z(10)-MF*(D1*L4*z(9)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(L3*z(9)+L4*z(10))-D1*z(10)*(L4*z(10)+2*L3*z(9))-(L3*z(9)+L4*z(10))*(L3*z(10)-L4*z(9)))-MG*(D1*z(15)*z(10)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(D3*z(9)-z(10)*z(15))-D1*z(9)*(z(9)*z(15)+2*D3*z(10))-(D3*z(10)+z(9)*z(15))*(D3*z(9)-z(10)*z(15)))))/(D1*z(9)+D3*z(9)-z(10)*z(15)))*(D3*MG*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))-IF31*z(9)-IF33*z(10)-IG33*z(10)-D3*MG*(RR+D1*z(10)+D3*z(10)+z(9)*z(15))-L3*MF*(RR+D1*z(10)+L3*z(10)-L4*z(9))-((RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(IF31*z(10)-IF33*z(9)-IG33*z(9)-L3*MF*(D1*z(9)+L3*z(9)+L4*z(10)))+(D3-RF*z(10))*(IF11*z(9)*z(10)+IF31*z(10)^2+IG11*z(9)*z(10)+MF*(D1*L4*z(9)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(L3*z(9)+L4*z(10))-D1*z(10)*(L4*z(10)+2*L3*z(9))-(L3*z(9)+L4*z(10))*(L3*z(10)-L4*z(9)))+MG*(D1*z(15)*z(10)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(D3*z(9)-z(10)*z(15))-D1*z(9)*(z(9)*z(15)+2*D3*z(10))-(D3*z(10)+z(9)*z(15))*(D3*z(9)-z(10)*z(15)))-ID31-IF31*z(9)^2-IF33*z(9)*z(10)-IG33*z(9)*z(10)-L1*MD*(RR-L2)-(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(2*IF31*z(9)*z(10)-IC33-ID33-IF11*z(10)^2-IF33*z(9)^2-IG11*z(10)^2-IG33*z(9)^2-MD*L1^2-MF*(D1^2*z(9)^2+(L3*z(9)+L4*z(10))^2+2*D1*z(9)*(L3*z(9)+L4*z(10)))-MG*(D1^2*z(9)^2+(D3*z(9)-z(10)*z(15))^2+2*D1*z(9)*(D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15))));
A(4,1) = ((L1*z(517)*(D3-RF*z(10))/(D1*z(9)+D3*z(9)-z(10)*z(15))-RF*z(10)*z(519)-z(518)*(L3-(D3-RF*z(10))*(D1*z(9)+L3*z(9)+L4*z(10))/(D1*z(9)+D3*z(9)-z(10)*z(15))))*(MD*(2*RR*L2-RR^2-L2^2)-IC11-ID11-2*IF31*z(9)*z(10)-RR^2*MC-IF11*z(9)^2-IF33*z(10)^2-IG11*z(9)^2-IG33*z(10)^2-MG*(RR^2+2*RR*D1*z(10)+D1^2*z(10)^2+(D3*z(10)+z(9)*z(15))^2+2*RR*(D3*z(10)+z(9)*z(15))+2*D1*z(10)*(D3*z(10)+z(9)*z(15)))-MF*(RR^2+2*RR*D1*z(10)+D1^2*z(10)^2+(L3*z(10)-L4*z(9))^2+2*RR*(L3*z(10)-L4*z(9))+2*D1*z(10)*(L3*z(10)-L4*z(9)))-(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(2*IF11*z(9)*z(10)+2*IF31*z(10)^2+2*IG11*z(9)*z(10)+2*MF*(D1*L4*z(9)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(L3*z(9)+L4*z(10))-D1*z(10)*(L4*z(10)+2*L3*z(9))-(L3*z(9)+L4*z(10))*(L3*z(10)-L4*z(9)))+2*MG*(D1*z(15)*z(10)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(D3*z(9)-z(10)*z(15))-D1*z(9)*(z(9)*z(15)+2*D3*z(10))-(D3*z(10)+z(9)*z(15))*(D3*z(9)-z(10)*z(15)))-2*ID31-2*IF31*z(9)^2-2*IF33*z(9)*z(10)-2*IG33*z(9)*z(10)-2*L1*MD*(RR-L2)-(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(2*IF31*z(9)*z(10)-IC33-ID33-IF11*z(10)^2-IF33*z(9)^2-IG11*z(10)^2-IG33*z(9)^2-MD*L1^2-MF*(D1^2*z(9)^2+(L3*z(9)+L4*z(10))^2+2*D1*z(9)*(L3*z(9)+L4*z(10)))-MG*(D1^2*z(9)^2+(D3*z(9)-z(10)*z(15))^2+2*D1*z(9)*(D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15)))/(D1*z(9)+D3*z(9)-z(10)*z(15)))+(RR*z(516)+z(519)*(RR+D3*z(10)+z(9)*z(15)+D1*(z(10)-z(9)*(RR+D1*z(10)+D3*z(10)+z(9)*z(15))/(D1*z(9)+D3*z(9)-z(10)*z(15)))+(RF*D1*z(9)+RF*D3*z(9)*(z(9)^2+z(10)^2)*(1-z(9)^2-z(10)^2)-RF*z(10)*z(15)*(z(9)^2+z(10)^2)*(1-z(9)^2-z(10)^2)-(D3*z(9)-z(10)*z(15))*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2)))/(D1*z(9)+D3*z(9)-z(10)*z(15)))-z(517)*(L2-RR-L1*(RF-RR-D1*z(10)-D3*z(10)-z(9)*z(15))/(D1*z(9)+D3*z(9)-z(10)*z(15)))-z(518)*(L4*z(9)-RR-L3*z(10)-D1*(z(10)-z(9)*(RR+D1*z(10)+D3*z(10)+z(9)*z(15))/(D1*z(9)+D3*z(9)-z(10)*z(15)))-(RF*D1*z(9)+RF*L3*z(9)*(z(9)^2+z(10)^2)*(1-z(9)^2-z(10)^2)-RF*L4*z(10)*(z(9)^2+z(10)^2)*(-1+z(9)^2+z(10)^2)-(L3*z(9)+L4*z(10))*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2)))/(D1*z(9)+D3*z(9)-z(10)*z(15))))*(D3*MG*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))-IF31*z(9)-IF33*z(10)-IG33*z(10)-D3*MG*(RR+D1*z(10)+D3*z(10)+z(9)*z(15))-L3*MF*(RR+D1*z(10)+L3*z(10)-L4*z(9))-((RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(IF31*z(10)-IF33*z(9)-IG33*z(9)-L3*MF*(D1*z(9)+L3*z(9)+L4*z(10))-(D3-RF*z(10))*(2*IF31*z(9)*z(10)-IC33-ID33-IF11*z(10)^2-IF33*z(9)^2-IG11*z(10)^2-IG33*z(9)^2-MD*L1^2-MF*(D1^2*z(9)^2+(L3*z(9)+L4*z(10))^2+2*D1*z(9)*(L3*z(9)+L4*z(10)))-MG*(D1^2*z(9)^2+(D3*z(9)-z(10)*z(15))^2+2*D1*z(9)*(D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15)))-(D3-RF*z(10))*(ID31+IF31*z(9)^2+IF33*z(9)*z(10)+IG33*z(9)*z(10)+L1*MD*(RR-L2)-IF11*z(9)*z(10)-IF31*z(10)^2-IG11*z(9)*z(10)-MF*(D1*L4*z(9)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(L3*z(9)+L4*z(10))-D1*z(10)*(L4*z(10)+2*L3*z(9))-(L3*z(9)+L4*z(10))*(L3*z(10)-L4*z(9)))-MG*(D1*z(15)*z(10)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(D3*z(9)-z(10)*z(15))-D1*z(9)*(z(9)*z(15)+2*D3*z(10))-(D3*z(10)+z(9)*z(15))*(D3*z(9)-z(10)*z(15)))))/(D1*z(9)+D3*z(9)-z(10)*z(15))))/((2*D3*MG*(D3-RF*z(10))-IF33-IG33-MF*L3^2-MG*D3^2-(D3-RF*z(10))*(2*IF31*z(10)-2*IF33*z(9)-2*IG33*z(9)-2*L3*MF*(D1*z(9)+L3*z(9)+L4*z(10))-(D3-RF*z(10))*(2*IF31*z(9)*z(10)-IC33-ID33-IF11*z(10)^2-IF33*z(9)^2-IG11*z(10)^2-IG33*z(9)^2-MD*L1^2-MF*(D1^2*z(9)^2+(L3*z(9)+L4*z(10))^2+2*D1*z(9)*(L3*z(9)+L4*z(10)))-MG*(D1^2*z(9)^2+(D3*z(9)-z(10)*z(15))^2+2*D1*z(9)*(D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15)))/(D1*z(9)+D3*z(9)-z(10)*z(15)))*(MD*(2*RR*L2-RR^2-L2^2)-IC11-ID11-2*IF31*z(9)*z(10)-RR^2*MC-IF11*z(9)^2-IF33*z(10)^2-IG11*z(9)^2-IG33*z(10)^2-MG*(RR^2+2*RR*D1*z(10)+D1^2*z(10)^2+(D3*z(10)+z(9)*z(15))^2+2*RR*(D3*z(10)+z(9)*z(15))+2*D1*z(10)*(D3*z(10)+z(9)*z(15)))-MF*(RR^2+2*RR*D1*z(10)+D1^2*z(10)^2+(L3*z(10)-L4*z(9))^2+2*RR*(L3*z(10)-L4*z(9))+2*D1*z(10)*(L3*z(10)-L4*z(9)))-(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(2*IF11*z(9)*z(10)+2*IF31*z(10)^2+2*IG11*z(9)*z(10)+2*MF*(D1*L4*z(9)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(L3*z(9)+L4*z(10))-D1*z(10)*(L4*z(10)+2*L3*z(9))-(L3*z(9)+L4*z(10))*(L3*z(10)-L4*z(9)))+2*MG*(D1*z(15)*z(10)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(D3*z(9)-z(10)*z(15))-D1*z(9)*(z(9)*z(15)+2*D3*z(10))-(D3*z(10)+z(9)*z(15))*(D3*z(9)-z(10)*z(15)))-2*ID31-2*IF31*z(9)^2-2*IF33*z(9)*z(10)-2*IG33*z(9)*z(10)-2*L1*MD*(RR-L2)-(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(2*IF31*z(9)*z(10)-IC33-ID33-IF11*z(10)^2-IF33*z(9)^2-IG11*z(10)^2-IG33*z(9)^2-MD*L1^2-MF*(D1^2*z(9)^2+(L3*z(9)+L4*z(10))^2+2*D1*z(9)*(L3*z(9)+L4*z(10)))-MG*(D1^2*z(9)^2+(D3*z(9)-z(10)*z(15))^2+2*D1*z(9)*(D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15)))/(D1*z(9)+D3*z(9)-z(10)*z(15)))-(D3*MG*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))-IF31*z(9)-IF33*z(10)-IG33*z(10)-D3*MG*(RR+D1*z(10)+D3*z(10)+z(9)*z(15))-L3*MF*(RR+D1*z(10)+L3*z(10)-L4*z(9))-((RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(IF31*z(10)-IF33*z(9)-IG33*z(9)-L3*MF*(D1*z(9)+L3*z(9)+L4*z(10))-(D3-RF*z(10))*(2*IF31*z(9)*z(10)-IC33-ID33-IF11*z(10)^2-IF33*z(9)^2-IG11*z(10)^2-IG33*z(9)^2-MD*L1^2-MF*(D1^2*z(9)^2+(L3*z(9)+L4*z(10))^2+2*D1*z(9)*(L3*z(9)+L4*z(10)))-MG*(D1^2*z(9)^2+(D3*z(9)-z(10)*z(15))^2+2*D1*z(9)*(D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15)))-(D3-RF*z(10))*(ID31+IF31*z(9)^2+IF33*z(9)*z(10)+IG33*z(9)*z(10)+L1*MD*(RR-L2)-IF11*z(9)*z(10)-IF31*z(10)^2-IG11*z(9)*z(10)-MF*(D1*L4*z(9)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(L3*z(9)+L4*z(10))-D1*z(10)*(L4*z(10)+2*L3*z(9))-(L3*z(9)+L4*z(10))*(L3*z(10)-L4*z(9)))-MG*(D1*z(15)*z(10)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(D3*z(9)-z(10)*z(15))-D1*z(9)*(z(9)*z(15)+2*D3*z(10))-(D3*z(10)+z(9)*z(15))*(D3*z(9)-z(10)*z(15)))))/(D1*z(9)+D3*z(9)-z(10)*z(15)))*(D3*MG*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))-IF31*z(9)-IF33*z(10)-IG33*z(10)-D3*MG*(RR+D1*z(10)+D3*z(10)+z(9)*z(15))-L3*MF*(RR+D1*z(10)+L3*z(10)-L4*z(9))-((RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(IF31*z(10)-IF33*z(9)-IG33*z(9)-L3*MF*(D1*z(9)+L3*z(9)+L4*z(10)))+(D3-RF*z(10))*(IF11*z(9)*z(10)+IF31*z(10)^2+IG11*z(9)*z(10)+MF*(D1*L4*z(9)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(L3*z(9)+L4*z(10))-D1*z(10)*(L4*z(10)+2*L3*z(9))-(L3*z(9)+L4*z(10))*(L3*z(10)-L4*z(9)))+MG*(D1*z(15)*z(10)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(D3*z(9)-z(10)*z(15))-D1*z(9)*(z(9)*z(15)+2*D3*z(10))-(D3*z(10)+z(9)*z(15))*(D3*z(9)-z(10)*z(15)))-ID31-IF31*z(9)^2-IF33*z(9)*z(10)-IG33*z(9)*z(10)-L1*MD*(RR-L2)-(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(2*IF31*z(9)*z(10)-IC33-ID33-IF11*z(10)^2-IF33*z(9)^2-IG11*z(10)^2-IG33*z(9)^2-MD*L1^2-MF*(D1^2*z(9)^2+(L3*z(9)+L4*z(10))^2+2*D1*z(9)*(L3*z(9)+L4*z(10)))-MG*(D1^2*z(9)^2+(D3*z(9)-z(10)*z(15))^2+2*D1*z(9)*(D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15))));
A(4,2) = -((V^2*z(9)*(z(10)*((IG11-IG22-IG33)/(RF*(D1*z(9)+D3*z(9)-z(10)*z(15)))-2*z(9)*(IG11-IG33)*(D3-RF*z(10))/(RF*(D1*z(9)+D3*z(9)-z(10)*z(15))^2))-(D3*MG+L3*MF-MG*(D3-RF*z(10))-(D3-RF*z(10))*(L1*MD+MF*(D1*z(9)+L3*z(9)+L4*z(10)))/(D1*z(9)+D3*z(9)-z(10)*z(15)))/(D1*z(9)+D3*z(9)-z(10)*z(15)))-z(10)*(L1*z(517)*(D3-RF*z(10))/(D1*z(9)+D3*z(9)-z(10)*z(15))-z(518)*(L3-(D3-RF*z(10))*(D1*z(9)+L3*z(9)+L4*z(10)+L4*z(10)*(-1+z(9)^2+z(10)^2)-L3*z(9)*(1-z(9)^2-z(10)^2))/(D1*z(9)+D3*z(9)-z(10)*z(15)))-z(519)*(D3+(D3-RF*z(10))*(z(10)*z(15)+D3*z(9)*(1-z(9)^2-z(10)^2)-D1*z(9)-D3*z(9)-z(10)*z(15)*(1-z(9)^2-z(10)^2))/(D1*z(9)+D3*z(9)-z(10)*z(15)))))*(MD*(2*RR*L2-RR^2-L2^2)-IC11-ID11-2*IF31*z(9)*z(10)-RR^2*MC-IF11*z(9)^2-IF33*z(10)^2-IG11*z(9)^2-IG33*z(10)^2-MG*(RR^2+2*RR*D1*z(10)+D1^2*z(10)^2+(D3*z(10)+z(9)*z(15))^2+2*RR*(D3*z(10)+z(9)*z(15))+2*D1*z(10)*(D3*z(10)+z(9)*z(15)))-MF*(RR^2+2*RR*D1*z(10)+D1^2*z(10)^2+(L3*z(10)-L4*z(9))^2+2*RR*(L3*z(10)-L4*z(9))+2*D1*z(10)*(L3*z(10)-L4*z(9)))-(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(2*IF11*z(9)*z(10)+2*IF31*z(10)^2+2*IG11*z(9)*z(10)+2*MF*(D1*L4*z(9)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(L3*z(9)+L4*z(10))-D1*z(10)*(L4*z(10)+2*L3*z(9))-(L3*z(9)+L4*z(10))*(L3*z(10)-L4*z(9)))+2*MG*(D1*z(15)*z(10)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(D3*z(9)-z(10)*z(15))-D1*z(9)*(z(9)*z(15)+2*D3*z(10))-(D3*z(10)+z(9)*z(15))*(D3*z(9)-z(10)*z(15)))-2*ID31-2*IF31*z(9)^2-2*IF33*z(9)*z(10)-2*IG33*z(9)*z(10)-2*L1*MD*(RR-L2)-(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(2*IF31*z(9)*z(10)-IC33-ID33-IF11*z(10)^2-IF33*z(9)^2-IG11*z(10)^2-IG33*z(9)^2-MD*L1^2-MF*(D1^2*z(9)^2+(L3*z(9)+L4*z(10))^2+2*D1*z(9)*(L3*z(9)+L4*z(10)))-MG*(D1^2*z(9)^2+(D3*z(9)-z(10)*z(15))^2+2*D1*z(9)*(D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15)))/(D1*z(9)+D3*z(9)-z(10)*z(15)))+(D3*MG*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))-IF31*z(9)-IF33*z(10)-IG33*z(10)-D3*MG*(RR+D1*z(10)+D3*z(10)+z(9)*z(15))-L3*MF*(RR+D1*z(10)+L3*z(10)-L4*z(9))-((RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(IF31*z(10)-IF33*z(9)-IG33*z(9)-L3*MF*(D1*z(9)+L3*z(9)+L4*z(10))-(D3-RF*z(10))*(2*IF31*z(9)*z(10)-IC33-ID33-IF11*z(10)^2-IF33*z(9)^2-IG11*z(10)^2-IG33*z(9)^2-MD*L1^2-MF*(D1^2*z(9)^2+(L3*z(9)+L4*z(10))^2+2*D1*z(9)*(L3*z(9)+L4*z(10)))-MG*(D1^2*z(9)^2+(D3*z(9)-z(10)*z(15))^2+2*D1*z(9)*(D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15)))-(D3-RF*z(10))*(ID31+IF31*z(9)^2+IF33*z(9)*z(10)+IG33*z(9)*z(10)+L1*MD*(RR-L2)-IF11*z(9)*z(10)-IF31*z(10)^2-IG11*z(9)*z(10)-MF*(D1*L4*z(9)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(L3*z(9)+L4*z(10))-D1*z(10)*(L4*z(10)+2*L3*z(9))-(L3*z(9)+L4*z(10))*(L3*z(10)-L4*z(9)))-MG*(D1*z(15)*z(10)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(D3*z(9)-z(10)*z(15))-D1*z(9)*(z(9)*z(15)+2*D3*z(10))-(D3*z(10)+z(9)*z(15))*(D3*z(9)-z(10)*z(15)))))/(D1*z(9)+D3*z(9)-z(10)*z(15)))*(L1*z(517)*(D3-RF*z(10))/(D1*z(9)+D3*z(9)-z(10)*z(15))+z(519)*((D1*z(9)*(D3-RF*z(10))+z(15)*z(10)^2*(RR+D1*z(10)+z(9)*z(15)))/(D1*z(9)+D3*z(9)-z(10)*z(15))-D3*z(9)*(z(9)-(D3*z(9)^2-z(10)*(RR+D1*z(10)+2*z(9)*z(15)))/(D1*z(9)+D3*z(9)-z(10)*z(15)))-z(10)*(D3*z(10)-(D3*z(9)-z(10)*z(15))*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))/(D1*z(9)+D3*z(9)-z(10)*z(15))))+z(518)*((D1*z(9)*(D3-RF*z(10))+L4*z(10)*(D3*z(9)^2-z(10)*(RR+D1*z(10)+z(9)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15))-L3*z(9)*(z(9)-(D3*z(9)^2-z(10)*(RR+D1*z(10)+z(9)*z(15)))/(D1*z(9)+D3*z(9)-z(10)*z(15)))-z(10)*(L3*z(10)-(L3*z(9)+L4*z(10))*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))/(D1*z(9)+D3*z(9)-z(10)*z(15))))+V^2*z(9)*((IC11+IC22-IC33)/(RR*(D1*z(9)+D3*z(9)-z(10)*z(15)))+(IG11*z(9)^2+IG33*z(10)^2+z(9)^2*(IG22-IG33)-z(10)^2*(IG11-IG22))/(RF*(D1*z(9)+D3*z(9)-z(10)*z(15)))+2*z(9)*z(10)*(IG11-IG33)*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))/(RF*(D1*z(9)+D3*z(9)-z(10)*z(15))^2)+(RR*MC+MD*(RR-L2)+MG*(RR+D1*z(10)+D3*z(10)+z(9)*z(15))+MF*(RR+D1*z(10)+L3*z(10)-L4*z(9))-MG*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))-(L1*MD+MF*(D1*z(9)+L3*z(9)+L4*z(10)))*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))/(D1*z(9)+D3*z(9)-z(10)*z(15)))/(D1*z(9)+D3*z(9)-z(10)*z(15)))))/((2*D3*MG*(D3-RF*z(10))-IF33-IG33-MF*L3^2-MG*D3^2-(D3-RF*z(10))*(2*IF31*z(10)-2*IF33*z(9)-2*IG33*z(9)-2*L3*MF*(D1*z(9)+L3*z(9)+L4*z(10))-(D3-RF*z(10))*(2*IF31*z(9)*z(10)-IC33-ID33-IF11*z(10)^2-IF33*z(9)^2-IG11*z(10)^2-IG33*z(9)^2-MD*L1^2-MF*(D1^2*z(9)^2+(L3*z(9)+L4*z(10))^2+2*D1*z(9)*(L3*z(9)+L4*z(10)))-MG*(D1^2*z(9)^2+(D3*z(9)-z(10)*z(15))^2+2*D1*z(9)*(D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15)))/(D1*z(9)+D3*z(9)-z(10)*z(15)))*(MD*(2*RR*L2-RR^2-L2^2)-IC11-ID11-2*IF31*z(9)*z(10)-RR^2*MC-IF11*z(9)^2-IF33*z(10)^2-IG11*z(9)^2-IG33*z(10)^2-MG*(RR^2+2*RR*D1*z(10)+D1^2*z(10)^2+(D3*z(10)+z(9)*z(15))^2+2*RR*(D3*z(10)+z(9)*z(15))+2*D1*z(10)*(D3*z(10)+z(9)*z(15)))-MF*(RR^2+2*RR*D1*z(10)+D1^2*z(10)^2+(L3*z(10)-L4*z(9))^2+2*RR*(L3*z(10)-L4*z(9))+2*D1*z(10)*(L3*z(10)-L4*z(9)))-(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(2*IF11*z(9)*z(10)+2*IF31*z(10)^2+2*IG11*z(9)*z(10)+2*MF*(D1*L4*z(9)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(L3*z(9)+L4*z(10))-D1*z(10)*(L4*z(10)+2*L3*z(9))-(L3*z(9)+L4*z(10))*(L3*z(10)-L4*z(9)))+2*MG*(D1*z(15)*z(10)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(D3*z(9)-z(10)*z(15))-D1*z(9)*(z(9)*z(15)+2*D3*z(10))-(D3*z(10)+z(9)*z(15))*(D3*z(9)-z(10)*z(15)))-2*ID31-2*IF31*z(9)^2-2*IF33*z(9)*z(10)-2*IG33*z(9)*z(10)-2*L1*MD*(RR-L2)-(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(2*IF31*z(9)*z(10)-IC33-ID33-IF11*z(10)^2-IF33*z(9)^2-IG11*z(10)^2-IG33*z(9)^2-MD*L1^2-MF*(D1^2*z(9)^2+(L3*z(9)+L4*z(10))^2+2*D1*z(9)*(L3*z(9)+L4*z(10)))-MG*(D1^2*z(9)^2+(D3*z(9)-z(10)*z(15))^2+2*D1*z(9)*(D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15)))/(D1*z(9)+D3*z(9)-z(10)*z(15)))-(D3*MG*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))-IF31*z(9)-IF33*z(10)-IG33*z(10)-D3*MG*(RR+D1*z(10)+D3*z(10)+z(9)*z(15))-L3*MF*(RR+D1*z(10)+L3*z(10)-L4*z(9))-((RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(IF31*z(10)-IF33*z(9)-IG33*z(9)-L3*MF*(D1*z(9)+L3*z(9)+L4*z(10))-(D3-RF*z(10))*(2*IF31*z(9)*z(10)-IC33-ID33-IF11*z(10)^2-IF33*z(9)^2-IG11*z(10)^2-IG33*z(9)^2-MD*L1^2-MF*(D1^2*z(9)^2+(L3*z(9)+L4*z(10))^2+2*D1*z(9)*(L3*z(9)+L4*z(10)))-MG*(D1^2*z(9)^2+(D3*z(9)-z(10)*z(15))^2+2*D1*z(9)*(D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15)))-(D3-RF*z(10))*(ID31+IF31*z(9)^2+IF33*z(9)*z(10)+IG33*z(9)*z(10)+L1*MD*(RR-L2)-IF11*z(9)*z(10)-IF31*z(10)^2-IG11*z(9)*z(10)-MF*(D1*L4*z(9)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(L3*z(9)+L4*z(10))-D1*z(10)*(L4*z(10)+2*L3*z(9))-(L3*z(9)+L4*z(10))*(L3*z(10)-L4*z(9)))-MG*(D1*z(15)*z(10)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(D3*z(9)-z(10)*z(15))-D1*z(9)*(z(9)*z(15)+2*D3*z(10))-(D3*z(10)+z(9)*z(15))*(D3*z(9)-z(10)*z(15)))))/(D1*z(9)+D3*z(9)-z(10)*z(15)))*(D3*MG*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))-IF31*z(9)-IF33*z(10)-IG33*z(10)-D3*MG*(RR+D1*z(10)+D3*z(10)+z(9)*z(15))-L3*MF*(RR+D1*z(10)+L3*z(10)-L4*z(9))-((RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(IF31*z(10)-IF33*z(9)-IG33*z(9)-L3*MF*(D1*z(9)+L3*z(9)+L4*z(10)))+(D3-RF*z(10))*(IF11*z(9)*z(10)+IF31*z(10)^2+IG11*z(9)*z(10)+MF*(D1*L4*z(9)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(L3*z(9)+L4*z(10))-D1*z(10)*(L4*z(10)+2*L3*z(9))-(L3*z(9)+L4*z(10))*(L3*z(10)-L4*z(9)))+MG*(D1*z(15)*z(10)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(D3*z(9)-z(10)*z(15))-D1*z(9)*(z(9)*z(15)+2*D3*z(10))-(D3*z(10)+z(9)*z(15))*(D3*z(9)-z(10)*z(15)))-ID31-IF31*z(9)^2-IF33*z(9)*z(10)-IG33*z(9)*z(10)-L1*MD*(RR-L2)-(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(2*IF31*z(9)*z(10)-IC33-ID33-IF11*z(10)^2-IF33*z(9)^2-IG11*z(10)^2-IG33*z(9)^2-MD*L1^2-MF*(D1^2*z(9)^2+(L3*z(9)+L4*z(10))^2+2*D1*z(9)*(L3*z(9)+L4*z(10)))-MG*(D1^2*z(9)^2+(D3*z(9)-z(10)*z(15))^2+2*D1*z(9)*(D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15))));
A(4,3) = V*(((IG11-IG22-IG33)*(z(9)+z(10)*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))/(D1*z(9)+D3*z(9)-z(10)*z(15)))/RF-((RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(D3*MG+L3*MF-MG*(D3-RF*z(10))-(D3-RF*z(10))*(L1*MD+MF*(D1*z(9)+L3*z(9)+L4*z(10)))/(D1*z(9)+D3*z(9)-z(10)*z(15)))+(D3-RF*z(10))*((IC11-IC22-IC33)/RR+(z(9)*(IG11-IG22)*(z(9)+z(10)*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))/(D1*z(9)+D3*z(9)-z(10)*z(15)))-IG33*z(9)*(z(9)+z(10)*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))/(D1*z(9)+D3*z(9)-z(10)*z(15)))-IG11*z(10)*(z(10)-z(9)*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))/(D1*z(9)+D3*z(9)-z(10)*z(15)))-z(10)*(IG22-IG33)*(z(10)-z(9)*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))/(D1*z(9)+D3*z(9)-z(10)*z(15))))/RF))/(D1*z(9)+D3*z(9)-z(10)*z(15)))*(MD*(2*RR*L2-RR^2-L2^2)-IC11-ID11-2*IF31*z(9)*z(10)-RR^2*MC-IF11*z(9)^2-IF33*z(10)^2-IG11*z(9)^2-IG33*z(10)^2-MG*(RR^2+2*RR*D1*z(10)+D1^2*z(10)^2+(D3*z(10)+z(9)*z(15))^2+2*RR*(D3*z(10)+z(9)*z(15))+2*D1*z(10)*(D3*z(10)+z(9)*z(15)))-MF*(RR^2+2*RR*D1*z(10)+D1^2*z(10)^2+(L3*z(10)-L4*z(9))^2+2*RR*(L3*z(10)-L4*z(9))+2*D1*z(10)*(L3*z(10)-L4*z(9)))-(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(2*IF11*z(9)*z(10)+2*IF31*z(10)^2+2*IG11*z(9)*z(10)+2*MF*(D1*L4*z(9)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(L3*z(9)+L4*z(10))-D1*z(10)*(L4*z(10)+2*L3*z(9))-(L3*z(9)+L4*z(10))*(L3*z(10)-L4*z(9)))+2*MG*(D1*z(15)*z(10)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(D3*z(9)-z(10)*z(15))-D1*z(9)*(z(9)*z(15)+2*D3*z(10))-(D3*z(10)+z(9)*z(15))*(D3*z(9)-z(10)*z(15)))-2*ID31-2*IF31*z(9)^2-2*IF33*z(9)*z(10)-2*IG33*z(9)*z(10)-2*L1*MD*(RR-L2)-(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(2*IF31*z(9)*z(10)-IC33-ID33-IF11*z(10)^2-IF33*z(9)^2-IG11*z(10)^2-IG33*z(9)^2-MD*L1^2-MF*(D1^2*z(9)^2+(L3*z(9)+L4*z(10))^2+2*D1*z(9)*(L3*z(9)+L4*z(10)))-MG*(D1^2*z(9)^2+(D3*z(9)-z(10)*z(15))^2+2*D1*z(9)*(D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15)))/(D1*z(9)+D3*z(9)-z(10)*z(15)))+(D3*MG*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))-IF31*z(9)-IF33*z(10)-IG33*z(10)-D3*MG*(RR+D1*z(10)+D3*z(10)+z(9)*z(15))-L3*MF*(RR+D1*z(10)+L3*z(10)-L4*z(9))-((RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(IF31*z(10)-IF33*z(9)-IG33*z(9)-L3*MF*(D1*z(9)+L3*z(9)+L4*z(10))-(D3-RF*z(10))*(2*IF31*z(9)*z(10)-IC33-ID33-IF11*z(10)^2-IF33*z(9)^2-IG11*z(10)^2-IG33*z(9)^2-MD*L1^2-MF*(D1^2*z(9)^2+(L3*z(9)+L4*z(10))^2+2*D1*z(9)*(L3*z(9)+L4*z(10)))-MG*(D1^2*z(9)^2+(D3*z(9)-z(10)*z(15))^2+2*D1*z(9)*(D3*z(9)-z(10)*z(15))))/(D1*z(9)+D3*z(9)-z(10)*z(15)))-(D3-RF*z(10))*(ID31+IF31*z(9)^2+IF33*z(9)*z(10)+IG33*z(9)*z(10)+L1*MD*(RR-L2)-IF11*z(9)*z(10)-IF31*z(10)^2-IG11*z(9)*z(10)-MF*(D1*L4*z(9)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(L3*z(9)+L4*z(10))-D1*z(10)*(L4*z(10)+2*L3*z(9))-(L3*z(9)+L4*z(10))*(L3*z(10)-L4*z(9)))-MG*(D1*z(15)*z(10)^2-RR*D1*z(9)-z(9)*z(10)*D1^2-RR*(D3*z(9)-z(10)*z(15))-D1*z(9)*(z(9)*z(15)+2*D3*z(10))-(D3*z(10)+z(9)*z(15))*(D3*z(9)-z(10)*z(15)))))/(D1*z(9)+D3*z(9)-z(10)*z(15)))*(2*IC11*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))/(RR*(D1*z(9)+D3*z(9)-z(10)*z(15)))+(IG33*z(10)*(z(9)+z(10)*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))/(D1*z(9)+D3*z(9)-z(10)*z(15)))-IG11*z(9)*(z(10)-z(9)*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))/(D1*z(9)+D3*z(9)-z(10)*z(15)))-z(10)*(IG11-IG22)*(z(9)+z(10)*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))/(D1*z(9)+D3*z(9)-z(10)*z(15)))-z(9)*(IG22-IG33)*(z(10)-z(9)*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))/(D1*z(9)+D3*z(9)-z(10)*z(15))))/RF-(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))*(2*IC33/RR+MG*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))+(L1*MD+MF*(D1*z(9)+L3*z(9)+L4*z(10)))*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))/(D1*z(9)+D3*z(9)-z(10)*z(15))-RR*MC-MD*(RR-L2)-MG*(RR+D1*z(10)+D3*z(10)+z(9)*z(15))-MF*(RR+D1*z(10)+L3*z(10)-L4*z(9))-(z(9)*(IG11-IG22)*(z(9)+z(10)*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))/(D1*z(9)+D3*z(9)-z(10)*z(15)))-IG33*z(9)*(z(9)+z(10)*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))/(D1*z(9)+D3*z(9)-z(10)*z(15)))-IG11*z(10)*(z(10)-z(9)*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9)^2+z(10)^2))/(D1*z(9)+D3*z(9)-z(10)*z(15)))-z(10)*(IG22-IG33)*(z(10)-z(9)*(RR+D1*z(10)+D3*z(10)+z(9)*z(15)-RF*(z(9…