In[8002]:= (*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*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*(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_] := (y-y0d)^mdyy*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*mdyy-1)!!])
(*First primitive of the dyz-type function*)
ldyz:=0
mdyz:=1
ndyz:=1
d1yz[x_,y_,z_] := (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:=0
mdzz:=0
ndzz:=2
d1zz[x_,y_,z_] := (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*ndzz-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*mdyy-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*ndzz-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],z];
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],z];
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],z];
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],z];
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],z];
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],z];
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["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["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["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["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["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["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["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["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["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]
(*First derivative d/dx*)
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["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["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["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["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["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["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["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["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["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]
(*First derivative d/dy*)
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["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["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["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["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["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["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["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["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["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]
(*First derivative d/dz*)
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["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["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["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["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["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["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["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["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["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]
(*Second derivative d^2/dx^2*)
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["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["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["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["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["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["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["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["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["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]
(*Second derivative d^2/dxdy*)
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["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["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["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["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["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["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["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["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["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]
(*Second derivative d^2/dxdz*)
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["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["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["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["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["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["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["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["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["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]
(*Second derivative d^2/dy^2*)
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["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["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["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["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["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["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["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["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["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]
(*Second derivative d^2/dydz*)
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["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["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["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["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["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["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["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["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["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]
(*Second derivative d^2/dz^2*)
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 - 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 - 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 - 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 - 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 - 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 - 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 - 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 - 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 - 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[8002]:= S-Function
Out[8177]= 0.832567765386
Out[8178]= 0.278411287801
Out[8179]= 0.00176735669536
Out[8180]= 0.0000251996341204
During evaluation of In[8002]:= Px-Function
Out[8182]= -0.570260987101
Out[8183]= 0
Out[8184]= -0.0000714206448118
Out[8185]= 0.0149784111869
During evaluation of In[8002]:= Py-Function
Out[8187]= 0
Out[8188]= 0
Out[8189]= 0.0000238068816039
Out[8190]= 0.00374460279674
During evaluation of In[8002]:= Pz-Function
Out[8192]= 0
Out[8193]= 0
Out[8194]= 0.0000238068816039
Out[8195]= -0.00748920559347
During evaluation of In[8002]:= Dxx-Function
Out[8197]= 0.00637294953921
Out[8198]= 0.00565239243008
Out[8199]= 0.0954200464484
Out[8200]= 0.0000384842654340
During evaluation of In[8002]:= Dxy-Function
Out[8202]= -0.0275956809900
Out[8203]= -0.0122377885915
Out[8204]= 0.247908552764
Out[8205]= -0.0000333283515118
During evaluation of In[8002]:= Dxz-Function
Out[8207]= -0.0110382723960
Out[8208]= -0.00489511543661
Out[8209]= 0
Out[8210]= -0.0000333283515118
During evaluation of In[8002]:= Dyy-Function
Out[8212]= 0.0398309346201
Out[8213]= 0.00883186317200
Out[8214]= 0.214695104509
Out[8215]= 9.62106635849*10^-6
During evaluation of In[8002]:= Dyz-Function
Out[8217]= 0.0275956809900
Out[8218]= 0.00611889429576
Out[8219]= 0
Out[8220]= 0.0000166641757559
During evaluation of In[8002]:= Dzz-Function
Out[8222]= 0.00637294953921
Out[8223]= 0.00141309810752
Out[8224]= 0
Out[8225]= 9.62106635849*10^-6
During evaluation of In[8002]:= S-Function - d/dx
Out[8227]= 0
Out[8228]= -0.589265770897
Out[8229]= 0.00707237261901
Out[8230]= -0.000151198703799
During evaluation of In[8002]:= Px-Function - d/dx
Out[8232]= -0.662024200686
Out[8233]= 1.76443270721
Out[8234]= -0.000404718689509
Out[8235]= -0.0524991227999
During evaluation of In[8002]:= Py-Function - d/dx
Out[8237]= 0
Out[8238]= 0
Out[8239]= 0.000142841857038
Out[8240]= -0.0149970820983
During evaluation of In[8002]:= Pz-Function - d/dx
Out[8242]= 0
Out[8243]= 0
Out[8244]= 0.000142841857038
Out[8245]= 0.0299941641967
During evaluation of In[8002]:= Dxx-Function - d/dx
Out[8247]= 0.00630932138969
Out[8248]= -0.00566630089326
Out[8249]= -0.0697506188560
Out[8250]= -0.000134694954234
During evaluation of In[8002]:= Dxy-Function - d/dx
Out[8252]= 0.000275517969394
Out[8253]= 0.0183867955934
Out[8254]= 0.0666911291870
Out[8255]= 0.000124981340006
During evaluation of In[8002]:= Dxz-Function - d/dx
Out[8257]= 0.000110207187758
Out[8258]= 0.00735471823735
Out[8259]= 0
Out[8260]= 0.000124981340006
During evaluation of In[8002]:= Dyy-Function - d/dx
Out[8262]= -0.0402286105545
Out[8263]= -0.0176854583177
Out[8264]= 0.272451316592
Out[8265]= -0.0000384842717378
During evaluation of In[8002]:= Dyz-Function - d/dx
Out[8267]= -0.0278711989593
Out[8268]= -0.0122528449446
Out[8269]= 0
Out[8270]= -0.0000666567139421
During evaluation of In[8002]:= Dzz-Function - d/dx
Out[8272]= -0.00643657768873
Out[8273]= -0.00282967333083
Out[8274]= 0
Out[8275]= -0.0000384842717378
During evaluation of In[8002]:= S-Function - d/dy
Out[8277]= 0
Out[8278]= 0
Out[8279]= -0.00353618630950
Out[8280]= -0.0000251997839666
During evaluation of In[8002]:= Px-Function - d/dy
Out[8282]= 0
Out[8283]= 0
Out[8284]= 0.000142841857038
Out[8285]= -0.0149970820983
During evaluation of In[8002]:= Py-Function - d/dy
Out[8287]= 0.570260987101
Out[8288]= 1.76443270721
Out[8289]= -0.0000238070707420
Out[8290]= 0.00373993506889
During evaluation of In[8002]:= Pz-Function - d/dy
Out[8292]= 0
Out[8293]= 0
Out[8294]= -0.0000476139523459
Out[8295]= 0.00749854104917
During evaluation of In[8002]:= Dxx-Function - d/dy
Out[8297]= 0.0160914442218
Out[8298]= 0.0141483666542
Out[8299]= 0.181634211061
Out[8300]= 0.0000769685434756
During evaluation of In[8002]:= Dxy-Function - d/dy
Out[8302]= -0.0586397250024
Out[8303]= -0.0257369969248
Out[8304]= 0.306627154417
Out[8305]= -0.0000499925381862
During evaluation of In[8002]:= Dxz-Function - d/dy
Out[8307]= -0.0278711989593
Out[8308]= -0.0122528449446
Out[8309]= 0
Out[8310]= -0.0000666567139421
During evaluation of In[8002]:= Dyy-Function - d/dy
Out[8312]= 0.0687067786903
Out[8313]= 0.0150413323595
Out[8314]= 0.122416835543
Out[8315]= 9.62106951040*10^-6
During evaluation of In[8002]:= Dyz-Function - d/dy
Out[8317]= 0.0586397250024
Out[8318]= 0.0128684984624
Out[8319]= 0
Out[8320]= 0.0000249962690931
During evaluation of In[8002]:= Dzz-Function - d/dy
Out[8322]= 0.0160914442218
Out[8323]= 0.00353709166354
Out[8324]= 0
Out[8325]= 0.0000192421358689
During evaluation of In[8002]:= S-Function - d/dz
Out[8327]= 0
Out[8328]= 0
Out[8329]= -0.00353618630950
Out[8330]= 0.0000503995679331
During evaluation of In[8002]:= Px-Function - d/dz
Out[8332]= 0
Out[8333]= 0
Out[8334]= 0.000142841857038
Out[8335]= 0.0299941641967
During evaluation of In[8002]:= Py-Function - d/dz
Out[8337]= 0
Out[8338]= 0
Out[8339]= -0.0000476139523459
Out[8340]= 0.00749854104917
During evaluation of In[8002]:= Pz-Function - d/dz
Out[8342]= 0.570260987101
Out[8343]= 1.76443270721
Out[8344]= -0.0000238070707420
Out[8345]= -0.00750787650488
During evaluation of In[8002]:= Dxx-Function - d/dz
Out[8347]= 0.00643657768873
Out[8348]= 0.00565934666167
Out[8349]= 0
Out[8350]= 0.0000769685434756
During evaluation of In[8002]:= Dxy-Function - d/dz
Out[8352]= -0.0278711989593
Out[8353]= -0.0122528449446
Out[8354]= 0
Out[8355]= -0.0000666567139421
During evaluation of In[8002]:= Dxz-Function - d/dz
Out[8357]= -0.000110207187758
Out[8358]= -6.02254121998*10^-6
Out[8359]= -0.165272368509
Out[8360]= -0.0000499925381862
During evaluation of In[8002]:= Dyy-Function - d/dz
Out[8362]= 0.0402286105545
Out[8363]= 0.00884272915885
Out[8364]= 0
Out[8365]= 0.0000192421358689
During evaluation of In[8002]:= Dyz-Function - d/dz
Out[8367]= 0.000275517969394
Out[8368]= 7.52817652497*10^-6
Out[8369]= -0.247908552764
Out[8370]= 0.0000249962690931
During evaluation of In[8002]:= Dzz-Function - d/dz
Out[8372]= -0.00630932138969
Out[8373]= -0.00141135954962
Out[8374]= 0
Out[8375]= 9.62106951040*10^-6
During evaluation of In[8002]:= S-Function - d^2/dx^2
Out[8377]= -1.90486012083
Out[8378]= 0.719038552076
Out[8379]= 0.0247768708671
Out[8380]= 0.000856803443789
During evaluation of In[8002]:= Px-Function - d^2/dx^2
Out[8382]= 0.865232333446
Out[8383]= 0
Out[8384]= -0.00214264828248
Out[8385]= 0.150568290148
During evaluation of In[8002]:= Py-Function - d^2/dx^2
Out[8387]= 0
Out[8388]= 0
Out[8389]= 0.000809443998852
Out[8390]= 0.0526391546354
During evaluation of In[8002]:= Pz-Function - d^2/dx^2
Out[8392]= 0
Out[8393]= 0
Out[8394]= 0.000809443998852
Out[8395]= -0.105278309271
During evaluation of In[8002]:= Dxx-Function - d^2/dx^2
Out[8397]= -0.0128477041176
Out[8398]= -0.00276638982337
Out[8399]= -0.231911177045
Out[8400]= 0.000428137764361
During evaluation of In[8002]:= Dxy-Function - d^2/dx^2
Out[8402]= 0.0550811547921
Out[8403]= -0.0123973859338
Out[8404]= -0.469140653853
Out[8405]= -0.000433268850260
During evaluation of In[8002]:= Dxz-Function - d^2/dx^2
Out[8407]= 0.0220324619169
Out[8408]= -0.00495895437354
Out[8409]= 0
Out[8410]= -0.000433268850260
During evaluation of In[8002]:= Dyy-Function - d^2/dx^2
Out[8412]= 0.000954422242778
Out[8413]= 0.0266325009504
Out[8414]= 0.138614908999
Out[8415]= 0.000144316079533
During evaluation of In[8002]:= Dyz-Function - d^2/dx^2
Out[8417]= 0.000661243126546
Out[8418]= 0.0184515379115
Out[8419]= 0
Out[8420]= 0.000249962782101
During evaluation of In[8002]:= Dzz-Function - d^2/dx^2
Out[8422]= 0.000152707558845
Out[8423]= 0.00426120015206
Out[8424]= 0
Out[8425]= 0.000144316079533
During evaluation of In[8002]:= S-Function - d^2/dxdy
Out[8427]= 0
Out[8428]= 0
Out[8429]= -0.0141565285883
Out[8430]= 0.000151200501954
During evaluation of In[8002]:= Px-Function - d^2/dxdy
Out[8432]= 0
Out[8433]= 0
Out[8434]= 0.000809443998852
Out[8435]= 0.0526391546354
During evaluation of In[8002]:= Py-Function - d^2/dxdy
Out[8437]= 1.23228518779
Out[8438]= 0
Out[8439]= -0.000142844126695
Out[8440]= -0.0149597402755
During evaluation of In[8002]:= Pz-Function - d^2/dxdy
Out[8442]= 0
Out[8443]= 0
Out[8444]= -0.000285685983732
Out[8445]= -0.0300688478423
During evaluation of In[8002]:= Dxx-Function - d^2/dxdy
Out[8447]= 0.0157096753247
Out[8448]= -0.0142318174332
Out[8449]= -0.0892242717285
Out[8450]= -0.000269390023198
During evaluation of In[8002]:= Dxy-Function - d^2/dxdy
Out[8452]= 0.00154290062861
Out[8453]= 0.0387741265414
Out[8454]= 0.195627312282
Out[8455]= 0.000187472112098
During evaluation of In[8002]:= Dxz-Function - d^2/dxdy
Out[8457]= 0.000661243126546
Out[8458]= 0.0184515379115
Out[8459]= 0
Out[8460]= 0.000249962782101
During evaluation of In[8002]:= Dyy-Function - d^2/dxdy
Out[8462]= -0.0707746935497
Out[8463]= -0.0301956709824
Out[8464]= 0.253330916264
Out[8465]= -0.0000384843019961
During evaluation of In[8002]:= Dyz-Function - d^2/dxdy
Out[8467]= -0.0601826256310
Out[8468]= -0.0258213125019
Out[8469]= 0
Out[8470]= -0.0000999851233221
During evaluation of In[8002]:= Dzz-Function - d^2/dxdy
Out[8472]= -0.0164732131189
Out[8473]= -0.00709504602185
Out[8474]= 0
Out[8475]= -0.0000769685737339
During evaluation of In[8002]:= S-Function - d^2/dxdz
Out[8477]= 0
Out[8478]= 0
Out[8479]= -0.0141565285883
Out[8480]= -0.000302401003907
During evaluation of In[8002]:= Px-Function - d^2/dxdz
Out[8482]= 0
Out[8483]= 0
Out[8484]= 0.000809443998852
Out[8485]= -0.105278309271
During evaluation of In[8002]:= Py-Function - d^2/dxdz
Out[8487]= 0
Out[8488]= 0
Out[8489]= -0.000285685983732
Out[8490]= -0.0300688478423
During evaluation of In[8002]:= Pz-Function - d^2/dxdz
Out[8492]= 1.23228518779
Out[8493]= 0
Out[8494]= -0.000142844126695
Out[8495]= 0.0301435314879
During evaluation of In[8002]:= Dxx-Function - d^2/dxdz
Out[8497]= 0.00628387012988
Out[8498]= -0.00569272697329
Out[8499]= 0
Out[8500]= -0.000269390023198
During evaluation of In[8002]:= Dxy-Function - d^2/dxdz
Out[8502]= 0.000661243126546
Out[8503]= 0.0184515379115
Out[8504]= 0
Out[8505]= 0.000249962782101
During evaluation of In[8002]:= Dxz-Function - d^2/dxdz
Out[8507]= 0.000154290062861
Out[8508]= 0.0000258969272459
Out[8509]= -0.0444607527913
Out[8510]= 0.000187472112098
During evaluation of In[8002]:= Dyy-Function - d^2/dxdz
Out[8512]= -0.0411830327973
Out[8513]= -0.0177376150546
Out[8514]= 0
Out[8515]= -0.0000769685737339
During evaluation of In[8002]:= Dyz-Function - d^2/dxdz
Out[8517]= -0.000661243126546
Out[8518]= -0.0000361352473199
Out[8519]= -0.314599681951
Out[8520]= -0.0000999851233221
During evaluation of In[8002]:= Dzz-Function - d^2/dxdz
Out[8522]= 0.00628387012988
Out[8523]= 0.00282132825293
Out[8524]= 0
Out[8525]= -0.0000384843019961
During evaluation of In[8002]:= S-Function - d^2/dy^2
Out[8527]= -1.90486012083
Out[8528]= -0.589265770897
Out[8529]= 0.00354207798465
Out[8530]= -0.0000251994842742
During evaluation of In[8002]:= Px-Function - d^2/dy^2
Out[8532]= 1.23228518779
Out[8533]= 0
Out[8534]= -0.000142844126695
Out[8535]= -0.0149597402755
During evaluation of In[8002]:= Py-Function - d^2/dy^2
Out[8537]= 0
Out[8538]= 0
Out[8539]= -0.0000476131957936
Out[8540]= -0.0187370171672
During evaluation of In[8002]:= Pz-Function - d^2/dy^2
Out[8542]= 0
Out[8543]= 0
Out[8544]= 0.0000476147088982
Out[8545]= 0.00747987013777
During evaluation of In[8002]:= Dxx-Function - d^2/dy^2
Out[8547]= 0.0347464551086
Out[8548]= 0.0298158834476
Out[8549]= 0.289976751550
Out[8550]= 0.000115452875730
During evaluation of In[8002]:= Dxy-Function - d^2/dy^2
Out[8552]= -0.0947141661588
Out[8553]= -0.0400475913656
Out[8554]= 0.124182336147
Out[8555]= -0.0000333284093800
During evaluation of In[8002]:= Dxz-Function - d^2/dy^2
Out[8557]= -0.0601826256310
Out[8558]= -0.0258213125019
Out[8559]= 0
Out[8560]= -0.0000999851233221
During evaluation of In[8002]:= Dyy-Function - d^2/dy^2
Out[8562]= 0.0689968012890
Out[8563]= 0.0140425974665
Out[8564]= -0.246517482483
Out[8565]= -4.81051962604*10^-6
During evaluation of In[8002]:= Dyz-Function - d^2/dy^2
Out[8567]= 0.0947141661588
Out[8568]= 0.0200237956828
Out[8569]= 0
Out[8570]= 0.0000166642046900
During evaluation of In[8002]:= Dzz-Function - d^2/dy^2
Out[8572]= 0.0347464551086
Out[8573]= 0.00745397086190
Out[8574]= 0
Out[8575]= 0.0000288632189325
During evaluation of In[8002]:= S-Function - d^2/dydz
Out[8577]= 0
Out[8578]= 0
Out[8579]= 0.00707826429416
Out[8580]= -0.0000504001673179
During evaluation of In[8002]:= Px-Function - d^2/dydz
Out[8582]= 0
Out[8583]= 0
Out[8584]= -0.000285685983732
Out[8585]= -0.0300688478423
During evaluation of In[8002]:= Py-Function - d^2/dydz
Out[8587]= 0
Out[8588]= 0
Out[8589]= 0.0000476147088982
Out[8590]= 0.00747987013777
During evaluation of In[8002]:= Pz-Function - d^2/dydz
Out[8592]= 0
Out[8593]= 0
Out[8594]= 0.0000476147088982
Out[8595]= 0.00753588287198
During evaluation of In[8002]:= Dxx-Function - d^2/dydz
Out[8597]= 0.0164732131189
Out[8598]= 0.0141900920437
Out[8599]= 0
Out[8600]= 0.000153937147468
During evaluation of In[8002]:= Dxy-Function - d^2/dydz
Out[8602]= -0.0601826256310
Out[8603]= -0.0258213125019
Out[8604]= 0
Out[8605]= -0.0000999851233221
During evaluation of In[8002]:= Dxz-Function - d^2/dydz
Out[8607]= -0.000661243126546
Out[8608]= -0.0000361352473199
Out[8609]= -0.314599681951
Out[8610]= -0.0000999851233221
During evaluation of In[8002]:= Dyy-Function - d^2/dydz
Out[8612]= 0.0707746935497
Out[8613]= 0.0150978354912
Out[8614]= 0
Out[8615]= 0.0000192421509981
During evaluation of In[8002]:= Dyz-Function - d^2/dydz
Out[8617]= 0.00154290062861
Out[8618]= 0.0000421577885398
Out[8619]= -0.306627154417
Out[8620]= 0.0000374944271145
During evaluation of In[8002]:= Dzz-Function - d^2/dydz
Out[8622]= -0.0157096753247
Out[8623]= -0.00352666031616
Out[8624]= 0
Out[8625]= 0.0000192421509981
During evaluation of In[8002]:= S-Function - d^2/dz^2
Out[8627]= -1.90486012083
Out[8628]= -0.589265770897
Out[8629]= 0.00354207798465
Out[8630]= 0.0000504007667026
During evaluation of In[8002]:= Px-Function - d^2/dz^2
Out[8632]= 1.23228518779
Out[8633]= 0
Out[8634]= -0.000142844126695
Out[8635]= 0.0301435314879
During evaluation of In[8002]:= Py-Function - d^2/dz^2
Out[8637]= 0
Out[8638]= 0
Out[8639]= 0.0000476147088982
Out[8640]= 0.00753588287198
During evaluation of In[8002]:= Pz-Function - d^2/dz^2
Out[8642]= 0
Out[8643]= 0
Out[8644]= -0.0000476131957936
Out[8645]= 0.0149223984527
During evaluation of In[8002]:= Dxx-Function - d^2/dz^2
Out[8647]= 0.000152707558845
Out[8648]= 0.0000166901558139
Out[8649]= -0.121089474041
Out[8650]= 0.000115452875730
During evaluation of In[8002]:= Dxy-Function - d^2/dz^2
Out[8652]= -0.000661243126546
Out[8653]= -0.0000361352473199
Out[8654]= -0.314599681951
Out[8655]= -0.0000999851233221
During evaluation of In[8002]:= Dxz-Function - d^2/dz^2
Out[8657]= 0.0220324619169
Out[8658]= 0.00978782185672
Out[8659]= 0
Out[8660]= -0.0000333284093800
During evaluation of In[8002]:= Dyy-Function - d^2/dz^2
Out[8662]= 0.000954422242778
Out[8663]= 0.0000260783684592
Out[8664]= -0.272451316592
Out[8665]= 0.0000288632189325
During evaluation of In[8002]:= Dyz-Function - d^2/dz^2
Out[8667]= -0.0550811547921
Out[8668]= -0.0122347773209
Out[8669]= 0
Out[8670]= 0.0000166642046900
During evaluation of In[8002]:= Dzz-Function - d^2/dz^2
Out[8672]= -0.0128477041176
Out[8673]= -0.00282897790767
Out[8674]= 0.190840092897
Out[8675]= -4.81051962604*10^-6