gh1h1h1=-3*(L4*R11*R12^2*v1 - 2*cb*IL5*R11*R12*R13*v1 + cb^2*L4*R11*R13^2*v1 + 
  R11*R12^2*RL5*v1 - cb^2*R11*R13^2*RL5*v1 - IL5*R12^2*R13*sb*v1 + 
  cb^2*IL5*R13^3*sb*v1 - 2*cb*R12*R13^2*RL5*sb*v1 + 
  L1*R11*(R11^2 + R13^2*sb^2)*v1 + L4*R11^2*R12*v2 + L2*R12^3*v2 - 
  cb*IL5*R11^2*R13*v2 + cb^2*L2*R12*R13^2*v2 + R11^2*R12*RL5*v2 - 
  2*IL5*R11*R12*R13*sb*v2 - 2*cb*R11*R13^2*RL5*sb*v2 + L4*R12*R13^2*sb^2*v2 + 
  cb*IL5*R13^3*sb^2*v2 - R12*R13^2*RL5*sb^2*v2 + 
  L3*(R11*(R12^2 + cb^2*R13^2)*v1 + R11^2*R12*v2 + R12*R13^2*sb^2*v2))
