Dp11	1+u1x
Dp21	u2x
Dp31	u3x
Dp12	u1y
Dp22	1+u2y
Dp32	u3y
Dp13	u1z
Dp23	u2z
Dp33	1+u3z
f11	Dp11/beta^(1/3)
f21	Dp21/beta^(1/3)
f31	Dp31/beta^(1/3)
f12	Dp12/beta^(1/3)
f22	Dp22/beta^(1/3)
f32	Dp32/beta^(1/3)
f13	Dp13/beta^(1/3)
f23	Dp23/beta^(1/3)
f33	Dp33/beta^(1/3)
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*beta^(2/3)
cfDp21	cfF21*beta^(2/3)
cfDp31	cfF31*beta^(2/3)
cfDp12	cfF12*beta^(2/3)
cfDp22	cfF22*beta^(2/3)
cfDp32	cfF32*beta^(2/3)
cfDp13	cfF13*beta^(2/3)
cfDp23	cfF23*beta^(2/3)
cfDp33	cfF33*beta^(2/3)
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/beta^(1/3)-(1/3)*(1+u1x)/beta^(4/3)*betat
f21t	u2xt/beta^(1/3)-(1/3)*u2x/beta^(4/3)*betat
f31t	u3xt/beta^(1/3)-(1/3)*u3x/beta^(4/3)*betat
f12t	u1yt/beta^(1/3)-(1/3)*u1y/beta^(4/3)*betat
f22t	u2yt/beta^(1/3)-(1/3)*(1+u2y)/beta^(4/3)*betat
f32t	u3yt/beta^(1/3)-(1/3)*u3y/beta^(4/3)*betat
f13t	u1zt/beta^(1/3)-(1/3)*u1z/beta^(4/3)*betat
f23t	u2zt/beta^(1/3)-(1/3)*u2z/beta^(4/3)*betat
f33t	u3zt/beta^(1/3)-(1/3)*(1+u3z)/beta^(4/3)*betat
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
So11	(Se11+Sdiss11)*beta^(2/3)
So21	(Se21+Sdiss21)*beta^(2/3)
So31	(Se31+Sdiss31)*beta^(2/3)
So12	(Se12+Sdiss12)*beta^(2/3)
So22	(Se22+Sdiss22)*beta^(2/3)
So32	(Se32+Sdiss32)*beta^(2/3)
So13	(Se13+Sdiss13)*beta^(2/3)
So23	(Se23+Sdiss23)*beta^(2/3)
So33	(Se33+Sdiss33)*beta^(2/3)
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
pEsh	phi-detF*TSph
DpEsh1	d(pEsh,x)
DpEsh2	d(pEsh,y)
DpEsh3	d(pEsh,z)
chmu	pEsh*(alpha/rho0)
Dchmu1	DpEsh1*(alpha/rho0)
Dchmu2	DpEsh2*(alpha/rho0)
Dchmu3	DpEsh3*(alpha/rho0)
ho1	-D*Dchmu1
ho2	-D*Dchmu2
ho3	-D*Dchmu3
beta	1+alpha*c
betat	alpha*ct
tchmu	test(c)*1[Pa*m^3/mol]
tDchmu1	test(cx)*1[Pa*m^3/mol]
tDchmu2	test(cy)*1[Pa*m^3/mol]
tDchmu3	test(cz)*1[Pa*m^3/mol]
tv1	test(u1)*1[1/s]
tv2	test(u2)*1[1/s]
tv3	test(u3)*1[1/s]
tDv11	test(u1x)*1[1/s]
tDv21	test(u2x)*1[1/s]
tDv31	test(u3x)*1[1/s]
tDv12	test(u1y)*1[1/s]
tDv22	test(u2y)*1[1/s]
tDv32	test(u3y)*1[1/s]
tDv13	test(u1z)*1[1/s]
tDv23	test(u2z)*1[1/s]
tDv33	test(u3z)*1[1/s]