In[5306]:= (*Create an s-type function*)
(*First primitive of the s-type function*)
x0s:=0.0
y0s:=0.0
z0s:=0.0
alphas1:=1.0
s1[x_,y_,z_] := Exp[-alphas1*((x-x0s)^2+(y-y0s)^2+(z-z0s)^2)]*(2.0*alphas1/Pi)^0.75
(*Second primitive of the s-type function*)
alphas2:=2.0
s2[x_,y_,z_] := Exp[-alphas2*((x-x0s)^2+(y-y0s)^2+(z-z0s)^2)]*(2.0*alphas2/Pi)^0.75
(*An s-type function as a contraction from the above primitives*)
s[x_, y_, z_] = s1[x, y, z] * 1.0 + s2[x, y, z] * 0.1;
(*Form derivatives*)
sX[x_,y_,z_] = D[s[x,y,z],x];
sY[x_,y_,z_] = D[s[x,y,z],y];
sZ[x_, y_, z_] = D[s[x,y,z],z];
sXX[x_,y_,z_] = D[sX[x,y,z],x];
sXY[x_,y_,z_] = D[sX[x,y,z],y];
sXZ[x_,y_,z_] = D[sX[x,y,z],z];
sYY[x_,y_,z_] = D[sY[x,y,z],y];
sYZ[x_,y_,z_] = D[sY[x,y,z],z];
sZZ[x_,y_,z_] = D[sZ[x,y,z],z];

(*Create a set of p-type function*)
(*First primitive of the px-type function*)
lpx:=1
mpx:=0
npx:=0
x0p:=1.0
y0p:=0.0
z0p:=0.0
alphap1:=1.0
p1x[x_,y_,z_] := (x-x0p)^lpx*Exp[-alphap1*((x-x0p)^2+(y-y0p)^2+(z-z0p)^2)]*(2.0*alphap1/Pi)^0.75*(4*alphap1)^(0.5*(lpx+mpx+npx))
(*First primitive of the py-type function*)
lpy:=0
mpy:=1
npy:=0
p1y[x_,y_,z_] := (y-y0p)^mpy*Exp[-alphap1*((x-x0p)^2+(y-y0p)^2+(z-z0p)^2)]*(2.0*alphap1/Pi)^0.75*(4*alphap1)^(0.5*(lpy+mpy+npy))
(*First primitive of the pz-type function*)
lpz:=0
mpz:=0
npz:=1
p1z[x_,y_,z_] := (z-z0p)^npz*Exp[-alphap1*((x-x0p)^2+(y-y0p)^2+(z-z0p)^2)]*(2.0*alphap1/Pi)^0.75*(4*alphap1)^(0.5*(lpz+mpz+npz))
(*Second primitive of the px-type function*)
alphap2:=2.0
p2x[x_,y_,z_] := (x-x0p)^lpx*Exp[-alphap2*((x-x0p)^2+(y-y0p)^2+(z-z0p)^2)]*(2.0*alphap2/Pi)^0.75*(4*alphap2)^(0.5*(lpx+mpx+npx))
(*Second primitive of the py-type function*)
p2y[x_,y_,z_] := (y-y0p)^mpy*Exp[-alphap2*((x-x0p)^2+(y-y0p)^2+(z-z0p)^2)]*(2.0*alphap2/Pi)^0.75*(4*alphap2)^(0.5*(lpy+mpy+npy))
(*Second primitive of the pz-type function*)
p2z[x_,y_,z_] := (z-z0p)^npz*Exp[-alphap2*((x-x0p)^2+(y-y0p)^2+(z-z0p)^2)]*(2.0*alphap2/Pi)^0.75*(4*alphap2)^(0.5*(lpz+mpz+npz))
(*A set of p-type functions as a contraction from the above primitives*)
px[x_, y_, z_] = p1x[x, y, z] * 1.0 + p2x[x, y, z] * 0.1;
py[x_, y_, z_] = p1y[x, y, z] * 1.0 + p2y[x, y, z] * 0.1;
pz[x_, y_, z_] = p1z[x, y, z] * 1.0 + p2z[x, y, z] * 0.1;
(*Form derivatives*)
pxX[x_,y_,z_] = D[px[x,y,z],x];
pxY[x_,y_,z_] = D[px[x,y,z],y];
pxZ[x_, y_, z_] = D[px[x,y,z],z];
pxXX[x_,y_,z_] = D[pxX[x,y,z],x];
pxXY[x_,y_,z_] = D[pxX[x,y,z],y];
pxXZ[x_,y_,z_] = D[pxX[x,y,z],z];
pxYY[x_,y_,z_] = D[pxY[x,y,z],y];
pxYZ[x_,y_,z_] = D[pxY[x,y,z],z];
pxZZ[x_,y_,z_] = D[pxZ[x,y,z],z];
pyX[x_,y_,z_] = D[py[x,y,z],x];
pyY[x_,y_,z_] = D[py[x,y,z],y];
pyZ[x_, y_, z_] = D[py[x,y,z],z];
pyXX[x_,y_,z_] = D[pyX[x,y,z],x];
pyXY[x_,y_,z_] = D[pyX[x,y,z],y];
pyXZ[x_,y_,z_] = D[pyX[x,y,z],z];
pyYY[x_,y_,z_] = D[pyY[x,y,z],y];
pyYZ[x_,y_,z_] = D[pyY[x,y,z],z];
pyZZ[x_,y_,z_] = D[pyZ[x,y,z],z];
pzX[x_,y_,z_] = D[pz[x,y,z],x];
pzY[x_,y_,z_] = D[pz[x,y,z],y];
pzZ[x_, y_, z_] = D[pz[x,y,z],z];
pzXX[x_,y_,z_] = D[pzX[x,y,z],x];
pzXY[x_,y_,z_] = D[pzX[x,y,z],y];
pzXZ[x_,y_,z_] = D[pzX[x,y,z],z];
pzYY[x_,y_,z_] = D[pzY[x,y,z],y];
pzYZ[x_,y_,z_] = D[pzY[x,y,z],z];
pzZZ[x_,y_,z_] = D[pzZ[x,y,z],z];
(*Create a set of d-type function*)
(*First primitive of the dyy-type function*)
ldxx:=2
mdxx:=0
ndxx:=0
x0d:=-1.0
y0d:=2.5
z0d:=1.0
alphad1:=0.5
d1xx[x_,y_,z_] := (x-x0d)^ldxx*Exp[-alphad1*((x-x0d)^2+(y-y0d)^2+(z-z0d)^2)]*(2.0*alphad1/Pi)^0.75*(4*alphad1)^(0.5*(ldxx+mdxx+ndxx))*(1/Sqrt[(2*ldxx-1)!!])
(*First primitive of the dxy-type function*)
ldxy:=1
mdxy:=1
ndxy:=0
d1xy[x_,y_,z_] := (x-x0d)^ldxy*(y-y0d)^mdxy*(z-z0d)^ndxy*Exp[-alphad1*((x-x0d)^2+(y-y0d)^2+(z-z0d)^2)]*(2.0*alphad1/Pi)^0.75*(4*alphad1)^(0.5*(ldxy+mdxy+ndxy))
(*First primitive of the dxz-type function*)
ldxz:=1
mdxz:=0
ndxz:=1
d1xz[x_,y_,z_] := (x-x0d)^ldxz*(y-y0d)^mdxz*(z-z0d)^ndxz*Exp[-alphad1*((x-x0d)^2+(y-y0d)^2+(z-z0d)^2)]*(2.0*alphad1/Pi)^0.75*(4*alphad1)^(0.5*(ldxz+mdxz+ndxz))
(*First primitive of the dyy-type function*)
ldyy:=0
mdyy:=2
ndyy:=0
d1yy[x_,y_,z_] := (x-x0d)^ldyy*(y-y0d)^mdyy*(z-z0d)^ndyy*Exp[-alphad1*((x-x0d)^2+(y-y0d)^2+(z-z0d)^2)]*(2.0*alphad1/Pi)^0.75*(4*alphad1)^(0.5*(ldyy+mdyy+ndyy))*(1/Sqrt[(2*ldyy-1)!!])
(*First primitive of the dyz-type function*)
ldyz:=0
mdyz:=1
ndyz:=1
d1yz[x_,y_,z_] := (x-x0d)^ldyz*(y-y0d)^mdyz*(z-z0d)^ndyz*Exp[-alphad1*((x-x0d)^2+(y-y0d)^2+(z-z0d)^2)]*(2.0*alphad1/Pi)^0.75*(4*alphad1)^(0.5*(ldyz+mdyz+ndyz))
(*First primitive of the dzz-type function*)
ldzz:=1
mdzz:=1
ndzz:=0
d1zz[x_,y_,z_] := (x-x0d)^ldzz*(y-y0d)^mdzz*(z-z0d)^ndzz*Exp[-alphad1*((x-x0d)^2+(y-y0d)^2+(z-z0d)^2)]*(2.0*alphad1/Pi)^0.75*(4*alphad1)^(0.5*(ldzz+mdzz+ndzz))*(1/Sqrt[(2*ldzz-1)!!])
(*Second primitive of the dxx-type function*)
alphad2:=1.2
d2xx[x_,y_,z_] := (x-x0d)^ldxx*Exp[-alphad2*((x-x0d)^2+(y-y0d)^2+(z-z0d)^2)]*(2.0*alphad2/Pi)^0.75*(4*alphad2)^(0.5*(ldxx+mdxx+ndxx))*(1/Sqrt[(2*ldxx-1)!!])
(*Second primitive of the dxy-type function*)
d2xy[x_,y_,z_] := (x-x0d)^ldxy*(y-y0d)^mdxy*(z-z0d)^ndxy*Exp[-alphad2*((x-x0d)^2+(y-y0d)^2+(z-z0d)^2)]*(2.0*alphad2/Pi)^0.75*(4*alphad2)^(0.5*(ldxy+mdxy+ndxy))
(*Second primitive of the dxz-type function*)
d2xz[x_,y_,z_] := (x-x0d)^ldxz*(y-y0d)^mdxz*(z-z0d)^ndxz*Exp[-alphad2*((x-x0d)^2+(y-y0d)^2+(z-z0d)^2)]*(2.0*alphad2/Pi)^0.75*(4*alphad2)^(0.5*(ldxz+mdxz+ndxz))
(*Second primitive of the dyy-type function*)
d2yy[x_,y_,z_] := (x-x0d)^ldyy*(y-y0d)^mdyy*(z-z0d)^ndyy*Exp[-alphad2*((x-x0d)^2+(y-y0d)^2+(z-z0d)^2)]*(2.0*alphad2/Pi)^0.75*(4*alphad2)^(0.5*(ldyy+mdyy+ndyy))*(1/Sqrt[(2*ldyy-1)!!])
(*Second primitive of the dyz-type function*)
d2yz[x_,y_,z_] := (x-x0d)^ldyz*(y-y0d)^mdyz*(z-z0d)^ndyz*Exp[-alphad2*((x-x0d)^2+(y-y0d)^2+(z-z0d)^2)]*(2.0*alphad2/Pi)^0.75*(4*alphad2)^(0.5*(ldyz+mdyz+ndyz))
(*Second primitive of the dzz-type function*)
d2zz[x_,y_,z_] := (x-x0d)^ldzz*(y-y0d)^mdzz*(z-z0d)^ndzz*Exp[-alphad2*((x-x0d)^2+(y-y0d)^2+(z-z0d)^2)]*(2.0*alphad2/Pi)^0.75*(4*alphad2)^(0.5*(ldzz+mdzz+ndzz))*(1/Sqrt[(2*ldzz-1)!!])
(*A set of d-type functions as a contraction from the above primitives*)
dxx[x_, y_, z_] = d1xx[x, y, z] * 0.8 + d2xx[x, y, z] * 0.4;
dxy[x_, y_, z_] = d1xy[x, y, z] * 0.8 + d2xy[x, y, z] * 0.4;
dxz[x_, y_, z_] = d1xz[x, y, z] * 0.8 + d2xz[x, y, z] * 0.4;
dyy[x_, y_, z_] = d1yy[x, y, z] * 0.8 + d2yy[x, y, z] * 0.4;
dyz[x_, y_, z_] = d1yz[x, y, z] * 0.8 + d2yz[x, y, z] * 0.4;
dzz[x_, y_, z_] = d1zz[x, y, z] * 0.8 + d2zz[x, y, z] * 0.4;
(*Form derivatives*)
dxxX[x_,y_,z_] = D[dxx[x,y,z],x];
dxxY[x_,y_,z_] = D[dxx[x,y,z],y];
dxxZ[x_,y_,z_] = D[dxx[x,y,z],z];
dxxXX[x_,y_,z_] = D[dxxX[x,y,z],x];
dxxXY[x_,y_,z_] = D[dxxX[x,y,z],y];
dxxXZ[x_,y_,z_] = D[dxxX[x,y,z],y];
dxxYY[x_,y_,z_] = D[dxxY[x,y,z],y];
dxxYZ[x_,y_,z_] = D[dxxY[x,y,z],z];
dxxZZ[x_,y_,z_] = D[dxxZ[x,y,z],z];
dxyX[x_,y_,z_] = D[dxy[x,y,z],x];
dxyY[x_,y_,z_] = D[dxy[x,y,z],y];
dxyZ[x_,y_,z_] = D[dxy[x,y,z],z];
dxyXX[x_,y_,z_] = D[dxyX[x,y,z],x];
dxyXY[x_,y_,z_] = D[dxyX[x,y,z],y];
dxyXZ[x_,y_,z_] = D[dxyX[x,y,z],y];
dxyYY[x_,y_,z_] = D[dxyY[x,y,z],y];
dxyYZ[x_,y_,z_] = D[dxyY[x,y,z],z];
dxyZZ[x_,y_,z_] = D[dxyZ[x,y,z],z];
dxzX[x_,y_,z_] = D[dxz[x,y,z],x];
dxzY[x_,y_,z_] = D[dxz[x,y,z],y];
dxzZ[x_,y_,z_] = D[dxz[x,y,z],z];
dxzXX[x_,y_,z_] = D[dxzX[x,y,z],x];
dxzXY[x_,y_,z_] = D[dxzX[x,y,z],y];
dxzXZ[x_,y_,z_] = D[dxzX[x,y,z],y];
dxzYY[x_,y_,z_] = D[dxzY[x,y,z],y];
dxzYZ[x_,y_,z_] = D[dxzY[x,y,z],z];
dxzZZ[x_,y_,z_] = D[dxzZ[x,y,z],z];
dyyX[x_,y_,z_] = D[dyy[x,y,z],x];
dyyY[x_,y_,z_] = D[dyy[x,y,z],y];
dyyZ[x_,y_,z_] = D[dyy[x,y,z],z];
dyyXX[x_,y_,z_] = D[dyyX[x,y,z],x];
dyyXY[x_,y_,z_] = D[dyyX[x,y,z],y];
dyyXZ[x_,y_,z_] = D[dyyX[x,y,z],y];
dyyYY[x_,y_,z_] = D[dyyY[x,y,z],y];
dyyYZ[x_,y_,z_] = D[dyyY[x,y,z],z];
dyyZZ[x_,y_,z_] = D[dyyZ[x,y,z],z];
dyzX[x_,y_,z_] = D[dyz[x,y,z],x];
dyzY[x_,y_,z_] = D[dyz[x,y,z],y];
dyzZ[x_,y_,z_] = D[dyz[x,y,z],z];
dyzXX[x_,y_,z_] = D[dyzX[x,y,z],x];
dyzXY[x_,y_,z_] = D[dyzX[x,y,z],y];
dyzXZ[x_,y_,z_] = D[dyzX[x,y,z],y];
dyzYY[x_,y_,z_] = D[dyzY[x,y,z],y];
dyzYZ[x_,y_,z_] = D[dyzY[x,y,z],z];
dyzZZ[x_,y_,z_] = D[dyzZ[x,y,z],z];
dzzX[x_,y_,z_] = D[dzz[x,y,z],x];
dzzY[x_,y_,z_] = D[dzz[x,y,z],y];
dzzZ[x_,y_,z_] = D[dzz[x,y,z],z];
dzzXX[x_,y_,z_] = D[dzzX[x,y,z],x];
dzzXY[x_,y_,z_] = D[dzzX[x,y,z],y];
dzzXZ[x_,y_,z_] = D[dzzX[x,y,z],y];
dzzYY[x_,y_,z_] = D[dzzY[x,y,z],y];
dzzYZ[x_,y_,z_] = D[dzzY[x,y,z],z];
dzzZZ[x_,y_,z_] = D[dzzZ[x,y,z],z];

(*Points to evaluate at*)
P0x:=0.0
P0y:=0.0
P0z:=0.0
P1x:=1.0
P1y:=0.0
P1z:=0.0
P2x:=-2.0
P2y:=1.0
P2z:=1.0
P3x:=3.0
P3y:=0.5
P3z:=-1.0
(*Produce results*)
Print["S-Function"]
SetPrecision[ s[P0x, P0y, P0z],12]
SetPrecision[ s[P1x, P1y, P1z],12]
SetPrecision[ s[P2x, P2y, P2z],12]
SetPrecision[ s[P3x, P3y, P3z],12]
Print["S-Function - d/dx"]
SetPrecision[ sX[P0x, P0y, P0z],12]
SetPrecision[ sX[P1x, P1y, P1z],12]
SetPrecision[ sX[P2x, P2y, P2z],12]
SetPrecision[ sX[P3x, P3y, P3z],12]
Print["S-Function - d/dy"]
SetPrecision[ sY[P0x, P0y, P0z],12]
SetPrecision[ sY[P1x, P1y, P1z],12]
SetPrecision[ sY[P2x, P2y, P2z],12]
SetPrecision[ sY[P3x, P3y, P3z],12]
Print["S-Function - d/dz"]
SetPrecision[ sZ[P0x, P0y, P0z],12]
SetPrecision[ sZ[P1x, P1y, P1z],12]
SetPrecision[ sZ[P2x, P2y, P2z],12]
SetPrecision[ sZ[P3x, P3y, P3z],12]
Print["S-Function - d^2/dx^2"]
SetPrecision[ sXX[P0x, P0y, P0z],12]
SetPrecision[ sXX[P1x, P1y, P1z],12]
SetPrecision[ sXX[P2x, P2y, P2z],12]
SetPrecision[ sXX[P3x, P3y, P3z],12]
Print["S-Function - d^2/dxdy"]
SetPrecision[ sXY[P0x, P0y, P0z],12]
SetPrecision[ sXY[P1x, P1y, P1z],12]
SetPrecision[ sXY[P2x, P2y, P2z],12]
SetPrecision[ sXY[P3x, P3y, P3z],12]
Print["S-Function - d^2/dxdz"]
SetPrecision[ sXZ[P0x, P0y, P0z],12]
SetPrecision[ sXZ[P1x, P1y, P1z],12]
SetPrecision[ sXZ[P2x, P2y, P2z],12]
SetPrecision[ sXZ[P3x, P3y, P3z],12]
Print["S-Function - d^2/dy^2"]
SetPrecision[ sYY[P0x, P0y, P0z],12]
SetPrecision[ sYY[P1x, P1y, P1z],12]
SetPrecision[ sYY[P2x, P2y, P2z],12]
SetPrecision[ sYY[P3x, P3y, P3z],12]
Print["S-Function - d^2/dydz"]
SetPrecision[ sYZ[P0x, P0y, P0z],12]
SetPrecision[ sYZ[P1x, P1y, P1z],12]
SetPrecision[ sYZ[P2x, P2y, P2z],12]
SetPrecision[ sYZ[P3x, P3y, P3z],12]
Print["S-Function - d^2/dz^2"]
SetPrecision[ sZZ[P0x, P0y, P0z],12]
SetPrecision[ sZZ[P1x, P1y, P1z],12]
SetPrecision[ sZZ[P2x, P2y, P2z],12]
SetPrecision[ sZZ[P3x, P3y, P3z],12]
Print["Px-Function"]
SetPrecision[ px[P0x, P0y, P0z],12]
SetPrecision[ px[P1x, P1y, P1z],12]
SetPrecision[ px[P2x, P2y, P2z],12]
SetPrecision[ px[P3x, P3y, P3z],12]
Print["Px-Function - d/dx"]
SetPrecision[ pxX[P0x, P0y, P0z],12]
SetPrecision[ pxX[P1x, P1y, P1z],12]
SetPrecision[ pxX[P2x, P2y, P2z],12]
SetPrecision[ pxX[P3x, P3y, P3z],12]
Print["Px-Function - d/dy"]
SetPrecision[ pxY[P0x, P0y, P0z],12]
SetPrecision[ pxY[P1x, P1y, P1z],12]
SetPrecision[ pxY[P2x, P2y, P2z],12]
SetPrecision[ pxY[P3x, P3y, P3z],12]
Print["Px-Function - d/dz"]
SetPrecision[ pxZ[P0x, P0y, P0z],12]
SetPrecision[ pxZ[P1x, P1y, P1z],12]
SetPrecision[ pxZ[P2x, P2y, P2z],12]
SetPrecision[ pxZ[P3x, P3y, P3z],12]
Print["Px-Function - d^2/dx^2"]
SetPrecision[ pxXX[P0x, P0y, P0z],12]
SetPrecision[ pxXX[P1x, P1y, P1z],12]
SetPrecision[ pxXX[P2x, P2y, P2z],12]
SetPrecision[ pxXX[P3x, P3y, P3z],12]
Print["Px-Function - d^2/dxdy"]
SetPrecision[ pxXY[P0x, P0y, P0z],12]
SetPrecision[ pxXY[P1x, P1y, P1z],12]
SetPrecision[ pxXY[P2x, P2y, P2z],12]
SetPrecision[ pxXY[P3x, P3y, P3z],12]
Print["Px-Function - d^2/dxdz"]
SetPrecision[ pxXZ[P0x, P0y, P0z],12]
SetPrecision[ pxXZ[P1x, P1y, P1z],12]
SetPrecision[ pxXZ[P2x, P2y, P2z],12]
SetPrecision[ pxXZ[P3x, P3y, P3z],12]
Print["Px-Function - d^2/dy^2"]
SetPrecision[ pxYY[P0x, P0y, P0z],12]
SetPrecision[ pxYY[P1x, P1y, P1z],12]
SetPrecision[ pxYY[P2x, P2y, P2z],12]
SetPrecision[ pxYY[P3x, P3y, P3z],12]
Print["Px-Function - d^2/dydz"]
SetPrecision[ pxYZ[P0x, P0y, P0z],12]
SetPrecision[ pxYZ[P1x, P1y, P1z],12]
SetPrecision[ pxYZ[P2x, P2y, P2z],12]
SetPrecision[ pxYZ[P3x, P3y, P3z],12]
Print["Px-Function - d^2/dz^2"]
SetPrecision[ pxZZ[P0x, P0y, P0z],12]
SetPrecision[ pxZZ[P1x, P1y, P1z],12]
SetPrecision[ pxZZ[P2x, P2y, P2z],12]
SetPrecision[ pxZZ[P3x, P3y, P3z],12]
Print["Py-Function"]
SetPrecision[ py[P0x, P0y, P0z],12]
SetPrecision[ py[P1x, P1y, P1z],12]
SetPrecision[ py[P2x, P2y, P2z],12]
SetPrecision[ py[P3x, P3y, P3z],12]
Print["Py-Function - d/dx"]
SetPrecision[ pyX[P0x, P0y, P0z],12]
SetPrecision[ pyX[P1x, P1y, P1z],12]
SetPrecision[ pyX[P2x, P2y, P2z],12]
SetPrecision[ pyX[P3x, P3y, P3z],12]
Print["Py-Function - d/dy"]
SetPrecision[ pyY[P0x, P0y, P0z],12]
SetPrecision[ pyY[P1x, P1y, P1z],12]
SetPrecision[ pyY[P2x, P2y, P2z],12]
SetPrecision[ pyY[P3x, P3y, P3z],12]
Print["Py-Function - d/dz"]
SetPrecision[ pyZ[P0x, P0y, P0z],12]
SetPrecision[ pyZ[P1x, P1y, P1z],12]
SetPrecision[ pyZ[P2x, P2y, P2z],12]
SetPrecision[ pyZ[P3x, P3y, P3z],12]
Print["Py-Function - d^2/dx^2"]
SetPrecision[ pyXX[P0x, P0y, P0z],12]
SetPrecision[ pyXX[P1x, P1y, P1z],12]
SetPrecision[ pyXX[P2x, P2y, P2z],12]
SetPrecision[ pyXX[P3x, P3y, P3z],12]
Print["Py-Function - d^2/dxdy"]
SetPrecision[ pyXY[P0x, P0y, P0z],12]
SetPrecision[ pyXY[P1x, P1y, P1z],12]
SetPrecision[ pyXY[P2x, P2y, P2z],12]
SetPrecision[ pyXY[P3x, P3y, P3z],12]
Print["Py-Function - d^2/dxdz"]
SetPrecision[ pyXZ[P0x, P0y, P0z],12]
SetPrecision[ pyXZ[P1x, P1y, P1z],12]
SetPrecision[ pyXZ[P2x, P2y, P2z],12]
SetPrecision[ pyXZ[P3x, P3y, P3z],12]
Print["Py-Function - d^2/dy^2"]
SetPrecision[ pyYY[P0x, P0y, P0z],12]
SetPrecision[ pyYY[P1x, P1y, P1z],12]
SetPrecision[ pyYY[P2x, P2y, P2z],12]
SetPrecision[ pyYY[P3x, P3y, P3z],12]
Print["Py-Function - d^2/dydz"]
SetPrecision[ pyYZ[P0x, P0y, P0z],12]
SetPrecision[ pyYZ[P1x, P1y, P1z],12]
SetPrecision[ pyYZ[P2x, P2y, P2z],12]
SetPrecision[ pyYZ[P3x, P3y, P3z],12]
Print["Py-Function - d^2/dz^2"]
SetPrecision[ pyZZ[P0x, P0y, P0z],12]
SetPrecision[ pyZZ[P1x, P1y, P1z],12]
SetPrecision[ pyZZ[P2x, P2y, P2z],12]
SetPrecision[ pyZZ[P3x, P3y, P3z],12]
Print["Pz-Function"]
SetPrecision[ pz[P0x, P0y, P0z],12]
SetPrecision[ pz[P1x, P1y, P1z],12]
SetPrecision[ pz[P2x, P2y, P2z],12]
SetPrecision[ pz[P3x, P3y, P3z],12]
Print["Pz-Function - d/dx"]
SetPrecision[ pzX[P0x, P0y, P0z],12]
SetPrecision[ pzX[P1x, P1y, P1z],12]
SetPrecision[ pzX[P2x, P2y, P2z],12]
SetPrecision[ pzX[P3x, P3y, P3z],12]
Print["Pz-Function - d/dy"]
SetPrecision[ pzY[P0x, P0y, P0z],12]
SetPrecision[ pzY[P1x, P1y, P1z],12]
SetPrecision[ pzY[P2x, P2y, P2z],12]
SetPrecision[ pzY[P3x, P3y, P3z],12]
Print["Pz-Function - d/dz"]
SetPrecision[ pzZ[P0x, P0y, P0z],12]
SetPrecision[ pzZ[P1x, P1y, P1z],12]
SetPrecision[ pzZ[P2x, P2y, P2z],12]
SetPrecision[ pzZ[P3x, P3y, P3z],12]
Print["Pz-Function - d^2/dx^2"]
SetPrecision[ pzXX[P0x, P0y, P0z],12]
SetPrecision[ pzXX[P1x, P1y, P1z],12]
SetPrecision[ pzXX[P2x, P2y, P2z],12]
SetPrecision[ pzXX[P3x, P3y, P3z],12]
Print["Pz-Function - d^2/dxdy"]
SetPrecision[ pzXY[P0x, P0y, P0z],12]
SetPrecision[ pzXY[P1x, P1y, P1z],12]
SetPrecision[ pzXY[P2x, P2y, P2z],12]
SetPrecision[ pzXY[P3x, P3y, P3z],12]
Print["Pz-Function - d^2/dxdz"]
SetPrecision[ pzXZ[P0x, P0y, P0z],12]
SetPrecision[ pzXZ[P1x, P1y, P1z],12]
SetPrecision[ pzXZ[P2x, P2y, P2z],12]
SetPrecision[ pzXZ[P3x, P3y, P3z],12]
Print["Pz-Function - d^2/dy^2"]
SetPrecision[ pzYY[P0x, P0y, P0z],12]
SetPrecision[ pzYY[P1x, P1y, P1z],12]
SetPrecision[ pzYY[P2x, P2y, P2z],12]
SetPrecision[ pzYY[P3x, P3y, P3z],12]
Print["Pz-Function - d^2/dydz"]
SetPrecision[ pzYZ[P0x, P0y, P0z],12]
SetPrecision[ pzYZ[P1x, P1y, P1z],12]
SetPrecision[ pzYZ[P2x, P2y, P2z],12]
SetPrecision[ pzYZ[P3x, P3y, P3z],12]
Print["Pz-Function - d^2/dz^2"]
SetPrecision[ pzZZ[P0x, P0y, P0z],12]
SetPrecision[ pzZZ[P1x, P1y, P1z],12]
SetPrecision[ pzZZ[P2x, P2y, P2z],12]
SetPrecision[ pzZZ[P3x, P3y, P3z],12]
Print["Dxx-Function"]
SetPrecision[ dxx[P0x, P0y, P0z],12]
SetPrecision[ dxx[P1x, P1y, P1z],12]
SetPrecision[ dxx[P2x, P2y, P2z],12]
SetPrecision[ dxx[P3x, P3y, P3z],12]
Print["Dxx-Function - d/dx"]
SetPrecision[ dxxX[P0x, P0y, P0z],12]
SetPrecision[ dxxX[P1x, P1y, P1z],12]
SetPrecision[ dxxX[P2x, P2y, P2z],12]
SetPrecision[ dxxX[P3x, P3y, P3z],12]
Print["Dxx-Function - d/dy"]
SetPrecision[ dxxY[P0x, P0y, P0z],12]
SetPrecision[ dxxY[P1x, P1y, P1z],12]
SetPrecision[ dxxY[P2x, P2y, P2z],12]
SetPrecision[ dxxY[P3x, P3y, P3z],12]
Print["Dxx-Function - d/dz"]
SetPrecision[ dxxZ[P0x, P0y, P0z],12]
SetPrecision[ dxxZ[P1x, P1y, P1z],12]
SetPrecision[ dxxZ[P2x, P2y, P2z],12]
SetPrecision[ dxxZ[P3x, P3y, P3z],12]
Print["Dxx-Function - d^2/dx^2"]
SetPrecision[ dxxXX[P0x, P0y, P0z],12]
SetPrecision[ dxxXX[P1x, P1y, P1z],12]
SetPrecision[ dxxXX[P2x, P2y, P2z],12]
SetPrecision[ dxxXX[P3x, P3y, P3z],12]
Print["Dxx-Function - d^2/dxdy"]
SetPrecision[ dxxXY[P0x, P0y, P0z],12]
SetPrecision[ dxxXY[P1x, P1y, P1z],12]
SetPrecision[ dxxXY[P2x, P2y, P2z],12]
SetPrecision[ dxxXY[P3x, P3y, P3z],12]
Print["Dxx-Function - d^2/dxdz"]
SetPrecision[ dxxXZ[P0x, P0y, P0z],12]
SetPrecision[ dxxXZ[P1x, P1y, P1z],12]
SetPrecision[ dxxXZ[P2x, P2y, P2z],12]
SetPrecision[ dxxXZ[P3x, P3y, P3z],12]
Print["Dxx-Function - d^2/dy^2"]
SetPrecision[ dxxYY[P0x, P0y, P0z],12]
SetPrecision[ dxxYY[P1x, P1y, P1z],12]
SetPrecision[ dxxYY[P2x, P2y, P2z],12]
SetPrecision[ dxxYY[P3x, P3y, P3z],12]
Print["Dxx-Function - d^2/dydz"]
SetPrecision[ dxxYZ[P0x, P0y, P0z],12]
SetPrecision[ dxxYZ[P1x, P1y, P1z],12]
SetPrecision[ dxxYZ[P2x, P2y, P2z],12]
SetPrecision[ dxxYZ[P3x, P3y, P3z],12]
Print["Dxx-Function - d^2/dz^2"]
SetPrecision[ dxxZZ[P0x, P0y, P0z],12]
SetPrecision[ dxxZZ[P1x, P1y, P1z],12]
SetPrecision[ dxxZZ[P2x, P2y, P2z],12]
SetPrecision[ dxxZZ[P3x, P3y, P3z],12]
Print["Dxy-Function"]
SetPrecision[ dxy[P0x, P0y, P0z],12]
SetPrecision[ dxy[P1x, P1y, P1z],12]
SetPrecision[ dxy[P2x, P2y, P2z],12]
SetPrecision[ dxy[P3x, P3y, P3z],12]
Print["Dxy-Function - d/dx"]
SetPrecision[ dxyX[P0x, P0y, P0z],12]
SetPrecision[ dxyX[P1x, P1y, P1z],12]
SetPrecision[ dxyX[P2x, P2y, P2z],12]
SetPrecision[ dxyX[P3x, P3y, P3z],12]
Print["Dxy-Function - d/dy"]
SetPrecision[ dxyY[P0x, P0y, P0z],12]
SetPrecision[ dxyY[P1x, P1y, P1z],12]
SetPrecision[ dxyY[P2x, P2y, P2z],12]
SetPrecision[ dxyY[P3x, P3y, P3z],12]
Print["Dxy-Function - d/dz"]
SetPrecision[ dxyZ[P0x, P0y, P0z],12]
SetPrecision[ dxyZ[P1x, P1y, P1z],12]
SetPrecision[ dxyZ[P2x, P2y, P2z],12]
SetPrecision[ dxyZ[P3x, P3y, P3z],12]
Print["Dxy-Function - d^2/dx^2"]
SetPrecision[ dxyXX[P0x, P0y, P0z],12]
SetPrecision[ dxyXX[P1x, P1y, P1z],12]
SetPrecision[ dxyXX[P2x, P2y, P2z],12]
SetPrecision[ dxyXX[P3x, P3y, P3z],12]
Print["Dxy-Function - d^2/dxdy"]
SetPrecision[ dxyXY[P0x, P0y, P0z],12]
SetPrecision[ dxyXY[P1x, P1y, P1z],12]
SetPrecision[ dxyXY[P2x, P2y, P2z],12]
SetPrecision[ dxyXY[P3x, P3y, P3z],12]
Print["Dxy-Function - d^2/dxdz"]
SetPrecision[ dxyXZ[P0x, P0y, P0z],12]
SetPrecision[ dxyXZ[P1x, P1y, P1z],12]
SetPrecision[ dxyXZ[P2x, P2y, P2z],12]
SetPrecision[ dxyXZ[P3x, P3y, P3z],12]
Print["Dxy-Function - d^2/dy^2"]
SetPrecision[ dxyYY[P0x, P0y, P0z],12]
SetPrecision[ dxyYY[P1x, P1y, P1z],12]
SetPrecision[ dxyYY[P2x, P2y, P2z],12]
SetPrecision[ dxyYY[P3x, P3y, P3z],12]
Print["Dxy-Function - d^2/dydz"]
SetPrecision[ dxyYZ[P0x, P0y, P0z],12]
SetPrecision[ dxyYZ[P1x, P1y, P1z],12]
SetPrecision[ dxyYZ[P2x, P2y, P2z],12]
SetPrecision[ dxyYZ[P3x, P3y, P3z],12]
Print["Dxy-Function - d^2/dz^2"]
SetPrecision[ dxyZZ[P0x, P0y, P0z],12]
SetPrecision[ dxyZZ[P1x, P1y, P1z],12]
SetPrecision[ dxyZZ[P2x, P2y, P2z],12]
SetPrecision[ dxyZZ[P3x, P3y, P3z],12]
Print["Dxz-Function"]
SetPrecision[ dxz[P0x, P0y, P0z],12]
SetPrecision[ dxz[P1x, P1y, P1z],12]
SetPrecision[ dxz[P2x, P2y, P2z],12]
SetPrecision[ dxz[P3x, P3y, P3z],12]
Print["Dxz-Function - d/dx"]
SetPrecision[ dxzX[P0x, P0y, P0z],12]
SetPrecision[ dxzX[P1x, P1y, P1z],12]
SetPrecision[ dxzX[P2x, P2y, P2z],12]
SetPrecision[ dxzX[P3x, P3y, P3z],12]
Print["Dxz-Function - d/dy"]
SetPrecision[ dxzY[P0x, P0y, P0z],12]
SetPrecision[ dxzY[P1x, P1y, P1z],12]
SetPrecision[ dxzY[P2x, P2y, P2z],12]
SetPrecision[ dxzY[P3x, P3y, P3z],12]
Print["Dxz-Function - d/dz"]
SetPrecision[ dxzZ[P0x, P0y, P0z],12]
SetPrecision[ dxzZ[P1x, P1y, P1z],12]
SetPrecision[ dxzZ[P2x, P2y, P2z],12]
SetPrecision[ dxzZ[P3x, P3y, P3z],12]
Print["Dxz-Function - d^2/dx^2"]
SetPrecision[ dxzXX[P0x, P0y, P0z],12]
SetPrecision[ dxzXX[P1x, P1y, P1z],12]
SetPrecision[ dxzXX[P2x, P2y, P2z],12]
SetPrecision[ dxzXX[P3x, P3y, P3z],12]
Print["Dxz-Function - d^2/dxdy"]
SetPrecision[ dxzXY[P0x, P0y, P0z],12]
SetPrecision[ dxzXY[P1x, P1y, P1z],12]
SetPrecision[ dxzXY[P2x, P2y, P2z],12]
SetPrecision[ dxzXY[P3x, P3y, P3z],12]
Print["Dxz-Function - d^2/dxdz"]
SetPrecision[ dxzXZ[P0x, P0y, P0z],12]
SetPrecision[ dxzXZ[P1x, P1y, P1z],12]
SetPrecision[ dxzXZ[P2x, P2y, P2z],12]
SetPrecision[ dxzXZ[P3x, P3y, P3z],12]
Print["Dxz-Function - d^2/dy^2"]
SetPrecision[ dxzYY[P0x, P0y, P0z],12]
SetPrecision[ dxzYY[P1x, P1y, P1z],12]
SetPrecision[ dxzYY[P2x, P2y, P2z],12]
SetPrecision[ dxzYY[P3x, P3y, P3z],12]
Print["Dxz-Function - d^2/dydz"]
SetPrecision[ dxzYZ[P0x, P0y, P0z],12]
SetPrecision[ dxzYZ[P1x, P1y, P1z],12]
SetPrecision[ dxzYZ[P2x, P2y, P2z],12]
SetPrecision[ dxzYZ[P3x, P3y, P3z],12]
Print["Dxz-Function - d^2/dz^2"]
SetPrecision[ dxzZZ[P0x, P0y, P0z],12]
SetPrecision[ dxzZZ[P1x, P1y, P1z],12]
SetPrecision[ dxzZZ[P2x, P2y, P2z],12]
SetPrecision[ dxzZZ[P3x, P3y, P3z],12]
Print["Dyy-Function"]
SetPrecision[ dyy[P0x, P0y, P0z],12]
SetPrecision[ dyy[P1x, P1y, P1z],12]
SetPrecision[ dyy[P2x, P2y, P2z],12]
SetPrecision[ dyy[P3x, P3y, P3z],12]
Print["Dyy-Function - d/dx"]
SetPrecision[ dyyX[P0x, P0y, P0z],12]
SetPrecision[ dyyX[P1x, P1y, P1z],12]
SetPrecision[ dyyX[P2x, P2y, P2z],12]
SetPrecision[ dyyX[P3x, P3y, P3z],12]
Print["Dyy-Function - d/dy"]
SetPrecision[ dyyY[P0x, P0y, P0z],12]
SetPrecision[ dyyY[P1x, P1y, P1z],12]
SetPrecision[ dyyY[P2x, P2y, P2z],12]
SetPrecision[ dyyY[P3x, P3y, P3z],12]
Print["Dyy-Function - d/dz"]
SetPrecision[ dyyZ[P0x, P0y, P0z],12]
SetPrecision[ dyyZ[P1x, P1y, P1z],12]
SetPrecision[ dyyZ[P2x, P2y, P2z],12]
SetPrecision[ dyyZ[P3x, P3y, P3z],12]
Print["Dyy-Function - d^2/dx^2"]
SetPrecision[ dyyXX[P0x, P0y, P0z],12]
SetPrecision[ dyyXX[P1x, P1y, P1z],12]
SetPrecision[ dyyXX[P2x, P2y, P2z],12]
SetPrecision[ dyyXX[P3x, P3y, P3z],12]
Print["Dyy-Function - d^2/dxdy"]
SetPrecision[ dyyXY[P0x, P0y, P0z],12]
SetPrecision[ dyyXY[P1x, P1y, P1z],12]
SetPrecision[ dyyXY[P2x, P2y, P2z],12]
SetPrecision[ dyyXY[P3x, P3y, P3z],12]
Print["Dyy-Function - d^2/dxdz"]
SetPrecision[ dyyXZ[P0x, P0y, P0z],12]
SetPrecision[ dyyXZ[P1x, P1y, P1z],12]
SetPrecision[ dyyXZ[P2x, P2y, P2z],12]
SetPrecision[ dyyXZ[P3x, P3y, P3z],12]
Print["Dyy-Function - d^2/dy^2"]
SetPrecision[ dyyYY[P0x, P0y, P0z],12]
SetPrecision[ dyyYY[P1x, P1y, P1z],12]
SetPrecision[ dyyYY[P2x, P2y, P2z],12]
SetPrecision[ dyyYY[P3x, P3y, P3z],12]
Print["Dyy-Function - d^2/dydz"]
SetPrecision[ dyyYZ[P0x, P0y, P0z],12]
SetPrecision[ dyyYZ[P1x, P1y, P1z],12]
SetPrecision[ dyyYZ[P2x, P2y, P2z],12]
SetPrecision[ dyyYZ[P3x, P3y, P3z],12]
Print["Dyy-Function - d^2/dz^2"]
SetPrecision[ dyyZZ[P0x, P0y, P0z],12]
SetPrecision[ dyyZZ[P1x, P1y, P1z],12]
SetPrecision[ dyyZZ[P2x, P2y, P2z],12]
SetPrecision[ dyyZZ[P3x, P3y, P3z],12]
Print["Dyz-Function"]
SetPrecision[ dyz[P0x, P0y, P0z],12]
SetPrecision[ dyz[P1x, P1y, P1z],12]
SetPrecision[ dyz[P2x, P2y, P2z],12]
SetPrecision[ dyz[P3x, P3y, P3z],12]
Print["Dyz-Function - d/dx"]
SetPrecision[ dyzX[P0x, P0y, P0z],12]
SetPrecision[ dyzX[P1x, P1y, P1z],12]
SetPrecision[ dyzX[P2x, P2y, P2z],12]
SetPrecision[ dyzX[P3x, P3y, P3z],12]
Print["Dyz-Function - d/dy"]
SetPrecision[ dyzY[P0x, P0y, P0z],12]
SetPrecision[ dyzY[P1x, P1y, P1z],12]
SetPrecision[ dyzY[P2x, P2y, P2z],12]
SetPrecision[ dyzY[P3x, P3y, P3z],12]
Print["Dyz-Function - d/dz"]
SetPrecision[ dyzZ[P0x, P0y, P0z],12]
SetPrecision[ dyzZ[P1x, P1y, P1z],12]
SetPrecision[ dyzZ[P2x, P2y, P2z],12]
SetPrecision[ dyzZ[P3x, P3y, P3z],12]
Print["Dyz-Function - d^2/dx^2"]
SetPrecision[ dyzXX[P0x, P0y, P0z],12]
SetPrecision[ dyzXX[P1x, P1y, P1z],12]
SetPrecision[ dyzXX[P2x, P2y, P2z],12]
SetPrecision[ dyzXX[P3x, P3y, P3z],12]
Print["Dyz-Function - d^2/dxdy"]
SetPrecision[ dyzXY[P0x, P0y, P0z],12]
SetPrecision[ dyzXY[P1x, P1y, P1z],12]
SetPrecision[ dyzXY[P2x, P2y, P2z],12]
SetPrecision[ dyzXY[P3x, P3y, P3z],12]
Print["Dyz-Function - d^2/dxdz"]
SetPrecision[ dyzXZ[P0x, P0y, P0z],12]
SetPrecision[ dyzXZ[P1x, P1y, P1z],12]
SetPrecision[ dyzXZ[P2x, P2y, P2z],12]
SetPrecision[ dyzXZ[P3x, P3y, P3z],12]
Print["Dyz-Function - d^2/dy^2"]
SetPrecision[ dyzYY[P0x, P0y, P0z],12]
SetPrecision[ dyzYY[P1x, P1y, P1z],12]
SetPrecision[ dyzYY[P2x, P2y, P2z],12]
SetPrecision[ dyzYY[P3x, P3y, P3z],12]
Print["Dyz-Function - d^2/dydz"]
SetPrecision[ dyzYZ[P0x, P0y, P0z],12]
SetPrecision[ dyzYZ[P1x, P1y, P1z],12]
SetPrecision[ dyzYZ[P2x, P2y, P2z],12]
SetPrecision[ dyzYZ[P3x, P3y, P3z],12]
Print["Dyz-Function - d^2/dz^2"]
SetPrecision[ dyzZZ[P0x, P0y, P0z],12]
SetPrecision[ dyzZZ[P1x, P1y, P1z],12]
SetPrecision[ dyzZZ[P2x, P2y, P2z],12]
SetPrecision[ dyzZZ[P3x, P3y, P3z],12]
Print["Dzz-Function"]
SetPrecision[ dzz[P0x, P0y, P0z],12]
SetPrecision[ dzz[P1x, P1y, P1z],12]
SetPrecision[ dzz[P2x, P2y, P2z],12]
SetPrecision[ dzz[P3x, P3y, P3z],12]
Print["Dzz-Function - d/dx"]
SetPrecision[ dzzX[P0x, P0y, P0z],12]
SetPrecision[ dzzX[P1x, P1y, P1z],12]
SetPrecision[ dzzX[P2x, P2y, P2z],12]
SetPrecision[ dzzX[P3x, P3y, P3z],12]
Print["Dzz-Function - d/dy"]
SetPrecision[ dzzY[P0x, P0y, P0z],12]
SetPrecision[ dzzY[P1x, P1y, P1z],12]
SetPrecision[ dzzY[P2x, P2y, P2z],12]
SetPrecision[ dzzY[P3x, P3y, P3z],12]
Print["Dzz-Function - d/dz"]
SetPrecision[ dzzZ[P0x, P0y, P0z],12]
SetPrecision[ dzzZ[P1x, P1y, P1z],12]
SetPrecision[ dzzZ[P2x, P2y, P2z],12]
SetPrecision[ dzzZ[P3x, P3y, P3z],12]
Print["Dzz-Function - d^2/dx^2"]
SetPrecision[ dzzXX[P0x, P0y, P0z],12]
SetPrecision[ dzzXX[P1x, P1y, P1z],12]
SetPrecision[ dzzXX[P2x, P2y, P2z],12]
SetPrecision[ dzzXX[P3x, P3y, P3z],12]
Print["Dzz-Function - d^2/dxdy"]
SetPrecision[ dzzXY[P0x, P0y, P0z],12]
SetPrecision[ dzzXY[P1x, P1y, P1z],12]
SetPrecision[ dzzXY[P2x, P2y, P2z],12]
SetPrecision[ dzzXY[P3x, P3y, P3z],12]
Print["Dzz-Function - d^2/dxdz"]
SetPrecision[ dzzXZ[P0x, P0y, P0z],12]
SetPrecision[ dzzXZ[P1x, P1y, P1z],12]
SetPrecision[ dzzXZ[P2x, P2y, P2z],12]
SetPrecision[ dzzXZ[P3x, P3y, P3z],12]
Print["Dzz-Function - d^2/dy^2"]
SetPrecision[ dzzYY[P0x, P0y, P0z],12]
SetPrecision[ dzzYY[P1x, P1y, P1z],12]
SetPrecision[ dzzYY[P2x, P2y, P2z],12]
SetPrecision[ dzzYY[P3x, P3y, P3z],12]
Print["Dzz-Function - d^2/dydz"]
SetPrecision[ dzzYZ[P0x, P0y, P0z],12]
SetPrecision[ dzzYZ[P1x, P1y, P1z],12]
SetPrecision[ dzzYZ[P2x, P2y, P2z],12]
SetPrecision[ dzzYZ[P3x, P3y, P3z],12]
Print["Dzz-Function - d^2/dz^2"]
SetPrecision[ dzzZZ[P0x, P0y, P0z],12]
SetPrecision[ dzzZZ[P1x, P1y, P1z],12]
SetPrecision[ dzzZZ[P2x, P2y, P2z],12]
SetPrecision[ dzzZZ[P3x, P3y, P3z],12]
During evaluation of In[5306]:= S-Function
Out[5481]= 0.832567765386
Out[5482]= 0.278411287801
Out[5483]= 0.00176735669536
Out[5484]= 0.0000251996341204
During evaluation of In[5306]:= S-Function - d/dx
Out[5486]= 0
Out[5487]= -0.589265770897
Out[5488]= 0.00707237261901
Out[5489]= -0.000151198703799
During evaluation of In[5306]:= S-Function - d/dy
Out[5491]= 0
Out[5492]= 0
Out[5493]= -0.00353618630950
Out[5494]= -0.0000251997839666
During evaluation of In[5306]:= S-Function - d/dz
Out[5496]= 0
Out[5497]= 0
Out[5498]= -0.00353618630950
Out[5499]= 0.0000503995679331
During evaluation of In[5306]:= S-Function - d^2/dx^2
Out[5501]= -1.90486012083
Out[5502]= 0.719038552076
Out[5503]= 0.0247768708671
Out[5504]= 0.000856803443789
During evaluation of In[5306]:= S-Function - d^2/dxdy
Out[5506]= 0
Out[5507]= 0
Out[5508]= -0.0141565285883
Out[5509]= 0.000151200501954
During evaluation of In[5306]:= S-Function - d^2/dxdz
Out[5511]= 0
Out[5512]= 0
Out[5513]= -0.0141565285883
Out[5514]= -0.000302401003907
During evaluation of In[5306]:= S-Function - d^2/dy^2
Out[5516]= -1.90486012083
Out[5517]= -0.589265770897
Out[5518]= 0.00354207798465
Out[5519]= -0.0000251994842742
During evaluation of In[5306]:= S-Function - d^2/dydz
Out[5521]= 0
Out[5522]= 0
Out[5523]= 0.00707826429416
Out[5524]= -0.0000504001673179
During evaluation of In[5306]:= S-Function - d^2/dz^2
Out[5526]= -1.90486012083
Out[5527]= -0.589265770897
Out[5528]= 0.00354207798465
Out[5529]= 0.0000504007667026
During evaluation of In[5306]:= Px-Function
Out[5531]= -0.570260987101
Out[5532]= 0
Out[5533]= -0.0000714206448118
Out[5534]= 0.0149784111869
During evaluation of In[5306]:= Px-Function - d/dx
Out[5536]= -0.662024200686
Out[5537]= 1.76443270721
Out[5538]= -0.000404718689509
Out[5539]= -0.0524991227999
During evaluation of In[5306]:= Px-Function - d/dy
Out[5541]= 0
Out[5542]= 0
Out[5543]= 0.000142841857038
Out[5544]= -0.0149970820983
During evaluation of In[5306]:= Px-Function - d/dz
Out[5546]= 0
Out[5547]= 0
Out[5548]= 0.000142841857038
Out[5549]= 0.0299941641967
During evaluation of In[5306]:= Px-Function - d^2/dx^2
Out[5551]= 0.865232333446
Out[5552]= 0
Out[5553]= -0.00214264828248
Out[5554]= 0.150568290148
During evaluation of In[5306]:= Px-Function - d^2/dxdy
Out[5556]= 0
Out[5557]= 0
Out[5558]= 0.000809443998852
Out[5559]= 0.0526391546354
During evaluation of In[5306]:= Px-Function - d^2/dxdz
Out[5561]= 0
Out[5562]= 0
Out[5563]= 0.000809443998852
Out[5564]= -0.105278309271
During evaluation of In[5306]:= Px-Function - d^2/dy^2
Out[5566]= 1.23228518779
Out[5567]= 0
Out[5568]= -0.000142844126695
Out[5569]= -0.0149597402755
During evaluation of In[5306]:= Px-Function - d^2/dydz
Out[5571]= 0
Out[5572]= 0
Out[5573]= -0.000285685983732
Out[5574]= -0.0300688478423
During evaluation of In[5306]:= Px-Function - d^2/dz^2
Out[5576]= 1.23228518779
Out[5577]= 0
Out[5578]= -0.000142844126695
Out[5579]= 0.0301435314879
During evaluation of In[5306]:= Py-Function
Out[5581]= 0
Out[5582]= 0
Out[5583]= 0.0000238068816039
Out[5584]= 0.00374460279674
During evaluation of In[5306]:= Py-Function - d/dx
Out[5586]= 0
Out[5587]= 0
Out[5588]= 0.000142841857038
Out[5589]= -0.0149970820983
During evaluation of In[5306]:= Py-Function - d/dy
Out[5591]= 0.570260987101
Out[5592]= 1.76443270721
Out[5593]= -0.0000238070707420
Out[5594]= 0.00373993506889
During evaluation of In[5306]:= Py-Function - d/dz
Out[5596]= 0
Out[5597]= 0
Out[5598]= -0.0000476139523459
Out[5599]= 0.00749854104917
During evaluation of In[5306]:= Py-Function - d^2/dx^2
Out[5601]= 0
Out[5602]= 0
Out[5603]= 0.000809443998852
Out[5604]= 0.0526391546354
During evaluation of In[5306]:= Py-Function - d^2/dxdy
Out[5606]= 1.23228518779
Out[5607]= 0
Out[5608]= -0.000142844126695
Out[5609]= -0.0149597402755
During evaluation of In[5306]:= Py-Function - d^2/dxdz
Out[5611]= 0
Out[5612]= 0
Out[5613]= -0.000285685983732
Out[5614]= -0.0300688478423
During evaluation of In[5306]:= Py-Function - d^2/dy^2
Out[5616]= 0
Out[5617]= 0
Out[5618]= -0.0000476131957936
Out[5619]= -0.0187370171672
During evaluation of In[5306]:= Py-Function - d^2/dydz
Out[5621]= 0
Out[5622]= 0
Out[5623]= 0.0000476147088982
Out[5624]= 0.00747987013777
During evaluation of In[5306]:= Py-Function - d^2/dz^2
Out[5626]= 0
Out[5627]= 0
Out[5628]= 0.0000476147088982
Out[5629]= 0.00753588287198
During evaluation of In[5306]:= Pz-Function
Out[5631]= 0
Out[5632]= 0
Out[5633]= 0.0000238068816039
Out[5634]= -0.00748920559347
During evaluation of In[5306]:= Pz-Function - d/dx
Out[5636]= 0
Out[5637]= 0
Out[5638]= 0.000142841857038
Out[5639]= 0.0299941641967
During evaluation of In[5306]:= Pz-Function - d/dy
Out[5641]= 0
Out[5642]= 0
Out[5643]= -0.0000476139523459
Out[5644]= 0.00749854104917
During evaluation of In[5306]:= Pz-Function - d/dz
Out[5646]= 0.570260987101
Out[5647]= 1.76443270721
Out[5648]= -0.0000238070707420
Out[5649]= -0.00750787650488
During evaluation of In[5306]:= Pz-Function - d^2/dx^2
Out[5651]= 0
Out[5652]= 0
Out[5653]= 0.000809443998852
Out[5654]= -0.105278309271
During evaluation of In[5306]:= Pz-Function - d^2/dxdy
Out[5656]= 0
Out[5657]= 0
Out[5658]= -0.000285685983732
Out[5659]= -0.0300688478423
During evaluation of In[5306]:= Pz-Function - d^2/dxdz
Out[5661]= 1.23228518779
Out[5662]= 0
Out[5663]= -0.000142844126695
Out[5664]= 0.0301435314879
During evaluation of In[5306]:= Pz-Function - d^2/dy^2
Out[5666]= 0
Out[5667]= 0
Out[5668]= 0.0000476147088982
Out[5669]= 0.00747987013777
During evaluation of In[5306]:= Pz-Function - d^2/dydz
Out[5671]= 0
Out[5672]= 0
Out[5673]= 0.0000476147088982
Out[5674]= 0.00753588287198
During evaluation of In[5306]:= Pz-Function - d^2/dz^2
Out[5676]= 0
Out[5677]= 0
Out[5678]= -0.0000476131957936
Out[5679]= 0.0149223984527
During evaluation of In[5306]:= Dxx-Function
Out[5681]= 0.00637294953921
Out[5682]= 0.00565239243008
Out[5683]= 0.0954200464484
Out[5684]= 0.0000384842654340
During evaluation of In[5306]:= Dxx-Function - d/dx
Out[5686]= 0.00630932138969
Out[5687]= -0.00566630089326
Out[5688]= -0.0697506188560
Out[5689]= -0.000134694954234
During evaluation of In[5306]:= Dxx-Function - d/dy
Out[5691]= 0.0160914442218
Out[5692]= 0.0141483666542
Out[5693]= 0.181634211061
Out[5694]= 0.0000769685434756
During evaluation of In[5306]:= Dxx-Function - d/dz
Out[5696]= 0.00643657768873
Out[5697]= 0.00565934666167
Out[5698]= 0
Out[5699]= 0.0000769685434756
During evaluation of In[5306]:= Dxx-Function - d^2/dx^2
Out[5701]= -0.0128477041176
Out[5702]= -0.00276638982337
Out[5703]= -0.231911177045
Out[5704]= 0.000428137764361
During evaluation of In[5306]:= Dxx-Function - d^2/dxdy
Out[5706]= 0.0157096753247
Out[5707]= -0.0142318174332
Out[5708]= -0.0892242717285
Out[5709]= -0.000269390023198
During evaluation of In[5306]:= Dxx-Function - d^2/dxdz
Out[5711]= 0.0157096753247
Out[5712]= -0.0142318174332
Out[5713]= -0.0892242717285
Out[5714]= -0.000269390023198
During evaluation of In[5306]:= Dxx-Function - d^2/dy^2
Out[5716]= 0.0347464551086
Out[5717]= 0.0298158834476
Out[5718]= 0.289976751550
Out[5719]= 0.000115452875730
During evaluation of In[5306]:= Dxx-Function - d^2/dydz
Out[5721]= 0.0164732131189
Out[5722]= 0.0141900920437
Out[5723]= 0
Out[5724]= 0.000153937147468
During evaluation of In[5306]:= Dxx-Function - d^2/dz^2
Out[5726]= 0.000152707558845
Out[5727]= 0.0000166901558139
Out[5728]= -0.121089474041
Out[5729]= 0.000115452875730
During evaluation of In[5306]:= Dxy-Function
Out[5731]= -0.0275956809900
Out[5732]= -0.0122377885915
Out[5733]= 0.247908552764
Out[5734]= -0.0000333283515118
During evaluation of In[5306]:= Dxy-Function - d/dx
Out[5736]= 0.000275517969394
Out[5737]= 0.0183867955934
Out[5738]= 0.0666911291870
Out[5739]= 0.000124981340006
During evaluation of In[5306]:= Dxy-Function - d/dy
Out[5741]= -0.0586397250024
Out[5742]= -0.0257369969248
Out[5743]= 0.306627154417
Out[5744]= -0.0000499925381862
During evaluation of In[5306]:= Dxy-Function - d/dz
Out[5746]= -0.0278711989593
Out[5747]= -0.0122528449446
Out[5748]= 0
Out[5749]= -0.0000666567139421
During evaluation of In[5306]:= Dxy-Function - d^2/dx^2
Out[5751]= 0.0550811547921
Out[5752]= -0.0123973859338
Out[5753]= -0.469140653853
Out[5754]= -0.000433268850260
During evaluation of In[5306]:= Dxy-Function - d^2/dxdy
Out[5756]= 0.00154290062861
Out[5757]= 0.0387741265414
Out[5758]= 0.195627312282
Out[5759]= 0.000187472112098
During evaluation of In[5306]:= Dxy-Function - d^2/dxdz
Out[5761]= 0.00154290062861
Out[5762]= 0.0387741265414
Out[5763]= 0.195627312282
Out[5764]= 0.000187472112098
During evaluation of In[5306]:= Dxy-Function - d^2/dy^2
Out[5766]= -0.0947141661588
Out[5767]= -0.0400475913656
Out[5768]= 0.124182336147
Out[5769]= -0.0000333284093800
During evaluation of In[5306]:= Dxy-Function - d^2/dydz
Out[5771]= -0.0601826256310
Out[5772]= -0.0258213125019
Out[5773]= 0
Out[5774]= -0.0000999851233221
During evaluation of In[5306]:= Dxy-Function - d^2/dz^2
Out[5776]= -0.000661243126546
Out[5777]= -0.0000361352473199
Out[5778]= -0.314599681951
Out[5779]= -0.0000999851233221
During evaluation of In[5306]:= Dxz-Function
Out[5781]= -0.0110382723960
Out[5782]= -0.00489511543661
Out[5783]= 0
Out[5784]= -0.0000333283515118
During evaluation of In[5306]:= Dxz-Function - d/dx
Out[5786]= 0.000110207187758
Out[5787]= 0.00735471823735
Out[5788]= 0
Out[5789]= 0.000124981340006
During evaluation of In[5306]:= Dxz-Function - d/dy
Out[5791]= -0.0278711989593
Out[5792]= -0.0122528449446
Out[5793]= 0
Out[5794]= -0.0000666567139421
During evaluation of In[5306]:= Dxz-Function - d/dz
Out[5796]= -0.000110207187758
Out[5797]= -6.02254121998*10^-6
Out[5798]= -0.165272368509
Out[5799]= -0.0000499925381862
During evaluation of In[5306]:= Dxz-Function - d^2/dx^2
Out[5801]= 0.0220324619169
Out[5802]= -0.00495895437354
Out[5803]= 0
Out[5804]= -0.000433268850260
During evaluation of In[5306]:= Dxz-Function - d^2/dxdy
Out[5806]= 0.000661243126546
Out[5807]= 0.0184515379115
Out[5808]= 0
Out[5809]= 0.000249962782101
During evaluation of In[5306]:= Dxz-Function - d^2/dxdz
Out[5811]= 0.000661243126546
Out[5812]= 0.0184515379115
Out[5813]= 0
Out[5814]= 0.000249962782101
During evaluation of In[5306]:= Dxz-Function - d^2/dy^2
Out[5816]= -0.0601826256310
Out[5817]= -0.0258213125019
Out[5818]= 0
Out[5819]= -0.0000999851233221
During evaluation of In[5306]:= Dxz-Function - d^2/dydz
Out[5821]= -0.000661243126546
Out[5822]= -0.0000361352473199
Out[5823]= -0.314599681951
Out[5824]= -0.0000999851233221
During evaluation of In[5306]:= Dxz-Function - d^2/dz^2
Out[5826]= 0.0220324619169
Out[5827]= 0.00978782185672
Out[5828]= 0
Out[5829]= -0.0000333284093800
During evaluation of In[5306]:= Dyy-Function
Out[5831]= 0.0689892024749
Out[5832]= 0.0152972357394
Out[5833]= 0.371862829146
Out[5834]= 0.0000166641757559
During evaluation of In[5306]:= Dyy-Function - d/dx
Out[5836]= -0.0696779973984
Out[5837]= -0.0306321123614
Out[5838]= 0.471899522926
Out[5839]= -0.0000666567139421
During evaluation of In[5306]:= Dyy-Function - d/dy
Out[5841]= 0.119003631516
Out[5842]= 0.0260523518603
Out[5843]= 0.212032178862
Out[5844]= 0.0000166641812152
During evaluation of In[5306]:= Dyy-Function - d/dz
Out[5846]= 0.0696779973984
Out[5847]= 0.0153160561807
Out[5848]= 0
Out[5849]= 0.0000333283569711
During evaluation of In[5306]:= Dyy-Function - d^2/dx^2
Out[5851]= 0.00165310781637
Out[5852]= 0.0461288447787
Out[5853]= 0.240088065073
Out[5854]= 0.000249962782101
During evaluation of In[5306]:= Dyy-Function - d^2/dxdy
Out[5856]= -0.122585365118
Out[5857]= -0.0523004363102
Out[5858]= 0.438782018097
Out[5859]= -0.0000666567663511
During evaluation of In[5306]:= Dyy-Function - d^2/dxdz
Out[5861]= -0.122585365118
Out[5862]= -0.0523004363102
Out[5863]= 0.438782018097
Out[5864]= -0.0000666567663511
During evaluation of In[5306]:= Dyy-Function - d^2/dy^2
Out[5866]= 0.119505965392
Out[5867]= 0.0243224922822
Out[5868]= -0.426980804614
Out[5869]= -8.33206440311*10^-6
During evaluation of In[5306]:= Dyy-Function - d^2/dydz
Out[5871]= 0.122585365118
Out[5872]= 0.0261502181551
Out[5873]= 0
Out[5874]= 0.0000333283831755
During evaluation of In[5306]:= Dyy-Function - d^2/dz^2
Out[5876]= 0.00165310781637
Out[5877]= 0.0000451690591498
Out[5878]= -0.471899522926
Out[5879]= 0.0000499925616611
During evaluation of In[5306]:= Dyz-Function
Out[5881]= 0.0275956809900
Out[5882]= 0.00611889429576
Out[5883]= 0
Out[5884]= 0.0000166641757559
During evaluation of In[5306]:= Dyz-Function - d/dx
Out[5886]= -0.0278711989593
Out[5887]= -0.0122528449446
Out[5888]= 0
Out[5889]= -0.0000666567139421
During evaluation of In[5306]:= Dyz-Function - d/dy
Out[5891]= 0.0586397250024
Out[5892]= 0.0128684984624
Out[5893]= 0
Out[5894]= 0.0000249962690931
During evaluation of In[5306]:= Dyz-Function - d/dz
Out[5896]= 0.000275517969394
Out[5897]= 7.52817652497*10^-6
Out[5898]= -0.247908552764
Out[5899]= 0.0000249962690931
During evaluation of In[5306]:= Dyz-Function - d^2/dx^2
Out[5901]= 0.000661243126546
Out[5902]= 0.0184515379115
Out[5903]= 0
Out[5904]= 0.000249962782101
During evaluation of In[5306]:= Dyz-Function - d^2/dxdy
Out[5906]= -0.0601826256310
Out[5907]= -0.0258213125019
Out[5908]= 0
Out[5909]= -0.0000999851233221
During evaluation of In[5306]:= Dyz-Function - d^2/dxdz
Out[5911]= -0.0601826256310
Out[5912]= -0.0258213125019
Out[5913]= 0
Out[5914]= -0.0000999851233221
During evaluation of In[5306]:= Dyz-Function - d^2/dy^2
Out[5916]= 0.0947141661588
Out[5917]= 0.0200237956828
Out[5918]= 0
Out[5919]= 0.0000166642046900
During evaluation of In[5306]:= Dyz-Function - d^2/dydz
Out[5921]= 0.00154290062861
Out[5922]= 0.0000421577885398
Out[5923]= -0.306627154417
Out[5924]= 0.0000374944271145
During evaluation of In[5306]:= Dyz-Function - d^2/dz^2
Out[5926]= -0.0550811547921
Out[5927]= -0.0122347773209
Out[5928]= 0
Out[5929]= 0.0000166642046900
During evaluation of In[5306]:= Dzz-Function
Out[5931]= -0.0275956809900
Out[5932]= -0.0122377885915
Out[5933]= 0.247908552764
Out[5934]= -0.0000333283515118
During evaluation of In[5306]:= Dzz-Function - d/dx
Out[5936]= 0.000275517969394
Out[5937]= 0.0183867955934
Out[5938]= 0.0666911291870
Out[5939]= 0.000124981340006
During evaluation of In[5306]:= Dzz-Function - d/dy
Out[5941]= -0.0586397250024
Out[5942]= -0.0257369969248
Out[5943]= 0.306627154417
Out[5944]= -0.0000499925381862
During evaluation of In[5306]:= Dzz-Function - d/dz
Out[5946]= -0.0278711989593
Out[5947]= -0.0122528449446
Out[5948]= 0
Out[5949]= -0.0000666567139421
During evaluation of In[5306]:= Dzz-Function - d^2/dx^2
Out[5951]= 0.0550811547921
Out[5952]= -0.0123973859338
Out[5953]= -0.469140653853
Out[5954]= -0.000433268850260
During evaluation of In[5306]:= Dzz-Function - d^2/dxdy
Out[5956]= 0.00154290062861
Out[5957]= 0.0387741265414
Out[5958]= 0.195627312282
Out[5959]= 0.000187472112098
During evaluation of In[5306]:= Dzz-Function - d^2/dxdz
Out[5961]= 0.00154290062861
Out[5962]= 0.0387741265414
Out[5963]= 0.195627312282
Out[5964]= 0.000187472112098
During evaluation of In[5306]:= Dzz-Function - d^2/dy^2
Out[5966]= -0.0947141661588
Out[5967]= -0.0400475913656
Out[5968]= 0.124182336147
Out[5969]= -0.0000333284093800
During evaluation of In[5306]:= Dzz-Function - d^2/dydz
Out[5971]= -0.0601826256310
Out[5972]= -0.0258213125019
Out[5973]= 0
Out[5974]= -0.0000999851233221
During evaluation of In[5306]:= Dzz-Function - d^2/dz^2
Out[5976]= -0.000661243126546
Out[5977]= -0.0000361352473199
Out[5978]= -0.314599681951
Out[5979]= -0.0000999851233221