List(List((a*b*c1*kt*pAu*pBu*qAu*qBu)/4.,0,0),
  List(0,(a*b*c1*kt*pAv*pBv*qAv*qBv)/4.,0),
  List(0,0,(a*b*c1*kt*(pAw*pBw*qAw*qBw + 
        (4*kr*pAwxi*pBwxi*qAw*qBw)/(a**2*kt) + 
        (4*kr*pAw*pBw*qAweta*qBweta)/(bb**2*kt)))/4.))
