Items tagged with matrix

I have a system of equations e.g.

A^2+B*A+C=0

where A,B,C are Matrices and I want to solve for A.

Sure I can write every equations in brakets [..=0], but isn'T it possible to just use the matrix notation?

I've been trying to numerically solve (and plot) this equation. Maple tells me that some matrix is singular - I have no idea, what I can do.

eq := diff(y(x), `$`(x, 3))-(diff(y(x), x))*(diff(y(x), x))+1 = 0;

cond := (D(y))(0) = 0, (D(y))(1) = 1, ((D@@2)(y))(0) = 0

de := dsolve({cond, eq}, y(x), numeric);

Error, (in dsolve/numeric/bvp) matrix is singular

Hi,

 

I am trying to solve a simple system of the form AX=0, where A is a N*N matrix, X is an N*1 vector (and the right-hand side of the equation is an N*1 vector of zeros, I apologize for the inexact notation). The difficulty comes from the fact that the values of A are parameterized by 2*N parameters (that I will write as the 2*N vector P), and I would like to get a solution in the form X=f(P).

 

One solution is to try to use LinearAlgebra[LinearSolve], but it only returns the trivial solution X=0, which I am not interested in.

Another solution is to compute analytically the Moore-Penrose pseudoinverse Ag of A, as the general solution is of the form

(I - Ag A)f ;

where f is a vector of free parameters. However, even for a small matrix size (N=4), Maple is still computing after 3 hours on my (fairly powerful) machine, and it is taking more and more memory over time. As the results are polynomial/rational equations in the parameters P, I was actually expecting Maple to be more powerful than other softwares, but for this particular problem, Matlab's symbolic toolbox (muPAD) gives quick solutions until N=6. I need, in the end, to solve additional polynomial/rational equations that are derived from the solutions X=f(P), where Matlab fails. This is why I would really like to be able to solve the above-mentioned problem AX=0 with Maple in order to try to solve the subsequent step of the problem (polynomial system) with Maple.

 

Any suggestions on how to do this would be highly appreciated! Thank you very much for your time and help.

 

Laureline

When I print 2 matrix it always try to multiply them. How to force just to show?

I want to create a matrix (B) from entries of other matrices (A) with a helper-function (helper). The helper function is defined such that it returns a certain matrix depending on the index variables. This is necessary because the inner matrices are constructed with another function.

Since the helper-function returns matrices, the big matrix is of datatype=matrix. Unfortunately, creating the big matrix with the correct size and forcing the datatype=float, does not yield the desired result. However, the manual definition using the constructor with a list of matrices does create the desired matrix.

How do I resolve a matrix of matrices?

Note: I know that I could write a convert function that copies the entries to a corresponding matrix, though this seems to be unnecessary effort to me.

This might not be minimal but shows the issue. (Compare B and test)

MWE_matrix_of_matrices.mw

restart;
with(LinearAlgebra):

size_A := 2;
size_B := 3;

2

 

3

(1)

helper2 := proc(i::integer,j::integer);
  if i=j then
    a;
  elif i=j-1 or i=j+1 then
    b;
  else
    c;
  end if;
end proc:

helper3 := proc(i::integer,j::integer);
  if i=j then
    Matrix(size_A,size_A,helper2);
  elif i=j-1 or i=j+1 then
    -IdentityMatrix(size_A);
  else
    Matrix(size_A);
  end if;
end proc:

A := Matrix(size_A, size_A, helper2);
B := Matrix(size_B, size_B, helper3);
B := Matrix(size_B,size_B, helper3, datatype = float);
B := Matrix(size_B*size_A, size_B*size_A,[Matrix(size_B,size_B,helper3)], datatype = float)

A := Matrix(2, 2, {(1, 1) = a, (1, 2) = b, (2, 1) = b, (2, 2) = a})

 

B := Matrix(3, 3, {(1, 1) = Matrix(2, 2, {(1, 1) = a, (1, 2) = b, (2, 1) = b, (2, 2) = a}), (1, 2) = Matrix(2, 2, {(1, 1) = -1, (1, 2) = 0, (2, 1) = 0, (2, 2) = -1}), (1, 3) = Matrix(2, 2, {(1, 1) = 0, (1, 2) = 0, (2, 1) = 0, (2, 2) = 0}), (2, 1) = Matrix(2, 2, {(1, 1) = -1, (1, 2) = 0, (2, 1) = 0, (2, 2) = -1}), (2, 2) = Matrix(2, 2, {(1, 1) = a, (1, 2) = b, (2, 1) = b, (2, 2) = a}), (2, 3) = Matrix(2, 2, {(1, 1) = -1, (1, 2) = 0, (2, 1) = 0, (2, 2) = -1}), (3, 1) = Matrix(2, 2, {(1, 1) = 0, (1, 2) = 0, (2, 1) = 0, (2, 2) = 0}), (3, 2) = Matrix(2, 2, {(1, 1) = -1, (1, 2) = 0, (2, 1) = 0, (2, 2) = -1}), (3, 3) = Matrix(2, 2, {(1, 1) = a, (1, 2) = b, (2, 1) = b, (2, 2) = a})})

 

Error, (in Matrix) unable to store 'Matrix(2, 2, {(1, 1) = a, (1, 2) = b, (2, 1) = b, (2, 2) = a})' when datatype=float[8]

 

Error, (in Matrix) unable to store 'Matrix(2, 2, {(1, 1) = a, (1, 2) = b, (2, 1) = b, (2, 2) = a})' when datatype=float[8]

 

test := Matrix(4, 4, [
                [Matrix([[1,2],[0,9]]), Matrix([[3,6],[0,9]])],
                [Matrix([[3,4],[7,8]]), Matrix([[7,6],[5,5]])]
               ]); # is converted to a matrix of floats

test := Matrix(4, 4, {(1, 1) = 1, (1, 2) = 2, (1, 3) = 3, (1, 4) = 6, (2, 1) = 0, (2, 2) = 9, (2, 3) = 0, (2, 4) = 9, (3, 1) = 3, (3, 2) = 4, (3, 3) = 7, (3, 4) = 6, (4, 1) = 7, (4, 2) = 8, (4, 3) = 5, (4, 4) = 5})

(2)


Hi

We know determinant of a square matrix A[ij] (i,j ∈ {1,2,3}) is equal to the following expression

det(A) = 1/6 * e[ijk] * e[pqr] * A[ip] * A[jq] * A[kr] 

in which e[ijk] is a third order Tensor (Permutation notation or Levi-Civita symbol) and has a simple form as follows:

e[mnr] = 1/2 * (m-n) * (n-r) * (r-m).

The (i,j) minor of A, denoted Mij, is the determinant of the (n − 1)×(n − 1) matrix that results from deleting row i and column j of A. The cofactor matrix of A is the n×n matrix C whose (i, j) entry is the (i, j) cofactor of A,

C[ij]= -1 i+j * M[ij]

A-1=CT/det(A)

The general form of Levi-Civita symbol is as bellow:

 

I want to write a program for finding inverse of (NxN) matrix:

N=2 →

restart;
N := 2:
with(LinearAlgebra):
f := (1/2)*(sum(sum(sum(sum((m-n)*(p-q)*A[m, p]*A[n, q], q = 1 .. 2), p = 1 .. 2), n = 1 .. 2), m = 1 .. 2)):
A := Matrix(N, N, proc (i, j) options operator, arrow; evalf((37*i^2+j^3)/(2*i+4*j)) end proc):
f/Determinant(A);

N=3 →

restart;
N := 3:
with(LinearAlgebra):
f := (1/24)*(sum(sum(sum(sum(sum(sum((m-n)*(n-r)*(r-m)*(p-q)*(q-z)*(z-p)*A[m, p]*A[n, q]*A[r, z], m = 1 .. N), n = 1 .. N), r = 1 .. N), p = 1 .. N), q = 1 .. N), z = 1 .. N)):
A := Matrix(N, N, proc (i, j) options operator, arrow; 10*i^2/(20*i+j) end proc):
f/Determinant(A);

The results of above programs are equal to 1 and validation of method is observed.

If we can write the general form of determinant then we can find the inverse of any square non-singular matrices.

Now I try to write the mentioned program.

restart;
with(linalg):
N := 7:
Digits := 40:
e := product(product(signum(a[j]-a[i]), j = i+1 .. N), i = 1 .. N):
ML := product(A[a[k], b[k]], k = 1 .. N):
s[0] := e*subs(`$`(a[q] = b[q], q = 1 .. N), e)*ML:
for i to N do
s[i] := sum(sum(s[i-1], a[i] = 1 .. N), b[i] = 1 .. N)
end do:
A := Matrix(N, N, proc (i, j) options operator, arrow; evalf((3*i+j)/(i+2*j)) end proc): # arbitrary matrix
CN:=simplify(s[N]/det(A));

Therefore det(A)= CN-1 * e[a1,a2,..an] * e [b1,b2,.., bn] * A[a1,b1] * A[a2,b2] * ... * A[an,bn].

The correction coefficient is CN(for N)/CN(for N-1) = N!/(N-1)! =N.

restart:
with(linalg): N := 4: Digits := 20:
e := product(product(signum(a[j]-a[i]), j = i+1 .. N), i = 1 .. N):
ML := product(A[a[k], b[k]], k = 1 .. N):
s[0] := e*subs(`$`(a[q] = b[q], q = 1 .. N), e)*ML:
for r to N do s[r] := sum(sum(s[r-1], a[r] = 1 .. N), b[r] = 1 .. N) end do:
A := Matrix(N, N, proc (i, j) options operator, arrow; evalf((3*i+2*j)/(i+2*j)) end proc):
DET:=S[N]:
for x to N do for y to N do
e := product(product(signum(a[j]-a[i]), j = i+1 .. N-1), i = 1 .. N-1):
ML := product(AA[a[k], b[k]], k = 1 .. N-1):
S[0, x, y] := e*subs(`$`(a[q] = b[q], q = 1 .. N-1), e)*ML:
for r to N-1 do S[r, x, y] := sum(sum(S[r-1, x, y], a[r] = 1 .. N-1), b[r] = 1 .. N-1) end do:
f[y, x] := (-1)^(x+y)*subs(seq(seq(AA[t, u] = delrows(delcols(A, y .. y), x .. x)[t, u], t = 1 .. N-1), u = 1 .. N-1), S[N-1, x, y])
end do: end do:
Matrix(N, N, f)/(DET)*(24/6);
A^(-1);

CN for N=4 and N=3 is 24 and 6 respectively i.e. CN(4)/CN(3)=24/6.

When I use bellow procedure the error "(in S) bad index into Matrix" is occurred.

Please help me to write this algorithm by using procedure.

Thank you 

restart; with(linalg): Digits := 40: n := 7:
S := proc (N) local e, ML, s, i:
e := product(product(signum(a[j]-a[i]), j = i+1 .. N), i = 1 .. N):
ML := product(A[a[k], b[k]], k = 1 .. N):
s[0] := e*subs(`$`(a[q] = b[q], q = 1 .. N), e)*ML:
for i to N do s[i] := sum(sum(s[i-1], a[i] = 1 .. N), b[i] = 1 .. N) end do
end proc:
A := Matrix(n, n, proc (i, j) options operator, arrow; evalf((3*i+j)/(i+2*j)) end proc): # arbitrary matrix
CN := simplify(S(n)/det(A))

How can I use a matrix in timeseriesanalysis?

 

m3 := Import("d:\\bok6.xlsx")

 

Areachart works fine:

AreaChart([m3], color = "pink" .. "niagara Navy", title = "Tidslinje Tilluft og avtrekk")

 

But when matrix m3 is entered in timeseriesanalysis it gives a error:

ts1 := TimeSeries([m3], startdate = "2016-08-1", header = "Tilluft VAV setpunkt", enddate = "2016-08-01");

Error, (in TimeSeriesAnalysis:-TimeSeries) invalid input: too many and/or wrong type of arguments passed to TimeSeries:-ModuleCopy; first unused argument is [Array(1..4, 1..1, {(1, 1) = 155.0, (2, 1) = 142.0, (3, 1) = 133.0, (4, 1) = 122.0})]

 

 

I have loaded the following packages:

with*matrix  with*linearalgebra  with(plots)  with(Statistics)  with(TimeSeriesAnalysis)

 

Regards

 

Kristian

I defined a procedure, OrderB(0,b,c), that essentially determines the order of a group that is defined by 2 parameters (the first parameter is always zero), so the output is an integer and the procedure is defined for every integer value of both parameters. I wanted to structure the outputs into a matrix with columns that represent one parameter and rows that represent the other. However, whenever I try to apply a method of defining the matrix, Maple returns the values of the procedure and then gives me an empty matrix.

>f:=proc (i, j) -> OrderB(0, i, j+3);
>Matrix(3, f);
                               1
                               4
                               1
                               1
                               4
                               1
                               1
                               4
                               1
               [NULL,NULL...]


I tried a second method where I defined a set s such that:

s:={(0,0)=OrderB(0,0,0) , (0,1)=OrderB(0,0,1).....}

but upon execution, maple returns:

1
4
1
....
s:={(0,0)=( ), (0,1)=( )...}

Please Help. I don't know what I'm mising in the code that is keeping Maple from putting the outputs of my procedure into the matrix.

Hello,

I cant find solution how to create matrix form from equations of motion. Equations looks like this:

My equations ar much more complicated and one of them looks something like this:

 http://i63.tinypic.com/21c5ctk.png

and I want form like this:

I tried to do it using the Generate Matrix but it does not work as I expected. How can you get this form?

How do I get the Matrix to recalculate?

restart

a := 5

5

(1)

``

M := simplify(Matrix(2, 2, {(1, 1) = a, (1, 2) = 2*a, (2, 1) = 3*a, (2, 2) = a^2}))

M := Matrix(2, 2, {(1, 1) = 5, (1, 2) = 10, (2, 1) = 15, (2, 2) = 25})

(2)

``

unassign('a')

a

a

(3)

 

M

Matrix(2, 2, {(1, 1) = 5, (1, 2) = 10, (2, 1) = 15, (2, 2) = 25})

(4)

expand(M)

Matrix(2, 2, {(1, 1) = 5, (1, 2) = 10, (2, 1) = 15, (2, 2) = 25})

(5)

``

 

Download Re-evaluate_a_Matrix.mw

I am using a stoiciometric matrix to generate a system of differential equations (this makes them easier to check). There are 13 chemical species and 16 reactions so I need to make and display vectors of length 13 and 16; as well as a 13x16 matrix.

When I make Vectors of this size they don't get displayed; I just get information about the vector (Fortran order etc) and this is the same for the matrix.

How do i make/edit/display Matrixes nd Vectors of this size?

Hi, here is the code I used to try to generate all invertible 4x4 matrices over the finite field F_2 = {0.1}. However, when I look at the elements of GROUP (see below) all the elements are 4x4 matrices with a 2 in each entry. I don't know why this is?

Also, I need help converting all of the invertible 4x4 matrices in the following way: I want the 4x4 matrices to each be written as a string of length 16 with no spaces, commas or brackets. So for example the matrix

a b c d

e f g h

i j k l

m n o p

becomes abcdefghijklmnop

restart:
with(LinearAlgebra):
PRIME := 2;
2

# the group of invertible 4 x 4 matrices over the field F_p

GROUP := []:
M := Matrix([[a,b,c,d],[x,f,g,h],[y,j,k,l],[m,n,o,p]]):
for a from 0 to PRIME-1 do
for b from 0 to PRIME-1 do
for c from 0 to PRIME-1 do
for d from 0 to PRIME-1 do
for x from 0 to PRIME-1 do
for f from 0 to PRIME-1 do
for g from 0 to PRIME-1 do
for h from 0 to PRIME-1 do
for y from 0 to PRIME-1 do
for j from 0 to PRIME-1 do
for k from 0 to PRIME-1 do
for l from 0 to PRIME-1 do
for m from 0 to PRIME-1 do
for n from 0 to PRIME-1 do
for o from 0 to PRIME-1 do
for p from 0 to PRIME-1 do
if Determinant(M) mod PRIME <> 0 then
GROUP := [ op(GROUP), M ]
fi
od od od od od od od od od od od od od od od od:


nops(GROUP);
20160

GROUP;


----------------------------------------------------------------------------------------------------------------------
Introduction

I have a matrix (named DC on the piece of code below) all the elements of which are complex numbers a+b*I with a and b floating point numbers.
I want to obtain the real part of DC.

Obviously, if you do something like :
DC := Matrix(2,2, [1.0+1.0*I, 1.0-1.0*I, -1.0+1.0*I, -1.0-1.0*I];
Re~(DC);

the result corresponds to the desired matrix

----------------------------------------------------------------------------------------------------------------------
Context

In fact this matrix DC comes from some computations  described in the piece of code below

# Purpose :
# Given N points in a plane (here  in [0,1]X[0,1]), compute the matrix DX of distances between these points
#
# Example : if Pi and Pj are two such points, DX[i,j]=DX[j,i] denotes the Euclidian distance between Pi and Pj
#
# As I did not be able to find any single function in MAPLE that would construct DX , I proceed that way :
#   1/ let X the (N,2) matrix that contains the coordinates of the N points
#   2/ I represent these  N points as N complex numbers (vector C)
#   3/ I construct the (N,N) matrix MC = <C | C …..| C>
#   4/ I put MC = C – Transpose(C) :
#   5/ I take the norm DX of each elements of DC : DX := abs~(DC)
#       At this point, DX should contain the desired distances
#       But, due to floating point arithmetics, each element of DX writes a+0.*I where a is some floating point number)
#   6/ Last stage : execute Re~(DX)
#

with(Statistics):
with(LinearAlgebra):

N := 4:
X := Matrix(N,2, convert(Sample(Uniform(0,1), 2*N), list)):  #just an example

C   := X[..,1] +~ X[..,2] *~I;
MC := Multiply(C, Vector[row](N, 1));
DC := MC - Transpose(MC);
DX := abs~(DC);
Re~(DX)


----------------------------------------------------------------------------------------------------------------------
 My observations : 

1/ Maple 2015, Windows XP, 64 bytes
Re~(DX) returns DX and does not remove the imaginary (0.*I) part
But  Matrix(N, N, Re~(convert(DC, list))) does (which is a satisfactory, even if not clever, stopgap)

Why (it is just a question to help me to understand correctly how MAPLE proceeds) Re~(DC) does not (seem) to work here ?


2/ Maple 2015.2, Mac OS X El Capitan
DX := abs~(DC) gives me this strange result :
If (for instance) DC[i, j] = -1 – 2*I, DX[i, j] = +1 + 2*I
According to the compatibility problems between Maple 2015.1 and “El Capitan” (fixed from February), could it remain a few other problems ?

Last but not lesat : Did I do any syntax error ?


 I look forward to your responses

Hi everybody!

I am trying to find explicitely the relations between the columns of a matrix

of non-maximal rank. For example, if I have the matrix

M := Matrix([<1,2,3>, <2,4,6>, <5,6,7>]);

I would like that Maple finds that the second column is twice the first one: v_2 = 2*v_1.

How can I do?

hi.please see attached file below and help me.one problem is apply differential operator on matrix and then caclute 3D integral?

maple2.mw

restart; x = zz/L; y = (2*r-b)/a; z = alpha/Pi-1; L := .1; a := 0.1e-1; b := .11; E; 207*10^9; upsilon := .3

NN1 := -((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1); NN2 := ((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1); NN3 := -((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1); NN4 := ((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1); NN5 := ((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1); NN6 := -((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1); NN7 := ((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1); NN8 := -((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1); NN9 := ((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1); NN10 := -((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1); NN11 := ((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1); NN12 := -((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1); NN13 := -((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1); NN14 := ((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1); NN15 := -((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1); NN16 := ((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1)

((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1)

(1)

``

 

N := Matrix([[NN1, 0, 0, NN2, 0, 0, NN3, 0, 0, NN4, 0, 0, NN5, 0, 0, NN6, 0, 0, NN7, 0, 0, NN8, 0, 0, NN9, 0, 0, NN10, 0, 0, NN11, 0, 0, NN12, 0, 0, NN13, 0, 0, NN14, 0, 0, NN15, 0, 0, NN16, 0, 0], [0, NN1, 0, 0, NN2, 0, 0, NN3, 0, 0, NN4, 0, 0, NN5, 0, 0, NN6, 0, 0, NN7, 0, 0, NN8, 0, 0, NN9, 0, 0, NN10, 0, 0, NN11, 0, 0, NN12, 0, 0, NN13, 0, 0, NN14, 0, 0, NN15, 0, 0, NN16, 0], [0, 0, NN1, 0, 0, NN2, 0, 0, NN3, 0, 0, NN4, 0, 0, NN5, 0, 0, NN6, 0, 0, NN7, 0, 0, NN8, 0, 0, NN9, 0, 0, NN10, 0, 0, NN11, 0, 0, NN12, 0, 0, NN13, 0, 0, NN14, 0, 0, NN15, 0, 0, NN16]])

RTABLE(18446744074182475774, anything, Matrix, rectangular, Fortran_order, [], 2, 1 .. 3, 1 .. 48)

(2)

"Q:=Matrix([[(2/(a))*(&PartialD;)/(&PartialD; y) , 0,0],[2/(a*y+b),2/(a*y+b)*1/(Pi)(&PartialD;)/(&PartialD;z ) ,0],[0,0,1/(L)*(&PartialD;)/(&PartialD; x)],[2/(a*y+b)*1/(Pi)(&PartialD;)/(&PartialD;z ),2/(a)(&PartialD;)/(&PartialD;y)-2/(a*y+b),0],[1/(L)*(&PartialD;)/(&PartialD; x),0,(2/(a))*(&PartialD;)/(&PartialD; y)],[0,1/(L)*(&PartialD;)/(&PartialD; x),2/(a*y+b)*1/(Pi)(&PartialD;)/(&PartialD;z )]])"

Error, invalid derivative

"Q:=Matrix([[(2/a)*(&PartialD;)/(&PartialD;y) , 0,0],[2/(a*y+b),2/(a*y+b)*1/Pi(&PartialD;)/(&PartialD;z ) ,0],[0,0,1/L*(&PartialD;)/(&PartialD; x)],[2/(a*y+b)*1/Pi(&PartialD;)/(&PartialD;z ),2/a(&PartialD;)/(&PartialD;y)-2/(a*y+b),0],[1/L*(&PartialD;)/(&PartialD; x),0,(2/a)*(&PartialD;)/(&PartialD; y)],[0,1/L*(&PartialD;)/(&PartialD; x),2/(a*y+b)*1/Pi(&PartialD;)/(&PartialD;z )]])"

 

NULL

Q := Matrix([[2*Y/a, 0, 0], [2/(a*y+b), 2*Z/((a*y+b)*Pi), 0], [0, 0, X/L], [2*Z/((a*y+b)*Pi), 2*Y/a-2/(a*y+b), 0], [X/L, 0, 2*Y/a], [0, X/L, 2*Z/((a*y+b)*Pi)]])

Matrix([[0.2e3*Y, 0, 0], [2/(0.1e-1*y+.11), 2*Z/((0.1e-1*y+.11)*Pi), 0], [0, 0, 0.1e2*X], [2*Z/((0.1e-1*y+.11)*Pi), 0.2e3*Y-2/(0.1e-1*y+.11), 0], [0.1e2*X, 0, 0.2e3*Y], [0, 0.1e2*X, 2*Z/((0.1e-1*y+.11)*Pi)]])

(3)

````

"Y :=(&PartialD;)/(&PartialD; y):X:=(&PartialD;)/(&PartialD; x):Z:=(&PartialD;)/(&PartialD; z):"

Error, Got internal error in Typesetting:-Parse : "invalid subscript selector"

"Y :=(&PartialD;)/(&PartialD; y):X:=(&PartialD;)/(&PartialD; x):Z:=(&PartialD;)/(&PartialD; z):"

 

0

(4)

````

B := Q.N

RTABLE(18446744074182476230, anything, Matrix, rectangular, Fortran_order, [], 2, 1 .. 6, 1 .. 48)

(5)

NULL

Vector(4, {(1) = ` 6 x 48 `*Matrix, (2) = `Data Type: `*anything, (3) = `Storage: `*rectangular, (4) = `Order: `*Fortran_order})

(6)

d := (1-upsilon)/(1-2*upsilon); e := upsilon/(1-2*upsilon); DD := E*Matrix([[d, e, e, 0, 0, 0], [e, d, e, 0, 0, 0], [e, e, d, 0, 0, 0], [0, 0, 0, 1/2, 0, 0], [0, 0, 0, 0, 1/2, 0], [0, 0, 0, 0, 0, 1/2]])/(1+upsilon)

Matrix([[1.346153846*E, .5769230769*E, .5769230769*E, 0, 0, 0], [.5769230769*E, 1.346153846*E, .5769230769*E, 0, 0, 0], [.5769230769*E, .5769230769*E, 1.346153846*E, 0, 0, 0], [0, 0, 0, .3846153846*E, 0, 0], [0, 0, 0, 0, .3846153846*E, 0], [0, 0, 0, 0, 0, .3846153846*E]])

(7)

T := Transpose(B).DD.B

Transpose(Matrix(6, 48, {(1, 1) = -0.2e3*Y*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (1, 2) = 0., (1, 3) = 0., (1, 4) = 0.2e3*Y*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (1, 5) = 0., (1, 6) = 0., (1, 7) = -0.2e3*Y*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (1, 8) = 0., (1, 9) = 0., (1, 10) = 0.2e3*Y*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (1, 11) = 0., (1, 12) = 0., (1, 13) = 0.2e3*Y*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (1, 14) = 0., (1, 15) = 0., (1, 16) = -0.2e3*Y*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (1, 17) = 0., (1, 18) = 0., (1, 19) = 0.2e3*Y*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (1, 20) = 0., (1, 21) = 0., (1, 22) = -0.2e3*Y*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (1, 23) = 0., (1, 24) = 0., (1, 25) = 0.2e3*Y*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (1, 26) = 0., (1, 27) = 0., (1, 28) = -0.2e3*Y*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (1, 29) = 0., (1, 30) = 0., (1, 31) = 0.2e3*Y*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (1, 32) = 0., (1, 33) = 0., (1, 34) = -0.2e3*Y*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (1, 35) = 0., (1, 36) = 0., (1, 37) = -0.2e3*Y*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (1, 38) = 0., (1, 39) = 0., (1, 40) = 0.2e3*Y*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (1, 41) = 0., (1, 42) = 0., (1, 43) = -0.2e3*Y*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (1, 44) = 0., (1, 45) = 0., (1, 46) = 0.2e3*Y*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (1, 47) = 0., (1, 48) = 0., (2, 1) = -2*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (2, 2) = -2*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1)/((0.1e-1*y+.11)*Pi), (2, 3) = 0, (2, 4) = 2*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (2, 5) = 2*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1)/((0.1e-1*y+.11)*Pi), (2, 6) = 0, (2, 7) = -2*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (2, 8) = -2*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1)/((0.1e-1*y+.11)*Pi), (2, 9) = 0, (2, 10) = 2*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (2, 11) = 2*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1)/((0.1e-1*y+.11)*Pi), (2, 12) = 0, (2, 13) = 2*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (2, 14) = 2*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1)/((0.1e-1*y+.11)*Pi), (2, 15) = 0, (2, 16) = -2*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (2, 17) = -2*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1)/((0.1e-1*y+.11)*Pi), (2, 18) = 0, (2, 19) = 2*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (2, 20) = 2*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1)/((0.1e-1*y+.11)*Pi), (2, 21) = 0, (2, 22) = -2*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (2, 23) = -2*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1)/((0.1e-1*y+.11)*Pi), (2, 24) = 0, (2, 25) = 2*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (2, 26) = 2*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1)/((0.1e-1*y+.11)*Pi), (2, 27) = 0, (2, 28) = -2*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (2, 29) = -2*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1)/((0.1e-1*y+.11)*Pi), (2, 30) = 0, (2, 31) = 2*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (2, 32) = 2*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1)/((0.1e-1*y+.11)*Pi), (2, 33) = 0, (2, 34) = -2*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (2, 35) = -2*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1)/((0.1e-1*y+.11)*Pi), (2, 36) = 0, (2, 37) = -2*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (2, 38) = -2*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1)/((0.1e-1*y+.11)*Pi), (2, 39) = 0, (2, 40) = 2*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (2, 41) = 2*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1)/((0.1e-1*y+.11)*Pi), (2, 42) = 0, (2, 43) = -2*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (2, 44) = -2*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1)/((0.1e-1*y+.11)*Pi), (2, 45) = 0, (2, 46) = 2*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (2, 47) = 2*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1)/((0.1e-1*y+.11)*Pi), (2, 48) = 0, (3, 1) = 0., (3, 2) = 0., (3, 3) = -0.1e2*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (3, 4) = 0., (3, 5) = 0., (3, 6) = 0.1e2*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (3, 7) = 0., (3, 8) = 0., (3, 9) = -0.1e2*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (3, 10) = 0., (3, 11) = 0., (3, 12) = 0.1e2*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (3, 13) = 0., (3, 14) = 0., (3, 15) = 0.1e2*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (3, 16) = 0., (3, 17) = 0., (3, 18) = -0.1e2*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (3, 19) = 0., (3, 20) = 0., (3, 21) = 0.1e2*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (3, 22) = 0., (3, 23) = 0., (3, 24) = -0.1e2*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (3, 25) = 0., (3, 26) = 0., (3, 27) = 0.1e2*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (3, 28) = 0., (3, 29) = 0., (3, 30) = -0.1e2*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (3, 31) = 0., (3, 32) = 0., (3, 33) = 0.1e2*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (3, 34) = 0., (3, 35) = 0., (3, 36) = -0.1e2*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (3, 37) = 0., (3, 38) = 0., (3, 39) = -0.1e2*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (3, 40) = 0., (3, 41) = 0., (3, 42) = 0.1e2*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (3, 43) = 0., (3, 44) = 0., (3, 45) = -0.1e2*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (3, 46) = 0., (3, 47) = 0., (3, 48) = 0.1e2*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (4, 1) = -2*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1)/((0.1e-1*y+.11)*Pi), (4, 2) = -(0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (4, 3) = 0., (4, 4) = 2*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1)/((0.1e-1*y+.11)*Pi), (4, 5) = (0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (4, 6) = 0., (4, 7) = -2*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1)/((0.1e-1*y+.11)*Pi), (4, 8) = -(0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (4, 9) = 0., (4, 10) = 2*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1)/((0.1e-1*y+.11)*Pi), (4, 11) = (0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (4, 12) = 0., (4, 13) = 2*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1)/((0.1e-1*y+.11)*Pi), (4, 14) = (0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (4, 15) = 0., (4, 16) = -2*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1)/((0.1e-1*y+.11)*Pi), (4, 17) = -(0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (4, 18) = 0., (4, 19) = 2*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1)/((0.1e-1*y+.11)*Pi), (4, 20) = (0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (4, 21) = 0., (4, 22) = -2*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1)/((0.1e-1*y+.11)*Pi), (4, 23) = -(0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (4, 24) = 0., (4, 25) = 2*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1)/((0.1e-1*y+.11)*Pi), (4, 26) = (0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (4, 27) = 0., (4, 28) = -2*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1)/((0.1e-1*y+.11)*Pi), (4, 29) = -(0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (4, 30) = 0., (4, 31) = 2*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1)/((0.1e-1*y+.11)*Pi), (4, 32) = (0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (4, 33) = 0., (4, 34) = -2*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1)/((0.1e-1*y+.11)*Pi), (4, 35) = -(0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (4, 36) = 0., (4, 37) = -2*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1)/((0.1e-1*y+.11)*Pi), (4, 38) = -(0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (4, 39) = 0., (4, 40) = 2*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1)/((0.1e-1*y+.11)*Pi), (4, 41) = (0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (4, 42) = 0., (4, 43) = -2*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1)/((0.1e-1*y+.11)*Pi), (4, 44) = -(0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (4, 45) = 0., (4, 46) = 2*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1)/((0.1e-1*y+.11)*Pi), (4, 47) = (0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (4, 48) = 0., (5, 1) = -0.1e2*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (5, 2) = 0., (5, 3) = -0.2e3*Y*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (5, 4) = 0.1e2*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (5, 5) = 0., (5, 6) = 0.2e3*Y*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (5, 7) = -0.1e2*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (5, 8) = 0., (5, 9) = -0.2e3*Y*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (5, 10) = 0.1e2*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (5, 11) = 0., (5, 12) = 0.2e3*Y*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (5, 13) = 0.1e2*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (5, 14) = 0., (5, 15) = 0.2e3*Y*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (5, 16) = -0.1e2*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (5, 17) = 0., (5, 18) = -0.2e3*Y*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (5, 19) = 0.1e2*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (5, 20) = 0., (5, 21) = 0.2e3*Y*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (5, 22) = -0.1e2*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (5, 23) = 0., (5, 24) = -0.2e3*Y*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (5, 25) = 0.1e2*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (5, 26) = 0., (5, 27) = 0.2e3*Y*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (5, 28) = -0.1e2*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (5, 29) = 0., (5, 30) = -0.2e3*Y*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (5, 31) = 0.1e2*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (5, 32) = 0., (5, 33) = 0.2e3*Y*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (5, 34) = -0.1e2*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (5, 35) = 0., (5, 36) = -0.2e3*Y*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (5, 37) = -0.1e2*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (5, 38) = 0., (5, 39) = -0.2e3*Y*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (5, 40) = 0.1e2*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (5, 41) = 0., (5, 42) = 0.2e3*Y*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (5, 43) = -0.1e2*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (5, 44) = 0., (5, 45) = -0.2e3*Y*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (5, 46) = 0.1e2*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (5, 47) = 0., (5, 48) = 0.2e3*Y*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (6, 1) = 0., (6, 2) = -0.1e2*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (6, 3) = -2*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1)/((0.1e-1*y+.11)*Pi), (6, 4) = 0., (6, 5) = 0.1e2*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (6, 6) = 2*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1)/((0.1e-1*y+.11)*Pi), (6, 7) = 0., (6, 8) = -0.1e2*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (6, 9) = -2*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1)/((0.1e-1*y+.11)*Pi), (6, 10) = 0., (6, 11) = 0.1e2*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (6, 12) = 2*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1)/((0.1e-1*y+.11)*Pi), (6, 13) = 0., (6, 14) = 0.1e2*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (6, 15) = 2*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1)/((0.1e-1*y+.11)*Pi), (6, 16) = 0., (6, 17) = -0.1e2*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (6, 18) = -2*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1)/((0.1e-1*y+.11)*Pi), (6, 19) = 0., (6, 20) = 0.1e2*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (6, 21) = 2*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1)/((0.1e-1*y+.11)*Pi), (6, 22) = 0., (6, 23) = -0.1e2*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (6, 24) = -2*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1)/((0.1e-1*y+.11)*Pi), (6, 25) = 0., (6, 26) = 0.1e2*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (6, 27) = 2*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1)/((0.1e-1*y+.11)*Pi), (6, 28) = 0., (6, 29) = -0.1e2*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (6, 30) = -2*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1)/((0.1e-1*y+.11)*Pi), (6, 31) = 0., (6, 32) = 0.1e2*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (6, 33) = 2*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1)/((0.1e-1*y+.11)*Pi), (6, 34) = 0., (6, 35) = -0.1e2*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (6, 36) = -2*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1)/((0.1e-1*y+.11)*Pi), (6, 37) = 0., (6, 38) = -0.1e2*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (6, 39) = -2*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1)/((0.1e-1*y+.11)*Pi), (6, 40) = 0., (6, 41) = 0.1e2*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (6, 42) = 2*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1)/((0.1e-1*y+.11)*Pi), (6, 43) = 0., (6, 44) = -0.1e2*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (6, 45) = -2*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1)/((0.1e-1*y+.11)*Pi), (6, 46) = 0., (6, 47) = 0.1e2*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (6, 48) = 2*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1)/((0.1e-1*y+.11)*Pi)})).(Matrix(6, 48, {(1, 1) = -269.2307692*E*Y*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1)-1.153846154*E*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (1, 2) = -.3672806379*E*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (1, 3) = -5.769230769*E*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (1, 4) = 269.2307692*E*Y*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1)+1.153846154*E*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (1, 5) = .3672806379*E*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (1, 6) = 5.769230769*E*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (1, 7) = -269.2307692*E*Y*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1)-1.153846154*E*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (1, 8) = -.3672806379*E*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (1, 9) = -5.769230769*E*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (1, 10) = 269.2307692*E*Y*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1)+1.153846154*E*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (1, 11) = .3672806379*E*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (1, 12) = 5.769230769*E*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (1, 13) = 269.2307692*E*Y*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1)+1.153846154*E*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (1, 14) = .3672806379*E*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (1, 15) = 5.769230769*E*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (1, 16) = -269.2307692*E*Y*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1)-1.153846154*E*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (1, 17) = -.3672806379*E*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (1, 18) = -5.769230769*E*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (1, 19) = 269.2307692*E*Y*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1)+1.153846154*E*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (1, 20) = .3672806379*E*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (1, 21) = 5.769230769*E*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (1, 22) = -269.2307692*E*Y*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1)-1.153846154*E*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (1, 23) = -.3672806379*E*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (1, 24) = -5.769230769*E*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (1, 25) = 269.2307692*E*Y*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1)+1.153846154*E*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (1, 26) = .3672806379*E*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (1, 27) = 5.769230769*E*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (1, 28) = -269.2307692*E*Y*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1)-1.153846154*E*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (1, 29) = -.3672806379*E*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (1, 30) = -5.769230769*E*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (1, 31) = 269.2307692*E*Y*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1)+1.153846154*E*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (1, 32) = .3672806379*E*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (1, 33) = 5.769230769*E*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (1, 34) = -269.2307692*E*Y*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1)-1.153846154*E*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (1, 35) = -.3672806379*E*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (1, 36) = -5.769230769*E*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (1, 37) = -269.2307692*E*Y*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1)-1.153846154*E*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (1, 38) = -.3672806379*E*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (1, 39) = -5.769230769*E*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (1, 40) = 269.2307692*E*Y*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1)+1.153846154*E*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (1, 41) = .3672806379*E*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (1, 42) = 5.769230769*E*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (1, 43) = -269.2307692*E*Y*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1)-1.153846154*E*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (1, 44) = -.3672806379*E*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (1, 45) = -5.769230769*E*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (1, 46) = 269.2307692*E*Y*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1)+1.153846154*E*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (1, 47) = .3672806379*E*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (1, 48) = 5.769230769*E*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (2, 1) = -115.3846154*E*Y*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1)-2.692307692*E*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (2, 2) = -.8569881549*E*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (2, 3) = -5.769230769*E*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (2, 4) = 115.3846154*E*Y*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1)+2.692307692*E*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (2, 5) = .8569881549*E*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (2, 6) = 5.769230769*E*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (2, 7) = -115.3846154*E*Y*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1)-2.692307692*E*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (2, 8) = -.8569881549*E*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (2, 9) = -5.769230769*E*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (2, 10) = 115.3846154*E*Y*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1)+2.692307692*E*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (2, 11) = .8569881549*E*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (2, 12) = 5.769230769*E*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (2, 13) = 115.3846154*E*Y*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1)+2.692307692*E*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (2, 14) = .8569881549*E*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (2, 15) = 5.769230769*E*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (2, 16) = -115.3846154*E*Y*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1)-2.692307692*E*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (2, 17) = -.8569881549*E*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (2, 18) = -5.769230769*E*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (2, 19) = 115.3846154*E*Y*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1)+2.692307692*E*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (2, 20) = .8569881549*E*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (2, 21) = 5.769230769*E*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (2, 22) = -115.3846154*E*Y*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1)-2.692307692*E*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (2, 23) = -.8569881549*E*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (2, 24) = -5.769230769*E*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (2, 25) = 115.3846154*E*Y*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1)+2.692307692*E*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (2, 26) = .8569881549*E*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (2, 27) = 5.769230769*E*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (2, 28) = -115.3846154*E*Y*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1)-2.692307692*E*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (2, 29) = -.8569881549*E*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (2, 30) = -5.769230769*E*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (2, 31) = 115.3846154*E*Y*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1)+2.692307692*E*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (2, 32) = .8569881549*E*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (2, 33) = 5.769230769*E*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (2, 34) = -115.3846154*E*Y*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1)-2.692307692*E*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (2, 35) = -.8569881549*E*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (2, 36) = -5.769230769*E*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (2, 37) = -115.3846154*E*Y*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1)-2.692307692*E*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (2, 38) = -.8569881549*E*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (2, 39) = -5.769230769*E*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (2, 40) = 115.3846154*E*Y*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1)+2.692307692*E*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (2, 41) = .8569881549*E*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (2, 42) = 5.769230769*E*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (2, 43) = -115.3846154*E*Y*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1)-2.692307692*E*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (2, 44) = -.8569881549*E*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (2, 45) = -5.769230769*E*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (2, 46) = 115.3846154*E*Y*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1)+2.692307692*E*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (2, 47) = .8569881549*E*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (2, 48) = 5.769230769*E*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (3, 1) = -115.3846154*E*Y*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1)-1.153846154*E*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (3, 2) = -.3672806379*E*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (3, 3) = -13.46153846*E*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (3, 4) = 115.3846154*E*Y*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1)+1.153846154*E*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (3, 5) = .3672806379*E*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (3, 6) = 13.46153846*E*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (3, 7) = -115.3846154*E*Y*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1)-1.153846154*E*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (3, 8) = -.3672806379*E*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (3, 9) = -13.46153846*E*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (3, 10) = 115.3846154*E*Y*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1)+1.153846154*E*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (3, 11) = .3672806379*E*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (3, 12) = 13.46153846*E*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (3, 13) = 115.3846154*E*Y*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1)+1.153846154*E*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (3, 14) = .3672806379*E*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (3, 15) = 13.46153846*E*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (3, 16) = -115.3846154*E*Y*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1)-1.153846154*E*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (3, 17) = -.3672806379*E*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (3, 18) = -13.46153846*E*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (3, 19) = 115.3846154*E*Y*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1)+1.153846154*E*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (3, 20) = .3672806379*E*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (3, 21) = 13.46153846*E*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (3, 22) = -115.3846154*E*Y*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1)-1.153846154*E*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (3, 23) = -.3672806379*E*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (3, 24) = -13.46153846*E*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (3, 25) = 115.3846154*E*Y*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1)+1.153846154*E*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (3, 26) = .3672806379*E*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (3, 27) = 13.46153846*E*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (3, 28) = -115.3846154*E*Y*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1)-1.153846154*E*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (3, 29) = -.3672806379*E*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (3, 30) = -13.46153846*E*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (3, 31) = 115.3846154*E*Y*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1)+1.153846154*E*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (3, 32) = .3672806379*E*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (3, 33) = 13.46153846*E*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (3, 34) = -115.3846154*E*Y*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1)-1.153846154*E*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (3, 35) = -.3672806379*E*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (3, 36) = -13.46153846*E*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (3, 37) = -115.3846154*E*Y*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1)-1.153846154*E*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (3, 38) = -.3672806379*E*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (3, 39) = -13.46153846*E*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (3, 40) = 115.3846154*E*Y*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1)+1.153846154*E*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (3, 41) = .3672806379*E*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (3, 42) = 13.46153846*E*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (3, 43) = -115.3846154*E*Y*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1)-1.153846154*E*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (3, 44) = -.3672806379*E*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (3, 45) = -13.46153846*E*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (3, 46) = 115.3846154*E*Y*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1)+1.153846154*E*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (3, 47) = .3672806379*E*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (3, 48) = 13.46153846*E*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (4, 1) = -.2448537586*E*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (4, 2) = -.3846153846*E*(0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (4, 3) = 0., (4, 4) = .2448537586*E*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (4, 5) = .3846153846*E*(0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (4, 6) = 0., (4, 7) = -.2448537586*E*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (4, 8) = -.3846153846*E*(0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (4, 9) = 0., (4, 10) = .2448537586*E*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (4, 11) = .3846153846*E*(0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (4, 12) = 0., (4, 13) = .2448537586*E*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (4, 14) = .3846153846*E*(0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (4, 15) = 0., (4, 16) = -.2448537586*E*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (4, 17) = -.3846153846*E*(0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (4, 18) = 0., (4, 19) = .2448537586*E*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (4, 20) = .3846153846*E*(0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (4, 21) = 0., (4, 22) = -.2448537586*E*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (4, 23) = -.3846153846*E*(0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (4, 24) = 0., (4, 25) = .2448537586*E*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (4, 26) = .3846153846*E*(0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (4, 27) = 0., (4, 28) = -.2448537586*E*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (4, 29) = -.3846153846*E*(0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (4, 30) = 0., (4, 31) = .2448537586*E*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (4, 32) = .3846153846*E*(0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (4, 33) = 0., (4, 34) = -.2448537586*E*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (4, 35) = -.3846153846*E*(0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (4, 36) = 0., (4, 37) = -.2448537586*E*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (4, 38) = -.3846153846*E*(0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (4, 39) = 0., (4, 40) = .2448537586*E*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (4, 41) = .3846153846*E*(0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (4, 42) = 0., (4, 43) = -.2448537586*E*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (4, 44) = -.3846153846*E*(0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (4, 45) = 0., (4, 46) = .2448537586*E*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (4, 47) = .3846153846*E*(0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (4, 48) = 0., (5, 1) = -3.846153846*E*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (5, 2) = 0., (5, 3) = -76.92307692*E*Y*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (5, 4) = 3.846153846*E*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (5, 5) = 0., (5, 6) = 76.92307692*E*Y*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (5, 7) = -3.846153846*E*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (5, 8) = 0., (5, 9) = -76.92307692*E*Y*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (5, 10) = 3.846153846*E*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (5, 11) = 0., (5, 12) = 76.92307692*E*Y*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (5, 13) = 3.846153846*E*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (5, 14) = 0., (5, 15) = 76.92307692*E*Y*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (5, 16) = -3.846153846*E*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (5, 17) = 0., (5, 18) = -76.92307692*E*Y*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (5, 19) = 3.846153846*E*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (5, 20) = 0., (5, 21) = 76.92307692*E*Y*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (5, 22) = -3.846153846*E*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (5, 23) = 0., (5, 24) = -76.92307692*E*Y*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (5, 25) = 3.846153846*E*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (5, 26) = 0., (5, 27) = 76.92307692*E*Y*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (5, 28) = -3.846153846*E*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (5, 29) = 0., (5, 30) = -76.92307692*E*Y*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (5, 31) = 3.846153846*E*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (5, 32) = 0., (5, 33) = 76.92307692*E*Y*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (5, 34) = -3.846153846*E*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (5, 35) = 0., (5, 36) = -76.92307692*E*Y*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (5, 37) = -3.846153846*E*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (5, 38) = 0., (5, 39) = -76.92307692*E*Y*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (5, 40) = 3.846153846*E*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (5, 41) = 0., (5, 42) = 76.92307692*E*Y*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (5, 43) = -3.846153846*E*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (5, 44) = 0., (5, 45) = -76.92307692*E*Y*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (5, 46) = 3.846153846*E*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (5, 47) = 0., (5, 48) = 76.92307692*E*Y*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (6, 1) = 0., (6, 2) = -3.846153846*E*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (6, 3) = -.2448537586*E*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (6, 4) = 0., (6, 5) = 3.846153846*E*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (6, 6) = .2448537586*E*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (6, 7) = 0., (6, 8) = -3.846153846*E*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (6, 9) = -.2448537586*E*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (6, 10) = 0., (6, 11) = 3.846153846*E*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (6, 12) = .2448537586*E*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (6, 13) = 0., (6, 14) = 3.846153846*E*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (6, 15) = .2448537586*E*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (6, 16) = 0., (6, 17) = -3.846153846*E*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (6, 18) = -.2448537586*E*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (6, 19) = 0., (6, 20) = 3.846153846*E*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (6, 21) = .2448537586*E*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (6, 22) = 0., (6, 23) = -3.846153846*E*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (6, 24) = -.2448537586*E*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (6, 25) = 0., (6, 26) = 3.846153846*E*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (6, 27) = .2448537586*E*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (6, 28) = 0., (6, 29) = -3.846153846*E*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (6, 30) = -.2448537586*E*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (6, 31) = 0., (6, 32) = 3.846153846*E*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (6, 33) = .2448537586*E*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (6, 34) = 0., (6, 35) = -3.846153846*E*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (6, 36) = -.2448537586*E*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (6, 37) = 0., (6, 38) = -3.846153846*E*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (6, 39) = -.2448537586*E*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (6, 40) = 0., (6, 41) = 3.846153846*E*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (6, 42) = .2448537586*E*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (6, 43) = 0., (6, 44) = -3.846153846*E*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (6, 45) = -.2448537586*E*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (6, 46) = 0., (6, 47) = 3.846153846*E*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (6, 48) = .2448537586*E*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11)}))

(8)

S := (1/4)*a*Pi*L*(a*y+b)*T

Typesetting[delayDotProduct](0.7853981635e-3*(0.1e-1*y+.11), Transpose(Matrix(6, 48, {(1, 1) = -0.2e3*Y*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (1, 2) = 0., (1, 3) = 0., (1, 4) = 0.2e3*Y*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (1, 5) = 0., (1, 6) = 0., (1, 7) = -0.2e3*Y*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (1, 8) = 0., (1, 9) = 0., (1, 10) = 0.2e3*Y*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (1, 11) = 0., (1, 12) = 0., (1, 13) = 0.2e3*Y*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (1, 14) = 0., (1, 15) = 0., (1, 16) = -0.2e3*Y*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (1, 17) = 0., (1, 18) = 0., (1, 19) = 0.2e3*Y*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (1, 20) = 0., (1, 21) = 0., (1, 22) = -0.2e3*Y*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (1, 23) = 0., (1, 24) = 0., (1, 25) = 0.2e3*Y*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (1, 26) = 0., (1, 27) = 0., (1, 28) = -0.2e3*Y*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (1, 29) = 0., (1, 30) = 0., (1, 31) = 0.2e3*Y*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (1, 32) = 0., (1, 33) = 0., (1, 34) = -0.2e3*Y*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (1, 35) = 0., (1, 36) = 0., (1, 37) = -0.2e3*Y*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (1, 38) = 0., (1, 39) = 0., (1, 40) = 0.2e3*Y*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (1, 41) = 0., (1, 42) = 0., (1, 43) = -0.2e3*Y*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (1, 44) = 0., (1, 45) = 0., (1, 46) = 0.2e3*Y*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (1, 47) = 0., (1, 48) = 0., (2, 1) = -2*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (2, 2) = -2*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1)/((0.1e-1*y+.11)*Pi), (2, 3) = 0, (2, 4) = 2*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (2, 5) = 2*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1)/((0.1e-1*y+.11)*Pi), (2, 6) = 0, (2, 7) = -2*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (2, 8) = -2*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1)/((0.1e-1*y+.11)*Pi), (2, 9) = 0, (2, 10) = 2*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (2, 11) = 2*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1)/((0.1e-1*y+.11)*Pi), (2, 12) = 0, (2, 13) = 2*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (2, 14) = 2*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1)/((0.1e-1*y+.11)*Pi), (2, 15) = 0, (2, 16) = -2*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (2, 17) = -2*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1)/((0.1e-1*y+.11)*Pi), (2, 18) = 0, (2, 19) = 2*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (2, 20) = 2*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1)/((0.1e-1*y+.11)*Pi), (2, 21) = 0, (2, 22) = -2*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (2, 23) = -2*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1)/((0.1e-1*y+.11)*Pi), (2, 24) = 0, (2, 25) = 2*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (2, 26) = 2*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1)/((0.1e-1*y+.11)*Pi), (2, 27) = 0, (2, 28) = -2*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (2, 29) = -2*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1)/((0.1e-1*y+.11)*Pi), (2, 30) = 0, (2, 31) = 2*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (2, 32) = 2*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1)/((0.1e-1*y+.11)*Pi), (2, 33) = 0, (2, 34) = -2*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (2, 35) = -2*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1)/((0.1e-1*y+.11)*Pi), (2, 36) = 0, (2, 37) = -2*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (2, 38) = -2*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1)/((0.1e-1*y+.11)*Pi), (2, 39) = 0, (2, 40) = 2*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (2, 41) = 2*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1)/((0.1e-1*y+.11)*Pi), (2, 42) = 0, (2, 43) = -2*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (2, 44) = -2*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1)/((0.1e-1*y+.11)*Pi), (2, 45) = 0, (2, 46) = 2*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (2, 47) = 2*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1)/((0.1e-1*y+.11)*Pi), (2, 48) = 0, (3, 1) = 0., (3, 2) = 0., (3, 3) = -0.1e2*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (3, 4) = 0., (3, 5) = 0., (3, 6) = 0.1e2*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (3, 7) = 0., (3, 8) = 0., (3, 9) = -0.1e2*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (3, 10) = 0., (3, 11) = 0., (3, 12) = 0.1e2*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (3, 13) = 0., (3, 14) = 0., (3, 15) = 0.1e2*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (3, 16) = 0., (3, 17) = 0., (3, 18) = -0.1e2*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (3, 19) = 0., (3, 20) = 0., (3, 21) = 0.1e2*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (3, 22) = 0., (3, 23) = 0., (3, 24) = -0.1e2*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (3, 25) = 0., (3, 26) = 0., (3, 27) = 0.1e2*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (3, 28) = 0., (3, 29) = 0., (3, 30) = -0.1e2*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (3, 31) = 0., (3, 32) = 0., (3, 33) = 0.1e2*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (3, 34) = 0., (3, 35) = 0., (3, 36) = -0.1e2*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (3, 37) = 0., (3, 38) = 0., (3, 39) = -0.1e2*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (3, 40) = 0., (3, 41) = 0., (3, 42) = 0.1e2*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (3, 43) = 0., (3, 44) = 0., (3, 45) = -0.1e2*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (3, 46) = 0., (3, 47) = 0., (3, 48) = 0.1e2*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (4, 1) = -2*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1)/((0.1e-1*y+.11)*Pi), (4, 2) = -(0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (4, 3) = 0., (4, 4) = 2*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1)/((0.1e-1*y+.11)*Pi), (4, 5) = (0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (4, 6) = 0., (4, 7) = -2*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1)/((0.1e-1*y+.11)*Pi), (4, 8) = -(0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (4, 9) = 0., (4, 10) = 2*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1)/((0.1e-1*y+.11)*Pi), (4, 11) = (0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (4, 12) = 0., (4, 13) = 2*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1)/((0.1e-1*y+.11)*Pi), (4, 14) = (0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (4, 15) = 0., (4, 16) = -2*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1)/((0.1e-1*y+.11)*Pi), (4, 17) = -(0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (4, 18) = 0., (4, 19) = 2*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1)/((0.1e-1*y+.11)*Pi), (4, 20) = (0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (4, 21) = 0., (4, 22) = -2*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1)/((0.1e-1*y+.11)*Pi), (4, 23) = -(0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (4, 24) = 0., (4, 25) = 2*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1)/((0.1e-1*y+.11)*Pi), (4, 26) = (0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (4, 27) = 0., (4, 28) = -2*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1)/((0.1e-1*y+.11)*Pi), (4, 29) = -(0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (4, 30) = 0., (4, 31) = 2*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1)/((0.1e-1*y+.11)*Pi), (4, 32) = (0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (4, 33) = 0., (4, 34) = -2*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1)/((0.1e-1*y+.11)*Pi), (4, 35) = -(0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (4, 36) = 0., (4, 37) = -2*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1)/((0.1e-1*y+.11)*Pi), (4, 38) = -(0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (4, 39) = 0., (4, 40) = 2*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1)/((0.1e-1*y+.11)*Pi), (4, 41) = (0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (4, 42) = 0., (4, 43) = -2*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1)/((0.1e-1*y+.11)*Pi), (4, 44) = -(0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (4, 45) = 0., (4, 46) = 2*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1)/((0.1e-1*y+.11)*Pi), (4, 47) = (0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (4, 48) = 0., (5, 1) = -0.1e2*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (5, 2) = 0., (5, 3) = -0.2e3*Y*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (5, 4) = 0.1e2*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (5, 5) = 0., (5, 6) = 0.2e3*Y*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (5, 7) = -0.1e2*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (5, 8) = 0., (5, 9) = -0.2e3*Y*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (5, 10) = 0.1e2*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (5, 11) = 0., (5, 12) = 0.2e3*Y*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (5, 13) = 0.1e2*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (5, 14) = 0., (5, 15) = 0.2e3*Y*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (5, 16) = -0.1e2*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (5, 17) = 0., (5, 18) = -0.2e3*Y*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (5, 19) = 0.1e2*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (5, 20) = 0., (5, 21) = 0.2e3*Y*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (5, 22) = -0.1e2*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (5, 23) = 0., (5, 24) = -0.2e3*Y*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (5, 25) = 0.1e2*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (5, 26) = 0., (5, 27) = 0.2e3*Y*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (5, 28) = -0.1e2*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (5, 29) = 0., (5, 30) = -0.2e3*Y*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (5, 31) = 0.1e2*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (5, 32) = 0., (5, 33) = 0.2e3*Y*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (5, 34) = -0.1e2*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (5, 35) = 0., (5, 36) = -0.2e3*Y*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (5, 37) = -0.1e2*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (5, 38) = 0., (5, 39) = -0.2e3*Y*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (5, 40) = 0.1e2*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (5, 41) = 0., (5, 42) = 0.2e3*Y*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (5, 43) = -0.1e2*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (5, 44) = 0., (5, 45) = -0.2e3*Y*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (5, 46) = 0.1e2*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (5, 47) = 0., (5, 48) = 0.2e3*Y*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (6, 1) = 0., (6, 2) = -0.1e2*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (6, 3) = -2*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1)/((0.1e-1*y+.11)*Pi), (6, 4) = 0., (6, 5) = 0.1e2*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (6, 6) = 2*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1)/((0.1e-1*y+.11)*Pi), (6, 7) = 0., (6, 8) = -0.1e2*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (6, 9) = -2*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1)/((0.1e-1*y+.11)*Pi), (6, 10) = 0., (6, 11) = 0.1e2*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (6, 12) = 2*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1)/((0.1e-1*y+.11)*Pi), (6, 13) = 0., (6, 14) = 0.1e2*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (6, 15) = 2*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1)/((0.1e-1*y+.11)*Pi), (6, 16) = 0., (6, 17) = -0.1e2*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (6, 18) = -2*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1)/((0.1e-1*y+.11)*Pi), (6, 19) = 0., (6, 20) = 0.1e2*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (6, 21) = 2*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1)/((0.1e-1*y+.11)*Pi), (6, 22) = 0., (6, 23) = -0.1e2*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (6, 24) = -2*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1)/((0.1e-1*y+.11)*Pi), (6, 25) = 0., (6, 26) = 0.1e2*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (6, 27) = 2*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1)/((0.1e-1*y+.11)*Pi), (6, 28) = 0., (6, 29) = -0.1e2*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (6, 30) = -2*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1)/((0.1e-1*y+.11)*Pi), (6, 31) = 0., (6, 32) = 0.1e2*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (6, 33) = 2*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1)/((0.1e-1*y+.11)*Pi), (6, 34) = 0., (6, 35) = -0.1e2*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (6, 36) = -2*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1)/((0.1e-1*y+.11)*Pi), (6, 37) = 0., (6, 38) = -0.1e2*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (6, 39) = -2*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1)/((0.1e-1*y+.11)*Pi), (6, 40) = 0., (6, 41) = 0.1e2*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (6, 42) = 2*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1)/((0.1e-1*y+.11)*Pi), (6, 43) = 0., (6, 44) = -0.1e2*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (6, 45) = -2*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1)/((0.1e-1*y+.11)*Pi), (6, 46) = 0., (6, 47) = 0.1e2*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (6, 48) = 2*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1)/((0.1e-1*y+.11)*Pi)})).(Matrix(6, 48, {(1, 1) = -269.2307692*E*Y*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1)-1.153846154*E*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (1, 2) = -.3672806379*E*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (1, 3) = -5.769230769*E*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (1, 4) = 269.2307692*E*Y*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1)+1.153846154*E*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (1, 5) = .3672806379*E*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (1, 6) = 5.769230769*E*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (1, 7) = -269.2307692*E*Y*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1)-1.153846154*E*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (1, 8) = -.3672806379*E*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (1, 9) = -5.769230769*E*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (1, 10) = 269.2307692*E*Y*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1)+1.153846154*E*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (1, 11) = .3672806379*E*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (1, 12) = 5.769230769*E*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (1, 13) = 269.2307692*E*Y*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1)+1.153846154*E*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (1, 14) = .3672806379*E*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (1, 15) = 5.769230769*E*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (1, 16) = -269.2307692*E*Y*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1)-1.153846154*E*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (1, 17) = -.3672806379*E*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (1, 18) = -5.769230769*E*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (1, 19) = 269.2307692*E*Y*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1)+1.153846154*E*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (1, 20) = .3672806379*E*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (1, 21) = 5.769230769*E*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (1, 22) = -269.2307692*E*Y*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1)-1.153846154*E*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (1, 23) = -.3672806379*E*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (1, 24) = -5.769230769*E*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (1, 25) = 269.2307692*E*Y*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1)+1.153846154*E*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (1, 26) = .3672806379*E*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (1, 27) = 5.769230769*E*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (1, 28) = -269.2307692*E*Y*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1)-1.153846154*E*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (1, 29) = -.3672806379*E*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (1, 30) = -5.769230769*E*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (1, 31) = 269.2307692*E*Y*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1)+1.153846154*E*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (1, 32) = .3672806379*E*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (1, 33) = 5.769230769*E*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (1, 34) = -269.2307692*E*Y*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1)-1.153846154*E*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (1, 35) = -.3672806379*E*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (1, 36) = -5.769230769*E*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (1, 37) = -269.2307692*E*Y*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1)-1.153846154*E*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (1, 38) = -.3672806379*E*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (1, 39) = -5.769230769*E*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (1, 40) = 269.2307692*E*Y*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1)+1.153846154*E*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (1, 41) = .3672806379*E*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (1, 42) = 5.769230769*E*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (1, 43) = -269.2307692*E*Y*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1)-1.153846154*E*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (1, 44) = -.3672806379*E*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (1, 45) = -5.769230769*E*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (1, 46) = 269.2307692*E*Y*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1)+1.153846154*E*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (1, 47) = .3672806379*E*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (1, 48) = 5.769230769*E*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (2, 1) = -115.3846154*E*Y*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1)-2.692307692*E*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (2, 2) = -.8569881549*E*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (2, 3) = -5.769230769*E*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (2, 4) = 115.3846154*E*Y*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1)+2.692307692*E*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (2, 5) = .8569881549*E*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (2, 6) = 5.769230769*E*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (2, 7) = -115.3846154*E*Y*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1)-2.692307692*E*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (2, 8) = -.8569881549*E*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (2, 9) = -5.769230769*E*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (2, 10) = 115.3846154*E*Y*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1)+2.692307692*E*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (2, 11) = .8569881549*E*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (2, 12) = 5.769230769*E*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (2, 13) = 115.3846154*E*Y*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1)+2.692307692*E*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (2, 14) = .8569881549*E*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (2, 15) = 5.769230769*E*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (2, 16) = -115.3846154*E*Y*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1)-2.692307692*E*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (2, 17) = -.8569881549*E*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (2, 18) = -5.769230769*E*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (2, 19) = 115.3846154*E*Y*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1)+2.692307692*E*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (2, 20) = .8569881549*E*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (2, 21) = 5.769230769*E*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (2, 22) = -115.3846154*E*Y*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1)-2.692307692*E*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (2, 23) = -.8569881549*E*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (2, 24) = -5.769230769*E*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (2, 25) = 115.3846154*E*Y*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1)+2.692307692*E*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (2, 26) = .8569881549*E*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (2, 27) = 5.769230769*E*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (2, 28) = -115.3846154*E*Y*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1)-2.692307692*E*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (2, 29) = -.8569881549*E*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (2, 30) = -5.769230769*E*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (2, 31) = 115.3846154*E*Y*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1)+2.692307692*E*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (2, 32) = .8569881549*E*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (2, 33) = 5.769230769*E*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (2, 34) = -115.3846154*E*Y*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1)-2.692307692*E*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (2, 35) = -.8569881549*E*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (2, 36) = -5.769230769*E*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (2, 37) = -115.3846154*E*Y*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1)-2.692307692*E*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (2, 38) = -.8569881549*E*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (2, 39) = -5.769230769*E*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (2, 40) = 115.3846154*E*Y*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1)+2.692307692*E*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (2, 41) = .8569881549*E*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (2, 42) = 5.769230769*E*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (2, 43) = -115.3846154*E*Y*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1)-2.692307692*E*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (2, 44) = -.8569881549*E*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (2, 45) = -5.769230769*E*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (2, 46) = 115.3846154*E*Y*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1)+2.692307692*E*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (2, 47) = .8569881549*E*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (2, 48) = 5.769230769*E*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (3, 1) = -115.3846154*E*Y*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1)-1.153846154*E*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (3, 2) = -.3672806379*E*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (3, 3) = -13.46153846*E*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (3, 4) = 115.3846154*E*Y*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1)+1.153846154*E*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (3, 5) = .3672806379*E*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (3, 6) = 13.46153846*E*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (3, 7) = -115.3846154*E*Y*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1)-1.153846154*E*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (3, 8) = -.3672806379*E*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (3, 9) = -13.46153846*E*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (3, 10) = 115.3846154*E*Y*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1)+1.153846154*E*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (3, 11) = .3672806379*E*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (3, 12) = 13.46153846*E*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (3, 13) = 115.3846154*E*Y*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1)+1.153846154*E*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (3, 14) = .3672806379*E*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (3, 15) = 13.46153846*E*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (3, 16) = -115.3846154*E*Y*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1)-1.153846154*E*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (3, 17) = -.3672806379*E*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (3, 18) = -13.46153846*E*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (3, 19) = 115.3846154*E*Y*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1)+1.153846154*E*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (3, 20) = .3672806379*E*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (3, 21) = 13.46153846*E*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (3, 22) = -115.3846154*E*Y*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1)-1.153846154*E*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (3, 23) = -.3672806379*E*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (3, 24) = -13.46153846*E*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (3, 25) = 115.3846154*E*Y*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1)+1.153846154*E*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (3, 26) = .3672806379*E*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (3, 27) = 13.46153846*E*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (3, 28) = -115.3846154*E*Y*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1)-1.153846154*E*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (3, 29) = -.3672806379*E*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (3, 30) = -13.46153846*E*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (3, 31) = 115.3846154*E*Y*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1)+1.153846154*E*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (3, 32) = .3672806379*E*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (3, 33) = 13.46153846*E*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (3, 34) = -115.3846154*E*Y*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1)-1.153846154*E*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (3, 35) = -.3672806379*E*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (3, 36) = -13.46153846*E*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (3, 37) = -115.3846154*E*Y*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1)-1.153846154*E*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (3, 38) = -.3672806379*E*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (3, 39) = -13.46153846*E*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (3, 40) = 115.3846154*E*Y*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1)+1.153846154*E*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (3, 41) = .3672806379*E*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (3, 42) = 13.46153846*E*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (3, 43) = -115.3846154*E*Y*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1)-1.153846154*E*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (3, 44) = -.3672806379*E*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (3, 45) = -13.46153846*E*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (3, 46) = 115.3846154*E*Y*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1)+1.153846154*E*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (3, 47) = .3672806379*E*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (3, 48) = 13.46153846*E*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (4, 1) = -.2448537586*E*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (4, 2) = -.3846153846*E*(0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (4, 3) = 0., (4, 4) = .2448537586*E*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (4, 5) = .3846153846*E*(0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (4, 6) = 0., (4, 7) = -.2448537586*E*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (4, 8) = -.3846153846*E*(0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (4, 9) = 0., (4, 10) = .2448537586*E*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (4, 11) = .3846153846*E*(0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (4, 12) = 0., (4, 13) = .2448537586*E*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (4, 14) = .3846153846*E*(0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (4, 15) = 0., (4, 16) = -.2448537586*E*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (4, 17) = -.3846153846*E*(0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (4, 18) = 0., (4, 19) = .2448537586*E*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (4, 20) = .3846153846*E*(0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (4, 21) = 0., (4, 22) = -.2448537586*E*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (4, 23) = -.3846153846*E*(0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (4, 24) = 0., (4, 25) = .2448537586*E*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (4, 26) = .3846153846*E*(0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (4, 27) = 0., (4, 28) = -.2448537586*E*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (4, 29) = -.3846153846*E*(0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (4, 30) = 0., (4, 31) = .2448537586*E*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (4, 32) = .3846153846*E*(0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (4, 33) = 0., (4, 34) = -.2448537586*E*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (4, 35) = -.3846153846*E*(0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (4, 36) = 0., (4, 37) = -.2448537586*E*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (4, 38) = -.3846153846*E*(0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (4, 39) = 0., (4, 40) = .2448537586*E*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (4, 41) = .3846153846*E*(0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (4, 42) = 0., (4, 43) = -.2448537586*E*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (4, 44) = -.3846153846*E*(0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (4, 45) = 0., (4, 46) = .2448537586*E*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (4, 47) = .3846153846*E*(0.2e3*Y-2/(0.1e-1*y+.11))*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (4, 48) = 0., (5, 1) = -3.846153846*E*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (5, 2) = 0., (5, 3) = -76.92307692*E*Y*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (5, 4) = 3.846153846*E*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (5, 5) = 0., (5, 6) = 76.92307692*E*Y*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (5, 7) = -3.846153846*E*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (5, 8) = 0., (5, 9) = -76.92307692*E*Y*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (5, 10) = 3.846153846*E*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (5, 11) = 0., (5, 12) = 76.92307692*E*Y*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (5, 13) = 3.846153846*E*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (5, 14) = 0., (5, 15) = 76.92307692*E*Y*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (5, 16) = -3.846153846*E*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (5, 17) = 0., (5, 18) = -76.92307692*E*Y*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (5, 19) = 3.846153846*E*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (5, 20) = 0., (5, 21) = 76.92307692*E*Y*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (5, 22) = -3.846153846*E*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (5, 23) = 0., (5, 24) = -76.92307692*E*Y*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (5, 25) = 3.846153846*E*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (5, 26) = 0., (5, 27) = 76.92307692*E*Y*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (5, 28) = -3.846153846*E*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (5, 29) = 0., (5, 30) = -76.92307692*E*Y*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (5, 31) = 3.846153846*E*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (5, 32) = 0., (5, 33) = 76.92307692*E*Y*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (5, 34) = -3.846153846*E*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (5, 35) = 0., (5, 36) = -76.92307692*E*Y*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (5, 37) = -3.846153846*E*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (5, 38) = 0., (5, 39) = -76.92307692*E*Y*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (5, 40) = 3.846153846*E*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (5, 41) = 0., (5, 42) = 76.92307692*E*Y*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (5, 43) = -3.846153846*E*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (5, 44) = 0., (5, 45) = -76.92307692*E*Y*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (5, 46) = 3.846153846*E*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (5, 47) = 0., (5, 48) = 76.92307692*E*Y*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (6, 1) = 0., (6, 2) = -3.846153846*E*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (6, 3) = -.2448537586*E*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (6, 4) = 0., (6, 5) = 3.846153846*E*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (6, 6) = .2448537586*E*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (6, 7) = 0., (6, 8) = -3.846153846*E*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (6, 9) = -.2448537586*E*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (6, 10) = 0., (6, 11) = 3.846153846*E*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (6, 12) = .2448537586*E*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (6, 13) = 0., (6, 14) = 3.846153846*E*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1), (6, 15) = .2448537586*E*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (6, 16) = 0., (6, 17) = -3.846153846*E*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1), (6, 18) = -.2448537586*E*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (6, 19) = 0., (6, 20) = 3.846153846*E*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1), (6, 21) = .2448537586*E*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y+1)/(0.1e-1*y+.11), (6, 22) = 0., (6, 23) = -3.846153846*E*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1), (6, 24) = -.2448537586*E*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y+1)/(0.1e-1*y+.11), (6, 25) = 0., (6, 26) = 3.846153846*E*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (6, 27) = .2448537586*E*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (6, 28) = 0., (6, 29) = -3.846153846*E*X*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (6, 30) = -.2448537586*E*Z*((1/8)*cos(pi*z)-(1/8)*cos(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (6, 31) = 0., (6, 32) = 3.846153846*E*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (6, 33) = .2448537586*E*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (6, 34) = 0., (6, 35) = -3.846153846*E*X*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (6, 36) = -.2448537586*E*Z*((1/8)*sin(pi*z)-(1/8)*sin(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (6, 37) = 0., (6, 38) = -3.846153846*E*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1), (6, 39) = -.2448537586*E*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (6, 40) = 0., (6, 41) = 3.846153846*E*X*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1), (6, 42) = .2448537586*E*Z*((1/8)*cos(pi*z)+(1/8)*cos(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11), (6, 43) = 0., (6, 44) = -3.846153846*E*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1), (6, 45) = -.2448537586*E*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x+1)*(y-1)/(0.1e-1*y+.11), (6, 46) = 0., (6, 47) = 3.846153846*E*X*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1), (6, 48) = .2448537586*E*Z*((1/8)*sin(pi*z)+(1/8)*sin(pi*z)^2)*(x-1)*(y-1)/(0.1e-1*y+.11)})), true)

(9)

int(int(int(S, z = -1 .. 1), y = -1 .. 1), x = -1 .. 1)

Warning,  computation interrupted

 

NULL

 

Download maple2.mw

3 4 5 6 7 8 9 Last Page 5 of 46