femwind/hexa_sym.m

   1 t=[1 1 1];
   2 % t=sym('t',[1,3])
   3 s =[
   4     -1    -1    -1
   5     -1    -1     1
   6     -1     1    -1
   7     -1     1     1
   8      1    -1    -1
   9      1    -1     1
  10      1     1    -1
  11      1     1     1];
  12 s = s .* (ones(8,1)*t)
  13 x=sym('x',[3,1]);
  14 b=sym('b',[8,1]);
  15 g=sym('g',[8,3]);
  16 for k=1:8
  17     b(k) = (1+x(1)/s(k,1))*(1+x(2)/s(k,2))*(1+x(3)/s(k,3));
  18     for l=1:3
  19         g(k,l)=diff(b(k),x(l));
  20     end
  21 end
  22 K=int(int(int(g*g.',x(1),-t(1),t(1)),x(2),-t(2),t(2)),x(3),-t(3),t(3))