Shat7(1,z6(1,1,1,1,1,1),ep)=
     &  - Csum0(ep)
     &  - 5.D0*Csum1(ep)
     &  - 5.D0*Csum2(ep)
     &  - 10.D0*Csum11(ep)
     &  - 20.D0*Csum12(ep)
     &  - 10.D0*Csum22(ep)
     &  - 10.D0*Csum111(ep)
     &  - 10.D0*Csum222(ep)
     &  - 30.D0*Csum112(ep)
     &  - 30.D0*Csum122(ep)
     &  - 5.D0*Csum1111(ep)
     &  - 5.D0*Csum2222(ep)
     &  - 20.D0*Csum1112(ep)
     &  - 30.D0*Csum1122(ep)
     &  - 20.D0*Csum1222(ep)
     &  - Csum11111(ep)
     &  - Csum22222(ep)
     &  - 5.D0*Csum11112(ep)
     &  - 5.D0*Csum12222(ep)
     &  - 10.D0*Csum11122(ep)
     &  - 10.D0*Csum11222(ep)
     &

      Shat7(2,z6(1,1,1,1,1,1),ep)=
     &  - Csum0(ep)
     &  - 5.D0*Csum1(ep)
     &  - 5.D0*Csum2(ep)
     &  - 10.D0*Csum11(ep)
     &  - 20.D0*Csum12(ep)
     &  - 10.D0*Csum22(ep)
     &  - 10.D0*Csum111(ep)
     &  - 10.D0*Csum222(ep)
     &  - 30.D0*Csum112(ep)
     &  - 30.D0*Csum122(ep)
     &  - 5.D0*Csum1111(ep)
     &  - 5.D0*Csum2222(ep)
     &  - 20.D0*Csum1112(ep)
     &  - 30.D0*Csum1122(ep)
     &  - 20.D0*Csum1222(ep)
     &  - Csum11111(ep)
     &  - Csum22222(ep)
     &  - 5.D0*Csum11112(ep)
     &  - 5.D0*Csum12222(ep)
     &  - 10.D0*Csum11122(ep)
     &  - 10.D0*Csum11222(ep)
     &  + Cv(cc111111 + C124,ep)
     &

      Shat7(3,z6(1,1,1,1,1,1),ep)=
     &  - Csum0(ep)
     &  - 5.D0*Csum1(ep)
     &  - 5.D0*Csum2(ep)
     &  - 10.D0*Csum11(ep)
     &  - 20.D0*Csum12(ep)
     &  - 10.D0*Csum22(ep)
     &  - 10.D0*Csum111(ep)
     &  - 10.D0*Csum222(ep)
     &  - 30.D0*Csum112(ep)
     &  - 30.D0*Csum122(ep)
     &  - 5.D0*Csum1111(ep)
     &  - 5.D0*Csum2222(ep)
     &  - 20.D0*Csum1112(ep)
     &  - 30.D0*Csum1122(ep)
     &  - 20.D0*Csum1222(ep)
     &  - Csum11111(ep)
     &  - Csum22222(ep)
     &  - 5.D0*Csum11112(ep)
     &  - 5.D0*Csum12222(ep)
     &  - 10.D0*Csum11122(ep)
     &  - 10.D0*Csum11222(ep)
     &  + Cv(cc111111 + C123,ep)
     &

      Shat7(1,z6(1,1,1,1,1,2),ep)=
     &  + Csum1(ep)
     &  + 4.D0*Csum11(ep)
     &  + 4.D0*Csum12(ep)
     &  + 6.D0*Csum111(ep)
     &  + 12.D0*Csum112(ep)
     &  + 6.D0*Csum122(ep)
     &  + 4.D0*Csum1111(ep)
     &  + 12.D0*Csum1112(ep)
     &  + 12.D0*Csum1122(ep)
     &  + 4.D0*Csum1222(ep)
     &  + Csum11111(ep)
     &  + 4.D0*Csum11112(ep)
     &  + Csum12222(ep)
     &  + 6.D0*Csum11122(ep)
     &  + 4.D0*Csum11222(ep)
     &

      Shat7(2,z6(1,1,1,1,1,2),ep)=
     &  + Csum1(ep)
     &  + 4.D0*Csum11(ep)
     &  + 4.D0*Csum12(ep)
     &  + 6.D0*Csum111(ep)
     &  + 12.D0*Csum112(ep)
     &  + 6.D0*Csum122(ep)
     &  + 4.D0*Csum1111(ep)
     &  + 12.D0*Csum1112(ep)
     &  + 12.D0*Csum1122(ep)
     &  + 4.D0*Csum1222(ep)
     &  + Csum11111(ep)
     &  + 4.D0*Csum11112(ep)
     &  + Csum12222(ep)
     &  + 6.D0*Csum11122(ep)
     &  + 4.D0*Csum11222(ep)
     &

      Shat7(3,z6(1,1,1,1,1,2),ep)=
     &  + Csum1(ep)
     &  + 4.D0*Csum11(ep)
     &  + 4.D0*Csum12(ep)
     &  + 6.D0*Csum111(ep)
     &  + 12.D0*Csum112(ep)
     &  + 6.D0*Csum122(ep)
     &  + 4.D0*Csum1111(ep)
     &  + 12.D0*Csum1112(ep)
     &  + 12.D0*Csum1122(ep)
     &  + 4.D0*Csum1222(ep)
     &  + Csum11111(ep)
     &  + 4.D0*Csum11112(ep)
     &  + Csum12222(ep)
     &  + 6.D0*Csum11122(ep)
     &  + 4.D0*Csum11222(ep)
     &  + Cv(cc111112 + C123,ep)
     &

      Shat7(1,z6(1,1,1,1,1,3),ep)=
     &  + Csum2(ep)
     &  + 4.D0*Csum12(ep)
     &  + 4.D0*Csum22(ep)
     &  + 6.D0*Csum222(ep)
     &  + 6.D0*Csum112(ep)
     &  + 12.D0*Csum122(ep)
     &  + 4.D0*Csum2222(ep)
     &  + 4.D0*Csum1112(ep)
     &  + 12.D0*Csum1122(ep)
     &  + 12.D0*Csum1222(ep)
     &  + Csum22222(ep)
     &  + Csum11112(ep)
     &  + 4.D0*Csum12222(ep)
     &  + 4.D0*Csum11122(ep)
     &  + 6.D0*Csum11222(ep)
     &

      Shat7(2,z6(1,1,1,1,1,3),ep)=
     &  + Csum2(ep)
     &  + 4.D0*Csum12(ep)
     &  + 4.D0*Csum22(ep)
     &  + 6.D0*Csum222(ep)
     &  + 6.D0*Csum112(ep)
     &  + 12.D0*Csum122(ep)
     &  + 4.D0*Csum2222(ep)
     &  + 4.D0*Csum1112(ep)
     &  + 12.D0*Csum1122(ep)
     &  + 12.D0*Csum1222(ep)
     &  + Csum22222(ep)
     &  + Csum11112(ep)
     &  + 4.D0*Csum12222(ep)
     &  + 4.D0*Csum11122(ep)
     &  + 6.D0*Csum11222(ep)
     &  + Cv(cc111112 + C124,ep)
     &

      Shat7(3,z6(1,1,1,1,1,3),ep)=
     &  + Csum2(ep)
     &  + 4.D0*Csum12(ep)
     &  + 4.D0*Csum22(ep)
     &  + 6.D0*Csum222(ep)
     &  + 6.D0*Csum112(ep)
     &  + 12.D0*Csum122(ep)
     &  + 4.D0*Csum2222(ep)
     &  + 4.D0*Csum1112(ep)
     &  + 12.D0*Csum1122(ep)
     &  + 12.D0*Csum1222(ep)
     &  + Csum22222(ep)
     &  + Csum11112(ep)
     &  + 4.D0*Csum12222(ep)
     &  + 4.D0*Csum11122(ep)
     &  + 6.D0*Csum11222(ep)
     &

      Shat7(1,z6(1,1,1,1,2,2),ep)=
     &  - Csum11(ep)
     &  - 3.D0*Csum111(ep)
     &  - 3.D0*Csum112(ep)
     &  - 3.D0*Csum1111(ep)
     &  - 6.D0*Csum1112(ep)
     &  - 3.D0*Csum1122(ep)
     &  - Csum11111(ep)
     &  - 3.D0*Csum11112(ep)
     &  - 3.D0*Csum11122(ep)
     &  - Csum11222(ep)
     &

      Shat7(2,z6(1,1,1,1,2,2),ep)=
     &  - Csum11(ep)
     &  - 3.D0*Csum111(ep)
     &  - 3.D0*Csum112(ep)
     &  - 3.D0*Csum1111(ep)
     &  - 6.D0*Csum1112(ep)
     &  - 3.D0*Csum1122(ep)
     &  - Csum11111(ep)
     &  - 3.D0*Csum11112(ep)
     &  - 3.D0*Csum11122(ep)
     &  - Csum11222(ep)
     &

      Shat7(3,z6(1,1,1,1,2,2),ep)=
     &  - Csum11(ep)
     &  - 3.D0*Csum111(ep)
     &  - 3.D0*Csum112(ep)
     &  - 3.D0*Csum1111(ep)
     &  - 6.D0*Csum1112(ep)
     &  - 3.D0*Csum1122(ep)
     &  - Csum11111(ep)
     &  - 3.D0*Csum11112(ep)
     &  - 3.D0*Csum11122(ep)
     &  - Csum11222(ep)
     &  + Cv(cc111122 + C123,ep)
     &

      Shat7(1,z6(1,1,1,1,2,3),ep)=
     &  - Csum12(ep)
     &  - 3.D0*Csum112(ep)
     &  - 3.D0*Csum122(ep)
     &  - 3.D0*Csum1112(ep)
     &  - 6.D0*Csum1122(ep)
     &  - 3.D0*Csum1222(ep)
     &  - Csum11112(ep)
     &  - Csum12222(ep)
     &  - 3.D0*Csum11122(ep)
     &  - 3.D0*Csum11222(ep)
     &

      Shat7(2,z6(1,1,1,1,2,3),ep)=
     &  - Csum12(ep)
     &  - 3.D0*Csum112(ep)
     &  - 3.D0*Csum122(ep)
     &  - 3.D0*Csum1112(ep)
     &  - 6.D0*Csum1122(ep)
     &  - 3.D0*Csum1222(ep)
     &  - Csum11112(ep)
     &  - Csum12222(ep)
     &  - 3.D0*Csum11122(ep)
     &  - 3.D0*Csum11222(ep)
     &

      Shat7(3,z6(1,1,1,1,2,3),ep)=
     &  - Csum12(ep)
     &  - 3.D0*Csum112(ep)
     &  - 3.D0*Csum122(ep)
     &  - 3.D0*Csum1112(ep)
     &  - 6.D0*Csum1122(ep)
     &  - 3.D0*Csum1222(ep)
     &  - Csum11112(ep)
     &  - Csum12222(ep)
     &  - 3.D0*Csum11122(ep)
     &  - 3.D0*Csum11222(ep)
     &

      Shat7(1,z6(1,1,1,1,3,3),ep)=
     &  - Csum22(ep)
     &  - 3.D0*Csum222(ep)
     &  - 3.D0*Csum122(ep)
     &  - 3.D0*Csum2222(ep)
     &  - 3.D0*Csum1122(ep)
     &  - 6.D0*Csum1222(ep)
     &  - Csum22222(ep)
     &  - 3.D0*Csum12222(ep)
     &  - Csum11122(ep)
     &  - 3.D0*Csum11222(ep)
     &

      Shat7(2,z6(1,1,1,1,3,3),ep)=
     &  - Csum22(ep)
     &  - 3.D0*Csum222(ep)
     &  - 3.D0*Csum122(ep)
     &  - 3.D0*Csum2222(ep)
     &  - 3.D0*Csum1122(ep)
     &  - 6.D0*Csum1222(ep)
     &  - Csum22222(ep)
     &  - 3.D0*Csum12222(ep)
     &  - Csum11122(ep)
     &  - 3.D0*Csum11222(ep)
     &  + Cv(cc111122 + C124,ep)
     &

      Shat7(3,z6(1,1,1,1,3,3),ep)=
     &  - Csum22(ep)
     &  - 3.D0*Csum222(ep)
     &  - 3.D0*Csum122(ep)
     &  - 3.D0*Csum2222(ep)
     &  - 3.D0*Csum1122(ep)
     &  - 6.D0*Csum1222(ep)
     &  - Csum22222(ep)
     &  - 3.D0*Csum12222(ep)
     &  - Csum11122(ep)
     &  - 3.D0*Csum11222(ep)
     &

      Shat7(1,z6(1,1,1,2,2,2),ep)=
     &  + Csum111(ep)
     &  + 2.D0*Csum1111(ep)
     &  + 2.D0*Csum1112(ep)
     &  + Csum11111(ep)
     &  + 2.D0*Csum11112(ep)
     &  + Csum11122(ep)
     &

      Shat7(2,z6(1,1,1,2,2,2),ep)=
     &  + Csum111(ep)
     &  + 2.D0*Csum1111(ep)
     &  + 2.D0*Csum1112(ep)
     &  + Csum11111(ep)
     &  + 2.D0*Csum11112(ep)
     &  + Csum11122(ep)
     &

      Shat7(3,z6(1,1,1,2,2,2),ep)=
     &  + Csum111(ep)
     &  + 2.D0*Csum1111(ep)
     &  + 2.D0*Csum1112(ep)
     &  + Csum11111(ep)
     &  + 2.D0*Csum11112(ep)
     &  + Csum11122(ep)
     &  + Cv(cc111222 + C123,ep)
     &

      Shat7(1,z6(1,1,1,2,2,3),ep)=
     &  + Csum112(ep)
     &  + 2.D0*Csum1112(ep)
     &  + 2.D0*Csum1122(ep)
     &  + Csum11112(ep)
     &  + 2.D0*Csum11122(ep)
     &  + Csum11222(ep)
     &

      Shat7(2,z6(1,1,1,2,2,3),ep)=
     &  + Csum112(ep)
     &  + 2.D0*Csum1112(ep)
     &  + 2.D0*Csum1122(ep)
     &  + Csum11112(ep)
     &  + 2.D0*Csum11122(ep)
     &  + Csum11222(ep)
     &

      Shat7(3,z6(1,1,1,2,2,3),ep)=
     &  + Csum112(ep)
     &  + 2.D0*Csum1112(ep)
     &  + 2.D0*Csum1122(ep)
     &  + Csum11112(ep)
     &  + 2.D0*Csum11122(ep)
     &  + Csum11222(ep)
     &

      Shat7(1,z6(1,1,1,2,3,3),ep)=
     &  + Csum122(ep)
     &  + 2.D0*Csum1122(ep)
     &  + 2.D0*Csum1222(ep)
     &  + Csum12222(ep)
     &  + Csum11122(ep)
     &  + 2.D0*Csum11222(ep)
     &

      Shat7(2,z6(1,1,1,2,3,3),ep)=
     &  + Csum122(ep)
     &  + 2.D0*Csum1122(ep)
     &  + 2.D0*Csum1222(ep)
     &  + Csum12222(ep)
     &  + Csum11122(ep)
     &  + 2.D0*Csum11222(ep)
     &

      Shat7(3,z6(1,1,1,2,3,3),ep)=
     &  + Csum122(ep)
     &  + 2.D0*Csum1122(ep)
     &  + 2.D0*Csum1222(ep)
     &  + Csum12222(ep)
     &  + Csum11122(ep)
     &  + 2.D0*Csum11222(ep)
     &

      Shat7(1,z6(1,1,1,3,3,3),ep)=
     &  + Csum222(ep)
     &  + 2.D0*Csum2222(ep)
     &  + 2.D0*Csum1222(ep)
     &  + Csum22222(ep)
     &  + 2.D0*Csum12222(ep)
     &  + Csum11222(ep)
     &

      Shat7(2,z6(1,1,1,3,3,3),ep)=
     &  + Csum222(ep)
     &  + 2.D0*Csum2222(ep)
     &  + 2.D0*Csum1222(ep)
     &  + Csum22222(ep)
     &  + 2.D0*Csum12222(ep)
     &  + Csum11222(ep)
     &  + Cv(cc111222 + C124,ep)
     &

      Shat7(3,z6(1,1,1,3,3,3),ep)=
     &  + Csum222(ep)
     &  + 2.D0*Csum2222(ep)
     &  + 2.D0*Csum1222(ep)
     &  + Csum22222(ep)
     &  + 2.D0*Csum12222(ep)
     &  + Csum11222(ep)
     &

      Shat7(1,z6(1,1,2,2,2,2),ep)=
     &  - Csum1111(ep)
     &  - Csum11111(ep)
     &  - Csum11112(ep)
     &

      Shat7(2,z6(1,1,2,2,2,2),ep)=
     &  - Csum1111(ep)
     &  - Csum11111(ep)
     &  - Csum11112(ep)
     &

      Shat7(3,z6(1,1,2,2,2,2),ep)=
     &  - Csum1111(ep)
     &  - Csum11111(ep)
     &  - Csum11112(ep)
     &  + Cv(cc112222 + C123,ep)
     &

      Shat7(1,z6(1,1,2,2,2,3),ep)=
     &  - Csum1112(ep)
     &  - Csum11112(ep)
     &  - Csum11122(ep)
     &

      Shat7(2,z6(1,1,2,2,2,3),ep)=
     &  - Csum1112(ep)
     &  - Csum11112(ep)
     &  - Csum11122(ep)
     &

      Shat7(3,z6(1,1,2,2,2,3),ep)=
     &  - Csum1112(ep)
     &  - Csum11112(ep)
     &  - Csum11122(ep)
     &

      Shat7(1,z6(1,1,2,2,3,3),ep)=
     &  - Csum1122(ep)
     &  - Csum11122(ep)
     &  - Csum11222(ep)
     &

      Shat7(2,z6(1,1,2,2,3,3),ep)=
     &  - Csum1122(ep)
     &  - Csum11122(ep)
     &  - Csum11222(ep)
     &

      Shat7(3,z6(1,1,2,2,3,3),ep)=
     &  - Csum1122(ep)
     &  - Csum11122(ep)
     &  - Csum11222(ep)
     &

      Shat7(1,z6(1,1,2,3,3,3),ep)=
     &  - Csum1222(ep)
     &  - Csum12222(ep)
     &  - Csum11222(ep)
     &

      Shat7(2,z6(1,1,2,3,3,3),ep)=
     &  - Csum1222(ep)
     &  - Csum12222(ep)
     &  - Csum11222(ep)
     &

      Shat7(3,z6(1,1,2,3,3,3),ep)=
     &  - Csum1222(ep)
     &  - Csum12222(ep)
     &  - Csum11222(ep)
     &

      Shat7(1,z6(1,1,3,3,3,3),ep)=
     &  - Csum2222(ep)
     &  - Csum22222(ep)
     &  - Csum12222(ep)
     &

      Shat7(2,z6(1,1,3,3,3,3),ep)=
     &  - Csum2222(ep)
     &  - Csum22222(ep)
     &  - Csum12222(ep)
     &  + Cv(cc112222 + C124,ep)
     &

      Shat7(3,z6(1,1,3,3,3,3),ep)=
     &  - Csum2222(ep)
     &  - Csum22222(ep)
     &  - Csum12222(ep)
     &

      Shat7(1,z6(1,2,2,2,2,2),ep)=
     &  + Csum11111(ep)
     &

      Shat7(2,z6(1,2,2,2,2,2),ep)=
     &  + Csum11111(ep)
     &

      Shat7(3,z6(1,2,2,2,2,2),ep)=
     &  + Csum11111(ep)
     &  + Cv(cc122222 + C123,ep)
     &

      Shat7(1,z6(1,2,2,2,2,3),ep)=
     &  + Csum11112(ep)
     &

      Shat7(2,z6(1,2,2,2,2,3),ep)=
     &  + Csum11112(ep)
     &

      Shat7(3,z6(1,2,2,2,2,3),ep)=
     &  + Csum11112(ep)
     &

      Shat7(1,z6(1,2,2,2,3,3),ep)=
     &  + Csum11122(ep)
     &

      Shat7(2,z6(1,2,2,2,3,3),ep)=
     &  + Csum11122(ep)
     &

      Shat7(3,z6(1,2,2,2,3,3),ep)=
     &  + Csum11122(ep)
     &

      Shat7(1,z6(1,2,2,3,3,3),ep)=
     &  + Csum11222(ep)
     &

      Shat7(2,z6(1,2,2,3,3,3),ep)=
     &  + Csum11222(ep)
     &

      Shat7(3,z6(1,2,2,3,3,3),ep)=
     &  + Csum11222(ep)
     &

      Shat7(1,z6(1,2,3,3,3,3),ep)=
     &  + Csum12222(ep)
     &

      Shat7(2,z6(1,2,3,3,3,3),ep)=
     &  + Csum12222(ep)
     &

      Shat7(3,z6(1,2,3,3,3,3),ep)=
     &  + Csum12222(ep)
     &

      Shat7(1,z6(1,3,3,3,3,3),ep)=
     &  + Csum22222(ep)
     &

      Shat7(2,z6(1,3,3,3,3,3),ep)=
     &  + Csum22222(ep)
     &  + Cv(cc122222 + C124,ep)
     &

      Shat7(3,z6(1,3,3,3,3,3),ep)=
     &  + Csum22222(ep)
     &

      Shat7(1,z6(2,2,2,2,2,2),ep)=
     &  + Cv(cc111111 + C134,ep)
     &  - Cv(cc111111 + C234,ep)
     &

      Shat7(2,z6(2,2,2,2,2,2),ep)=
     &  - Cv(cc111111 + C234,ep)
     &

      Shat7(3,z6(2,2,2,2,2,2),ep)=
     &  - Cv(cc111111 + C234,ep)
     &  + Cv(cc222222 + C123,ep)
     &

      Shat7(1,z6(2,2,2,2,2,3),ep)=
     &  + Cv(cc111112 + C134,ep)
     &  - Cv(cc111112 + C234,ep)
     &

      Shat7(2,z6(2,2,2,2,2,3),ep)=
     &  - Cv(cc111112 + C234,ep)
     &

      Shat7(3,z6(2,2,2,2,2,3),ep)=
     &  - Cv(cc111112 + C234,ep)
     &

      Shat7(1,z6(2,2,2,2,3,3),ep)=
     &  + Cv(cc111122 + C134,ep)
     &  - Cv(cc111122 + C234,ep)
     &

      Shat7(2,z6(2,2,2,2,3,3),ep)=
     &  - Cv(cc111122 + C234,ep)
     &

      Shat7(3,z6(2,2,2,2,3,3),ep)=
     &  - Cv(cc111122 + C234,ep)
     &

      Shat7(1,z6(2,2,2,3,3,3),ep)=
     &  + Cv(cc111222 + C134,ep)
     &  - Cv(cc111222 + C234,ep)
     &

      Shat7(2,z6(2,2,2,3,3,3),ep)=
     &  - Cv(cc111222 + C234,ep)
     &

      Shat7(3,z6(2,2,2,3,3,3),ep)=
     &  - Cv(cc111222 + C234,ep)
     &

      Shat7(1,z6(2,2,3,3,3,3),ep)=
     &  + Cv(cc112222 + C134,ep)
     &  - Cv(cc112222 + C234,ep)
     &

      Shat7(2,z6(2,2,3,3,3,3),ep)=
     &  - Cv(cc112222 + C234,ep)
     &

      Shat7(3,z6(2,2,3,3,3,3),ep)=
     &  - Cv(cc112222 + C234,ep)
     &

      Shat7(1,z6(2,3,3,3,3,3),ep)=
     &  + Cv(cc122222 + C134,ep)
     &  - Cv(cc122222 + C234,ep)
     &

      Shat7(2,z6(2,3,3,3,3,3),ep)=
     &  - Cv(cc122222 + C234,ep)
     &

      Shat7(3,z6(2,3,3,3,3,3),ep)=
     &  - Cv(cc122222 + C234,ep)
     &

      Shat7(1,z6(3,3,3,3,3,3),ep)=
     &  + Cv(cc222222 + C134,ep)
     &  - Cv(cc222222 + C234,ep)
     &

      Shat7(2,z6(3,3,3,3,3,3),ep)=
     &  + Cv(cc222222 + C124,ep)
     &  - Cv(cc222222 + C234,ep)
     &

      Shat7(3,z6(3,3,3,3,3,3),ep)=
     &  - Cv(cc222222 + C234,ep)


      Shat6(1,z5(1,1,1,1,1),ep)=
     &  + Csum0(ep)
     &  + 4.D0*Csum1(ep)
     &  + 4.D0*Csum2(ep)
     &  + 4.D0*Csum222(ep)
     &  + 4.D0*Csum111(ep)
     &  + 6.D0*Csum11(ep)
     &  + 6.D0*Csum22(ep)
     &  + 12.D0*Csum12(ep)
     &  + 12.D0*Csum112(ep)
     &  + 12.D0*Csum122(ep)
     &  + 6.D0*Csum1122(ep)
     &  + 4.D0*Csum1112(ep)
     &  + 4.D0*Csum1222(ep)
     &  + Csum1111(ep)
     &  + Csum2222(ep)
     &

      Shat6(2,z5(1,1,1,1,1),ep)=
     &  + Csum0(ep)
     &  + 4.D0*Csum1(ep)
     &  + 4.D0*Csum2(ep)
     &  + 4.D0*Csum222(ep)
     &  + 4.D0*Csum111(ep)
     &  + 6.D0*Csum11(ep)
     &  + 6.D0*Csum22(ep)
     &  + 12.D0*Csum12(ep)
     &  + 12.D0*Csum112(ep)
     &  + 12.D0*Csum122(ep)
     &  + 6.D0*Csum1122(ep)
     &  + 4.D0*Csum1112(ep)
     &  + 4.D0*Csum1222(ep)
     &  + Csum1111(ep)
     &  + Csum2222(ep)
     &  + Cv(cc11111 + C124,ep)
     &

      Shat6(3,z5(1,1,1,1,1),ep)=
     &  + Csum0(ep)
     &  + 4.D0*Csum1(ep)
     &  + 4.D0*Csum2(ep)
     &  + 4.D0*Csum222(ep)
     &  + 4.D0*Csum111(ep)
     &  + 6.D0*Csum11(ep)
     &  + 6.D0*Csum22(ep)
     &  + 12.D0*Csum12(ep)
     &  + 12.D0*Csum112(ep)
     &  + 12.D0*Csum122(ep)
     &  + 6.D0*Csum1122(ep)
     &  + 4.D0*Csum1112(ep)
     &  + 4.D0*Csum1222(ep)
     &  + Csum1111(ep)
     &  + Csum2222(ep)
     &  + Cv(cc11111 + C123,ep)
     &

      Shat6(1,z5(1,1,1,1,2),ep)=
     &  - Csum1(ep)
     &  - 3.D0*Csum111(ep)
     &  - 3.D0*Csum11(ep)
     &  - 3.D0*Csum12(ep)
     &  - 6.D0*Csum112(ep)
     &  - 3.D0*Csum122(ep)
     &  - 3.D0*Csum1122(ep)
     &  - 3.D0*Csum1112(ep)
     &  - Csum1222(ep)
     &  - Csum1111(ep)
     &

      Shat6(2,z5(1,1,1,1,2),ep)=
     &  - Csum1(ep)
     &  - 3.D0*Csum111(ep)
     &  - 3.D0*Csum11(ep)
     &  - 3.D0*Csum12(ep)
     &  - 6.D0*Csum112(ep)
     &  - 3.D0*Csum122(ep)
     &  - 3.D0*Csum1122(ep)
     &  - 3.D0*Csum1112(ep)
     &  - Csum1222(ep)
     &  - Csum1111(ep)
     &

      Shat6(3,z5(1,1,1,1,2),ep)=
     &  - Csum1(ep)
     &  - 3.D0*Csum111(ep)
     &  - 3.D0*Csum11(ep)
     &  - 3.D0*Csum12(ep)
     &  - 6.D0*Csum112(ep)
     &  - 3.D0*Csum122(ep)
     &  - 3.D0*Csum1122(ep)
     &  - 3.D0*Csum1112(ep)
     &  - Csum1222(ep)
     &  - Csum1111(ep)
     &  + Cv(cc11112 + C123,ep)
     &

      Shat6(1,z5(1,1,1,1,3),ep)=
     &  - Csum2(ep)
     &  - 3.D0*Csum222(ep)
     &  - 3.D0*Csum22(ep)
     &  - 3.D0*Csum12(ep)
     &  - 3.D0*Csum112(ep)
     &  - 6.D0*Csum122(ep)
     &  - 3.D0*Csum1122(ep)
     &  - Csum1112(ep)
     &  - 3.D0*Csum1222(ep)
     &  - Csum2222(ep)
     &

      Shat6(2,z5(1,1,1,1,3),ep)=
     &  - Csum2(ep)
     &  - 3.D0*Csum222(ep)
     &  - 3.D0*Csum22(ep)
     &  - 3.D0*Csum12(ep)
     &  - 3.D0*Csum112(ep)
     &  - 6.D0*Csum122(ep)
     &  - 3.D0*Csum1122(ep)
     &  - Csum1112(ep)
     &  - 3.D0*Csum1222(ep)
     &  - Csum2222(ep)
     &  + Cv(cc11112 + C124,ep)
     &

      Shat6(3,z5(1,1,1,1,3),ep)=
     &  - Csum2(ep)
     &  - 3.D0*Csum222(ep)
     &  - 3.D0*Csum22(ep)
     &  - 3.D0*Csum12(ep)
     &  - 3.D0*Csum112(ep)
     &  - 6.D0*Csum122(ep)
     &  - 3.D0*Csum1122(ep)
     &  - Csum1112(ep)
     &  - 3.D0*Csum1222(ep)
     &  - Csum2222(ep)
     &

      Shat6(1,z5(1,1,1,2,2),ep)=
     &  + 2.D0*Csum111(ep)
     &  + Csum11(ep)
     &  + 2.D0*Csum112(ep)
     &  + Csum1122(ep)
     &  + 2.D0*Csum1112(ep)
     &  + Csum1111(ep)
     &

      Shat6(2,z5(1,1,1,2,2),ep)=
     &  + 2.D0*Csum111(ep)
     &  + Csum11(ep)
     &  + 2.D0*Csum112(ep)
     &  + Csum1122(ep)
     &  + 2.D0*Csum1112(ep)
     &  + Csum1111(ep)
     &

      Shat6(3,z5(1,1,1,2,2),ep)=
     &  + 2.D0*Csum111(ep)
     &  + Csum11(ep)
     &  + 2.D0*Csum112(ep)
     &  + Csum1122(ep)
     &  + 2.D0*Csum1112(ep)
     &  + Csum1111(ep)
     &  + Cv(cc11122 + C123,ep)
     &

      Shat6(1,z5(1,1,1,2,3),ep)=
     &  + Csum12(ep)
     &  + 2.D0*Csum112(ep)
     &  + 2.D0*Csum122(ep)
     &  + 2.D0*Csum1122(ep)
     &  + Csum1112(ep)
     &  + Csum1222(ep)
     &

      Shat6(2,z5(1,1,1,2,3),ep)=
     &  + Csum12(ep)
     &  + 2.D0*Csum112(ep)
     &  + 2.D0*Csum122(ep)
     &  + 2.D0*Csum1122(ep)
     &  + Csum1112(ep)
     &  + Csum1222(ep)
     &

      Shat6(3,z5(1,1,1,2,3),ep)=
     &  + Csum12(ep)
     &  + 2.D0*Csum112(ep)
     &  + 2.D0*Csum122(ep)
     &  + 2.D0*Csum1122(ep)
     &  + Csum1112(ep)
     &  + Csum1222(ep)
     &

      Shat6(1,z5(1,1,1,3,3),ep)=
     &  + 2.D0*Csum222(ep)
     &  + Csum22(ep)
     &  + 2.D0*Csum122(ep)
     &  + Csum1122(ep)
     &  + 2.D0*Csum1222(ep)
     &  + Csum2222(ep)
     &

      Shat6(2,z5(1,1,1,3,3),ep)=
     &  + 2.D0*Csum222(ep)
     &  + Csum22(ep)
     &  + 2.D0*Csum122(ep)
     &  + Csum1122(ep)
     &  + 2.D0*Csum1222(ep)
     &  + Csum2222(ep)
     &  + Cv(cc11122 + C124,ep)
     &

      Shat6(3,z5(1,1,1,3,3),ep)=
     &  + 2.D0*Csum222(ep)
     &  + Csum22(ep)
     &  + 2.D0*Csum122(ep)
     &  + Csum1122(ep)
     &  + 2.D0*Csum1222(ep)
     &  + Csum2222(ep)
     &

      Shat6(1,z5(1,1,2,2,2),ep)=
     &  - Csum111(ep)
     &  - Csum1112(ep)
     &  - Csum1111(ep)
     &

      Shat6(2,z5(1,1,2,2,2),ep)=
     &  - Csum111(ep)
     &  - Csum1112(ep)
     &  - Csum1111(ep)
     &

      Shat6(3,z5(1,1,2,2,2),ep)=
     &  - Csum111(ep)
     &  - Csum1112(ep)
     &  - Csum1111(ep)
     &  + Cv(cc11222 + C123,ep)
     &

      Shat6(1,z5(1,1,2,2,3),ep)=
     &  - Csum112(ep)
     &  - Csum1122(ep)
     &  - Csum1112(ep)
     &

      Shat6(2,z5(1,1,2,2,3),ep)=
     &  - Csum112(ep)
     &  - Csum1122(ep)
     &  - Csum1112(ep)
     &

      Shat6(3,z5(1,1,2,2,3),ep)=
     &  - Csum112(ep)
     &  - Csum1122(ep)
     &  - Csum1112(ep)
     &

      Shat6(1,z5(1,1,2,3,3),ep)=
     &  - Csum122(ep)
     &  - Csum1122(ep)
     &  - Csum1222(ep)
     &

      Shat6(2,z5(1,1,2,3,3),ep)=
     &  - Csum122(ep)
     &  - Csum1122(ep)
     &  - Csum1222(ep)
     &

      Shat6(3,z5(1,1,2,3,3),ep)=
     &  - Csum122(ep)
     &  - Csum1122(ep)
     &  - Csum1222(ep)
     &

      Shat6(1,z5(1,1,3,3,3),ep)=
     &  - Csum222(ep)
     &  - Csum1222(ep)
     &  - Csum2222(ep)
     &

      Shat6(2,z5(1,1,3,3,3),ep)=
     &  - Csum222(ep)
     &  - Csum1222(ep)
     &  - Csum2222(ep)
     &  + Cv(cc11222 + C124,ep)
     &

      Shat6(3,z5(1,1,3,3,3),ep)=
     &  - Csum222(ep)
     &  - Csum1222(ep)
     &  - Csum2222(ep)
     &

      Shat6(1,z5(1,2,2,2,2),ep)=
     &  + Csum1111(ep)
     &

      Shat6(2,z5(1,2,2,2,2),ep)=
     &  + Csum1111(ep)
     &

      Shat6(3,z5(1,2,2,2,2),ep)=
     &  + Csum1111(ep)
     &  + Cv(cc12222 + C123,ep)
     &

      Shat6(1,z5(1,2,2,2,3),ep)=
     &  + Csum1112(ep)
     &

      Shat6(2,z5(1,2,2,2,3),ep)=
     &  + Csum1112(ep)
     &

      Shat6(3,z5(1,2,2,2,3),ep)=
     &  + Csum1112(ep)
     &

      Shat6(1,z5(1,2,2,3,3),ep)=
     &  + Csum1122(ep)
     &

      Shat6(2,z5(1,2,2,3,3),ep)=
     &  + Csum1122(ep)
     &

      Shat6(3,z5(1,2,2,3,3),ep)=
     &  + Csum1122(ep)
     &

      Shat6(1,z5(1,2,3,3,3),ep)=
     &  + Csum1222(ep)
     &

      Shat6(2,z5(1,2,3,3,3),ep)=
     &  + Csum1222(ep)
     &

      Shat6(3,z5(1,2,3,3,3),ep)=
     &  + Csum1222(ep)
     &

      Shat6(1,z5(1,3,3,3,3),ep)=
     &  + Csum2222(ep)
     &

      Shat6(2,z5(1,3,3,3,3),ep)=
     &  + Csum2222(ep)
     &  + Cv(cc12222 + C124,ep)
     &

      Shat6(3,z5(1,3,3,3,3),ep)=
     &  + Csum2222(ep)
     &

      Shat6(1,z5(2,2,2,2,2),ep)=
     &  + Cv(cc11111 + C134,ep)
     &  - Cv(cc11111 + C234,ep)
     &

      Shat6(2,z5(2,2,2,2,2),ep)=
     &  - Cv(cc11111 + C234,ep)
     &

      Shat6(3,z5(2,2,2,2,2),ep)=
     &  - Cv(cc11111 + C234,ep)
     &  + Cv(cc22222 + C123,ep)
     &

      Shat6(1,z5(2,2,2,2,3),ep)=
     &  + Cv(cc11112 + C134,ep)
     &  - Cv(cc11112 + C234,ep)
     &

      Shat6(2,z5(2,2,2,2,3),ep)=
     &  - Cv(cc11112 + C234,ep)
     &

      Shat6(3,z5(2,2,2,2,3),ep)=
     &  - Cv(cc11112 + C234,ep)
     &

      Shat6(1,z5(2,2,2,3,3),ep)=
     &  + Cv(cc11122 + C134,ep)
     &  - Cv(cc11122 + C234,ep)
     &

      Shat6(2,z5(2,2,2,3,3),ep)=
     &  - Cv(cc11122 + C234,ep)
     &

      Shat6(3,z5(2,2,2,3,3),ep)=
     &  - Cv(cc11122 + C234,ep)
     &

      Shat6(1,z5(2,2,3,3,3),ep)=
     &  + Cv(cc11222 + C134,ep)
     &  - Cv(cc11222 + C234,ep)
     &

      Shat6(2,z5(2,2,3,3,3),ep)=
     &  - Cv(cc11222 + C234,ep)
     &

      Shat6(3,z5(2,2,3,3,3),ep)=
     &  - Cv(cc11222 + C234,ep)
     &

      Shat6(1,z5(2,3,3,3,3),ep)=
     &  + Cv(cc12222 + C134,ep)
     &  - Cv(cc12222 + C234,ep)
     &

      Shat6(2,z5(2,3,3,3,3),ep)=
     &  - Cv(cc12222 + C234,ep)
     &

      Shat6(3,z5(2,3,3,3,3),ep)=
     &  - Cv(cc12222 + C234,ep)
     &

      Shat6(1,z5(3,3,3,3,3),ep)=
     &  + Cv(cc22222 + C134,ep)
     &  - Cv(cc22222 + C234,ep)
     &

      Shat6(2,z5(3,3,3,3,3),ep)=
     &  + Cv(cc22222 + C124,ep)
     &  - Cv(cc22222 + C234,ep)
     &

      Shat6(3,z5(3,3,3,3,3),ep)=
     &  - Cv(cc22222 + C234,ep)
     &

      Shat5(1,z4(1,1,1,1),ep)=
     &  - Csum0(ep)
     &  - 3.D0*Csum1(ep)
     &  - 3.D0*Csum2(ep)
     &  - Csum222(ep)
     &  - Csum111(ep)
     &  - 3.D0*Csum11(ep)
     &  - 3.D0*Csum22(ep)
     &  - 6.D0*Csum12(ep)
     &  - 3.D0*Csum112(ep)
     &  - 3.D0*Csum122(ep)
     &

      Shat5(2,z4(1,1,1,1),ep)=
     &  - Csum0(ep)
     &  - 3.D0*Csum1(ep)
     &  - 3.D0*Csum2(ep)
     &  - Csum222(ep)
     &  - Csum111(ep)
     &  - 3.D0*Csum11(ep)
     &  - 3.D0*Csum22(ep)
     &  - 6.D0*Csum12(ep)
     &  - 3.D0*Csum112(ep)
     &  - 3.D0*Csum122(ep)
     &  + Cv(cc1111 + C124,ep)
     &

      Shat5(3,z4(1,1,1,1),ep)=
     &  - Csum0(ep)
     &  - 3.D0*Csum1(ep)
     &  - 3.D0*Csum2(ep)
     &  - Csum222(ep)
     &  - Csum111(ep)
     &  - 3.D0*Csum11(ep)
     &  - 3.D0*Csum22(ep)
     &  - 6.D0*Csum12(ep)
     &  - 3.D0*Csum112(ep)
     &  - 3.D0*Csum122(ep)
     &  + Cv(cc1111 + C123,ep)
     &

      Shat5(1,z4(1,1,1,2),ep)=
     &  + Csum1(ep)
     &  + Csum111(ep)
     &  + 2.D0*Csum11(ep)
     &  + 2.D0*Csum12(ep)
     &  + 2.D0*Csum112(ep)
     &  + Csum122(ep)
     &

      Shat5(2,z4(1,1,1,2),ep)=
     &  + Csum1(ep)
     &  + Csum111(ep)
     &  + 2.D0*Csum11(ep)
     &  + 2.D0*Csum12(ep)
     &  + 2.D0*Csum112(ep)
     &  + Csum122(ep)
     &

      Shat5(3,z4(1,1,1,2),ep)=
     &  + Csum1(ep)
     &  + Csum111(ep)
     &  + 2.D0*Csum11(ep)
     &  + 2.D0*Csum12(ep)
     &  + 2.D0*Csum112(ep)
     &  + Csum122(ep)
     &  + Cv(cc1112 + C123,ep)
     &

      Shat5(1,z4(1,1,1,3),ep)=
     &  + Csum2(ep)
     &  + Csum222(ep)
     &  + 2.D0*Csum22(ep)
     &  + 2.D0*Csum12(ep)
     &  + Csum112(ep)
     &  + 2.D0*Csum122(ep)
     &

      Shat5(2,z4(1,1,1,3),ep)=
     &  + Csum2(ep)
     &  + Csum222(ep)
     &  + 2.D0*Csum22(ep)
     &  + 2.D0*Csum12(ep)
     &  + Csum112(ep)
     &  + 2.D0*Csum122(ep)
     &  + Cv(cc1112 + C124,ep)
     &

      Shat5(3,z4(1,1,1,3),ep)=
     &  + Csum2(ep)
     &  + Csum222(ep)
     &  + 2.D0*Csum22(ep)
     &  + 2.D0*Csum12(ep)
     &  + Csum112(ep)
     &  + 2.D0*Csum122(ep)
     &

      Shat5(1,z4(1,1,2,2),ep)=
     &  - Csum111(ep)
     &  - Csum11(ep)
     &  - Csum112(ep)
     &

      Shat5(2,z4(1,1,2,2),ep)=
     &  - Csum111(ep)
     &  - Csum11(ep)
     &  - Csum112(ep)
     &

      Shat5(3,z4(1,1,2,2),ep)=
     &  - Csum111(ep)
     &  - Csum11(ep)
     &  - Csum112(ep)
     &  + Cv(cc1122 + C123,ep)
     &

      Shat5(1,z4(1,1,2,3),ep)=
     &  - Csum12(ep)
     &  - Csum112(ep)
     &  - Csum122(ep)
     &

      Shat5(2,z4(1,1,2,3),ep)=
     &  - Csum12(ep)
     &  - Csum112(ep)
     &  - Csum122(ep)
     &

      Shat5(3,z4(1,1,2,3),ep)=
     &  - Csum12(ep)
     &  - Csum112(ep)
     &  - Csum122(ep)
     &

      Shat5(1,z4(1,1,3,3),ep)=
     &  - Csum222(ep)
     &  - Csum22(ep)
     &  - Csum122(ep)
     &

      Shat5(2,z4(1,1,3,3),ep)=
     &  - Csum222(ep)
     &  - Csum22(ep)
     &  - Csum122(ep)
     &  + Cv(cc1122 + C124,ep)
     &

      Shat5(3,z4(1,1,3,3),ep)=
     &  - Csum222(ep)
     &  - Csum22(ep)
     &  - Csum122(ep)
     &

      Shat5(1,z4(1,2,2,2),ep)=
     &  + Csum111(ep)
     &

      Shat5(2,z4(1,2,2,2),ep)=
     &  + Csum111(ep)
     &

      Shat5(3,z4(1,2,2,2),ep)=
     &  + Csum111(ep)
     &  + Cv(cc1222 + C123,ep)
     &

      Shat5(1,z4(1,2,2,3),ep)=
     &  + Csum112(ep)
     &

      Shat5(2,z4(1,2,2,3),ep)=
     &  + Csum112(ep)
     &

      Shat5(3,z4(1,2,2,3),ep)=
     &  + Csum112(ep)
     &

      Shat5(1,z4(1,2,3,3),ep)=
     &  + Csum122(ep)
     &

      Shat5(2,z4(1,2,3,3),ep)=
     &  + Csum122(ep)
     &

      Shat5(3,z4(1,2,3,3),ep)=
     &  + Csum122(ep)
     &

      Shat5(1,z4(1,3,3,3),ep)=
     &  + Csum222(ep)
     &

      Shat5(2,z4(1,3,3,3),ep)=
     &  + Csum222(ep)
     &  + Cv(cc1222 + C124,ep)
     &

      Shat5(3,z4(1,3,3,3),ep)=
     &  + Csum222(ep)
     &

      Shat5(1,z4(2,2,2,2),ep)=
     &  + Cv(cc1111 + C134,ep)
     &  - Cv(cc1111 + C234,ep)
     &

      Shat5(2,z4(2,2,2,2),ep)=
     &  - Cv(cc1111 + C234,ep)
     &

      Shat5(3,z4(2,2,2,2),ep)=
     &  - Cv(cc1111 + C234,ep)
     &  + Cv(cc2222 + C123,ep)
     &

      Shat5(1,z4(2,2,2,3),ep)=
     &  + Cv(cc1112 + C134,ep)
     &  - Cv(cc1112 + C234,ep)
     &

      Shat5(2,z4(2,2,2,3),ep)=
     &  - Cv(cc1112 + C234,ep)
     &

      Shat5(3,z4(2,2,2,3),ep)=
     &  - Cv(cc1112 + C234,ep)
     &

      Shat5(1,z4(2,2,3,3),ep)=
     &  + Cv(cc1122 + C134,ep)
     &  - Cv(cc1122 + C234,ep)
     &

      Shat5(2,z4(2,2,3,3),ep)=
     &  - Cv(cc1122 + C234,ep)
     &

      Shat5(3,z4(2,2,3,3),ep)=
     &  - Cv(cc1122 + C234,ep)
     &

      Shat5(1,z4(2,3,3,3),ep)=
     &  + Cv(cc1222 + C134,ep)
     &  - Cv(cc1222 + C234,ep)
     &

      Shat5(2,z4(2,3,3,3),ep)=
     &  - Cv(cc1222 + C234,ep)
     &

      Shat5(3,z4(2,3,3,3),ep)=
     &  - Cv(cc1222 + C234,ep)
     &

      Shat5(1,z4(3,3,3,3),ep)=
     &  + Cv(cc2222 + C134,ep)
     &  - Cv(cc2222 + C234,ep)
     &

      Shat5(2,z4(3,3,3,3),ep)=
     &  + Cv(cc2222 + C124,ep)
     &  - Cv(cc2222 + C234,ep)
     &

      Shat5(3,z4(3,3,3,3),ep)=
     &  - Cv(cc2222 + C234,ep)
     &

      Shat4(1,z3(1,1,1),ep)=
     &  + Csum0(ep)
     &  + 2.D0*Csum1(ep)
     &  + 2.D0*Csum2(ep)
     &  + Csum11(ep)
     &  + Csum22(ep)
     &  + 2.D0*Csum12(ep)
     &

      Shat4(2,z3(1,1,1),ep)=
     &  + Csum0(ep)
     &  + 2.D0*Csum1(ep)
     &  + 2.D0*Csum2(ep)
     &  + Csum11(ep)
     &  + Csum22(ep)
     &  + 2.D0*Csum12(ep)
     &  + Cv(cc111 + C124,ep)
     &

      Shat4(3,z3(1,1,1),ep)=
     &  + Csum0(ep)
     &  + 2.D0*Csum1(ep)
     &  + 2.D0*Csum2(ep)
     &  + Csum11(ep)
     &  + Csum22(ep)
     &  + 2.D0*Csum12(ep)
     &  + Cv(cc111 + C123,ep)
     &

      Shat4(1,z3(2,2,2),ep)=
     &  + Cv(cc111 + C134,ep)
     &  - Cv(cc111 + C234,ep)
     &

      Shat4(2,z3(2,2,2),ep)=
     &  - Cv(cc111 + C234,ep)
     &

      Shat4(3,z3(2,2,2),ep)=
     &  - Cv(cc111 + C234,ep)
     &  + Cv(cc222 + C123,ep)
     &

      Shat4(1,z3(3,3,3),ep)=
     &  + Cv(cc222 + C134,ep)
     &  - Cv(cc222 + C234,ep)
     &

      Shat4(2,z3(3,3,3),ep)=
     &  + Cv(cc222 + C124,ep)
     &  - Cv(cc222 + C234,ep)
     &

      Shat4(3,z3(3,3,3),ep)=
     &  - Cv(cc222 + C234,ep)
     &

      Shat4(1,z3(1,1,2),ep)=
     &  - Csum1(ep)
     &  - Csum11(ep)
     &  - Csum12(ep)
     &

      Shat4(2,z3(1,1,2),ep)=
     &  - Csum1(ep)
     &  - Csum11(ep)
     &  - Csum12(ep)
     &

      Shat4(3,z3(1,1,2),ep)=
     &  - Csum1(ep)
     &  - Csum11(ep)
     &  - Csum12(ep)
     &  + Cv(cc112 + C123,ep)
     &

      Shat4(1,z3(1,1,3),ep)=
     &  - Csum2(ep)
     &  - Csum22(ep)
     &  - Csum12(ep)
     &

      Shat4(2,z3(1,1,3),ep)=
     &  - Csum2(ep)
     &  - Csum22(ep)
     &  - Csum12(ep)
     &  + Cv(cc112 + C124,ep)
     &

      Shat4(3,z3(1,1,3),ep)=
     &  - Csum2(ep)
     &  - Csum22(ep)
     &  - Csum12(ep)
     &

      Shat4(1,z3(1,2,2),ep)=
     &  + Csum11(ep)
     &

      Shat4(2,z3(1,2,2),ep)=
     &  + Csum11(ep)
     &

      Shat4(3,z3(1,2,2),ep)=
     &  + Csum11(ep)
     &  + Cv(cc122 + C123,ep)
     &

      Shat4(1,z3(1,3,3),ep)=
     &  + Csum22(ep)
     &

      Shat4(2,z3(1,3,3),ep)=
     &  + Csum22(ep)
     &  + Cv(cc122 + C124,ep)
     &

      Shat4(3,z3(1,3,3),ep)=
     &  + Csum22(ep)
     &

      Shat4(1,z3(2,2,3),ep)=
     &  + Cv(cc112 + C134,ep)
     &  - Cv(cc112 + C234,ep)
     &

      Shat4(2,z3(2,2,3),ep)=
     &  - Cv(cc112 + C234,ep)
     &

      Shat4(3,z3(2,2,3),ep)=
     &  - Cv(cc112 + C234,ep)
     &

      Shat4(1,z3(2,3,3),ep)=
     &  + Cv(cc122 + C134,ep)
     &  - Cv(cc122 + C234,ep)
     &

      Shat4(2,z3(2,3,3),ep)=
     &  - Cv(cc122 + C234,ep)
     &

      Shat4(3,z3(2,3,3),ep)=
     &  - Cv(cc122 + C234,ep)
     &

      Shat4(1,z3(1,2,3),ep)=
     &  + Csum12(ep)
     &

      Shat4(2,z3(1,2,3),ep)=
     &  + Csum12(ep)
     &

      Shat4(3,z3(1,2,3),ep)=
     &  + Csum12(ep)
     &

      Shat3(1,z2(1,1),ep)=
     &  - Csum0(ep)
     &  - Csum1(ep)
     &  - Csum2(ep)
     &

      Shat3(2,z2(1,1),ep)=
     &  - Csum0(ep)
     &  - Csum1(ep)
     &  - Csum2(ep)
     &  + Cv(cc11 + C124,ep)
     &

      Shat3(3,z2(1,1),ep)=
     &  - Csum0(ep)
     &  - Csum1(ep)
     &  - Csum2(ep)
     &  + Cv(cc11 + C123,ep)
     &

      Shat3(1,z2(2,2),ep)=
     &  + Cv(cc11 + C134,ep)
     &  - Cv(cc11 + C234,ep)
     &

      Shat3(2,z2(2,2),ep)=
     &  - Cv(cc11 + C234,ep)
     &

      Shat3(3,z2(2,2),ep)=
     &  - Cv(cc11 + C234,ep)
     &  + Cv(cc22 + C123,ep)
     &

      Shat3(1,z2(3,3),ep)=
     &  + Cv(cc22 + C134,ep)
     &  - Cv(cc22 + C234,ep)
     &

      Shat3(2,z2(3,3),ep)=
     &  + Cv(cc22 + C124,ep)
     &  - Cv(cc22 + C234,ep)
     &

      Shat3(3,z2(3,3),ep)=
     &  - Cv(cc22 + C234,ep)
     &

      Shat3(1,z2(1,2),ep)=
     &  + Csum1(ep)
     &

      Shat3(2,z2(1,2),ep)=
     &  + Csum1(ep)
     &

      Shat3(3,z2(1,2),ep)=
     &  + Csum1(ep)
     &  + Cv(cc12 + C123,ep)
     &

      Shat3(1,z2(1,3),ep)=
     &  + Csum2(ep)
     &

      Shat3(2,z2(1,3),ep)=
     &  + Csum2(ep)
     &  + Cv(cc12 + C124,ep)
     &

      Shat3(3,z2(1,3),ep)=
     &  + Csum2(ep)
     &

      Shat3(1,z2(2,3),ep)=
     &  + Cv(cc12 + C134,ep)
     &  - Cv(cc12 + C234,ep)
     &

      Shat3(2,z2(2,3),ep)=
     &  - Cv(cc12 + C234,ep)
     &

      Shat3(3,z2(2,3),ep)=
     &  - Cv(cc12 + C234,ep)
     &

      Shat2(1,1,ep)=
     &  + Csum0(ep)
     &

      Shat2(2,1,ep)=
     &  + Csum0(ep)
     &  + Cv(cc1 + C124,ep)
     &

      Shat2(3,1,ep)=
     &  + Csum0(ep)
     &  + Cv(cc1 + C123,ep)
     &

      Shat2(1,2,ep)=
     &  + Cv(cc1 + C134,ep)
     &  - Cv(cc1 + C234,ep)
     &

      Shat2(2,2,ep)=
     &  - Cv(cc1 + C234,ep)
     &

      Shat2(3,2,ep)=
     &  - Cv(cc1 + C234,ep)
     &  + Cv(cc2 + C123,ep)
     &

      Shat2(1,3,ep)=
     &  + Cv(cc2 + C134,ep)
     &  - Cv(cc2 + C234,ep)
     &

      Shat2(2,3,ep)=
     &  + Cv(cc2 + C124,ep)
     &  - Cv(cc2 + C234,ep)
     &

      Shat2(3,3,ep)=
     &  - Cv(cc2 + C234,ep)
     &

      Shat1(1,ep) =
     &  + Cv(cc0 + C134,ep)
     &  - Cv(cc0 + C234,ep)
     &

      Shat1(2,ep) =
     &  + Cv(cc0 + C124,ep)
     &  - Cv(cc0 + C234,ep)
     &

      Shat1(3,ep) =
     &  + Cv(cc0 + C123,ep)
     &  - Cv(cc0 + C234,ep)
     &