Dp11	1+u1x
Dp21	u2x
Dp31	u3x 
Dp12	u1y
Dp22	1+u2y
Dp32	u3y 
Dp13	u1z
Dp23	u2z
Dp33	1+u3z
f11	Dp11
f21	Dp21
f31	Dp31
f12	Dp12
f22	Dp22
f32	Dp32
f13	Dp13
f23	Dp23
f33	Dp33 
B11	f11*f11+f12*f12+f13*f13
B12	f11*f21+f12*f22+f13*f23
B13	f11*f31+f12*f32+f13*f33
B21	f21*f11+f22*f12+f23*f13
B22	f21*f21+f22*f22+f23*f23
B23	f21*f31+f22*f32+f23*f33
B31	f31*f11+f32*f12+f33*f13
B32	f31*f21+f32*f22+f33*f23
B33	f31*f31+f32*f32+f33*f33
detF	-(f13*f22*f31)+f12*f23*f31+f13*f21*f32-f11*f23*f32-f12*f21*f33+f11*f22*f33
J	detF
i1	B11+B22+B33
i2	(i1^2-(B11^2+B12^2+B13^2+B21^2+B22^2+B23^2+B31^2+B32^2+B33^2))/2
cfF11	-(f23*f32)+f22*f33
cfF12	f23*f31-f21*f33
cfF13	-(f22*f31)+f21*f32
cfF21	f13*f32-f12*f33
cfF22	-(f13*f31)+f11*f33
cfF23	f12*f31-f11*f32
cfF31	-(f13*f22)+f12*f23
cfF32	f13*f21-f11*f23
cfF33	-(f12*f21)+f11*f22
cfDp11	cfF11 
cfDp21	cfF21
cfDp31	cfF31
cfDp12	cfF12
cfDp22	cfF22
cfDp32	cfF32
cfDp13	cfF13
cfDp23	cfF23
cfDp33	cfF33
finv11	cfF11/detF
finv21	cfF12/detF
finv31	cfF13/detF
finv12	cfF21/detF
finv22	cfF22/detF
finv32	cfF23/detF
finv13	cfF31/detF
finv23	cfF32/detF
finv33	cfF33/detF
f11t	u1xt
f21t	u2xt
f31t	u3xt
f12t	u1yt
f22t	u2yt
f32t	u3yt
f13t	u1zt
f23t	u2zt
f33t	u3zt
Dv11	f11t*finv11+f12t*finv21+f13t*finv31
Dv21	f21t*finv11+f22t*finv21+f23t*finv31
Dv31	f31t*finv11+f32t*finv21+f33t*finv31
Dv12	f11t*finv12+f12t*finv22+f13t*finv32
Dv22	f21t*finv12+f22t*finv22+f23t*finv32
Dv32	f31t*finv12+f32t*finv22+f33t*finv32
Dv13	f11t*finv13+f12t*finv23+f13t*finv33
Dv23	f21t*finv13+f22t*finv23+f23t*finv33
Dv33	f31t*finv13+f32t*finv23+f33t*finv33
S11	(Se11+Sdiss11)
S21	(Se21+Sdiss21)
S31	(Se31+Sdiss31)
S12	(Se12+Sdiss12)
S22	(Se22+Sdiss22)
S32	(Se32+Sdiss32)
S13	(Se13+Sdiss13)
S23	(Se23+Sdiss23)
S33	(Se33+Sdiss33)
Se11	2*kI*J^(-2/3)*(f11-(i1/3)*finv11)+2*kV*(J-1)*cfF11
Se12	2*kI*J^(-2/3)*(f12-(i1/3)*finv21)+2*kV*(J-1)*cfF12
Se13	2*kI*J^(-2/3)*(f13-(i1/3)*finv31)+2*kV*(J-1)*cfF13
Se21	2*kI*J^(-2/3)*(f21-(i1/3)*finv12)+2*kV*(J-1)*cfF21
Se22	2*kI*J^(-2/3)*(f22-(i1/3)*finv22)+2*kV*(J-1)*cfF22
Se23	2*kI*J^(-2/3)*(f23-(i1/3)*finv32)+2*kV*(J-1)*cfF23
Se31	2*kI*J^(-2/3)*(f31-(i1/3)*finv13)+2*kV*(J-1)*cfF31
Se32	2*kI*J^(-2/3)*(f32-(i1/3)*finv23)+2*kV*(J-1)*cfF32
Se33	2*kI*J^(-2/3)*(f33-(i1/3)*finv33)+2*kV*(J-1)*cfF33
Sdiss11	cfF11*Tdiss11+cfF21*Tdiss12+cfF31*Tdiss13
Sdiss12	cfF12*Tdiss11+cfF22*Tdiss12+cfF32*Tdiss13
Sdiss13	cfF13*Tdiss11+cfF23*Tdiss12+cfF33*Tdiss13
Sdiss21	cfF11*Tdiss21+cfF21*Tdiss22+cfF31*Tdiss23
Sdiss22	cfF12*Tdiss21+cfF22*Tdiss22+cfF32*Tdiss23
Sdiss23	cfF13*Tdiss21+cfF23*Tdiss22+cfF33*Tdiss23
Sdiss31	cfF11*Tdiss31+cfF21*Tdiss32+cfF31*Tdiss33
Sdiss32	cfF12*Tdiss31+cfF22*Tdiss32+cfF32*Tdiss33
Sdiss33	cfF13*Tdiss31+cfF23*Tdiss32+cfF33*Tdiss33
T11	(Te11+Tdiss11)
T21	(Te21+Tdiss21)
T31	(Te31+Tdiss31)
T12	(Te12+Tdiss12)
T22	(Te22+Tdiss22)
T32	(Te32+Tdiss32)
T13	(Te13+Tdiss13)
T23	(Te23+Tdiss23)
T33	(Te33+Tdiss33)
Te11	2*kI*J^(-5/3)*(B11-i1/3)+2*kV*(J-1)
Te12	2*kI*J^(-5/3)*B12
Te13	2*kI*J^(-5/3)*B13
Te21	2*kI*J^(-5/3)*B21
Te22	2*kI*J^(-5/3)*(B22-i1/3)+2*kV*(J-1)
Te23	2*kI*J^(-5/3)*B23
Te31	2*kI*J^(-5/3)*B31
Te32	2*kI*J^(-5/3)*B32
Te33	2*kI*J^(-5/3)*(B33-i1/3)+2*kV*(J-1) 
Tdiss11	mu*(Dv11+Dv11)
Tdiss12	mu*(Dv12+Dv21)
Tdiss13	mu*(Dv13+Dv31)
Tdiss21	mu*(Dv21+Dv12)
Tdiss22	mu*(Dv22+Dv22)
Tdiss23	mu*(Dv23+Dv32)
Tdiss31	mu*(Dv31+Dv13)
Tdiss32	mu*(Dv32+Dv23)
Tdiss33	mu*(Dv33+Dv33)
TSph	(T11+T22+T33)/3
phi	kI*(J^(-2/3)*i1-3)+kV*(J-1)^2