01 function E=holmberg2(S,a,b,Q1,Q2)
02
03
04
05
06
07
08
09
10 Q1=(Q1+Q1')/2;
11 Q2=(Q2+Q2')/2;
12 E=1/4*(4*holmquist2(S,Q1,Q2)*(a'*S*b)+8*holmquist2(S,Q1,Q2)*(a'*S*b)*(a'*S*a)+...
13 4*holmquist2(S,Q1,Q2)*(a'*S*b)*(a'*S*a)^2+4*holmquist1(S,Q1)*(2*a'*S*Q2*S*b)+...
14 8*holmquist1(S,Q1)*(2*a'*S*Q2*S*b)*(a'*S*a)+4*holmquist1(S,Q1)*(2*a'*S*Q2*S*b)*(a'*S*a)^2-...
15 2*holmquist1(S,Q1)*(a'*S*b)*(2*a'*S*Q2*S*a)-2*holmquist1(S,Q1)*(a'*S*b)*(2*a'*S*Q2*S*a)*(a'*S*a)+...
16 4*holmquist1(S,Q2)*(2*a'*S*Q1*S*b)+8*holmquist1(S,Q2)*(2*a'*S*Q1*S*b)*(a'*S*a)+...
17 4*holmquist1(S,Q2)*(2*a'*S*Q1*S*b)*(a'*S*a)^2+4*1*(4*a'*(S*Q1*S*Q2*S+S*Q2*S*Q1*S)*b)+...
18 8*1*(4*a'*(S*Q1*S*Q2*S+S*Q2*S*Q1*S)*b)*(a'*S*a)+4*1*(4*a'*(S*Q1*S*Q2*S+S*Q2*S*Q1*S)*b)*(a'*S*a)^2-...
19 2*1*(2*a'*S*Q1*S*b)*(2*a'*S*Q2*S*a)-2*1*(2*a'*S*Q1*S*b)*(2*a'*S*Q2*S*a)*(a'*S*a)-...
20 2*holmquist1(S,Q2)*(a'*S*b)*(2*a'*S*Q1*S*a)-2*holmquist1(S,Q2)*(a'*S*b)*(2*a'*S*Q1*S*a)*(a'*S*a)-...
21 2*1*(2*a'*S*Q2*S*b)*(2*a'*S*Q1*S*a)-2*1*(2*a'*S*Q2*S*b)*(2*a'*S*Q1*S*a)*(a'*S*a)+...
22 3*1*(a'*S*b)*(2*a'*S*Q1*S*a)*(2*a'*S*Q2*S*a)-2*1*(a'*S*b)*(4*a'*(S*Q1*S*Q2*S+S*Q2*S*Q1*S)*a)-...
23 2*1*(a'*S*b)*(4*a'*(S*Q1*S*Q2*S+S*Q2*S*Q1*S)*a)*(a'*S*a))/(1+(a'*S*a))^(5/2);
24
25
26
27
28
29
30