(*ポスト処理*)

(接点電位の表示)

Results=Point/@Table[Flatten[{Pts[[i]],v[i]*10}],{i,nPts}] //Simplify;(*グラフを適当な縮尺にするためにv[i]に係数をかける*)

Show[Graphics3D[{PointSize[0.05],Results} ],
ViewPoint->{2.980,1.502,0.558},Axes->True,FaceGrids->All]
[Graphics:indexgr2.gif][Graphics:indexgr30.gif]
[Graphics:indexgr2.gif][Graphics:indexgr31.gif]
[Graphics:indexgr2.gif][Graphics:indexgr32.gif]
[Graphics:indexgr2.gif][Graphics:indexgr33.gif]
[Graphics:indexgr2.gif][Graphics:indexgr34.gif]
[Graphics:indexgr2.gif][Graphics:indexgr35.gif]

(*電位を表示Sとして鳥瞰図にする*)

(*変数の準備*)

(*パラメータによる位置リストの作成*)

ParameterPosition=Table[
Flatten[{Apply[Plus,{ Pts[[ ele[[i]] ]] [[1]],
sP*tP*( Pts[[ ele[[i]] ]] [[2]] - Pts[[ ele[[i]] ]] [[1]] ),
sP*(1-tP)*( Pts[[ ele[[i]] ]] [[3]] - Pts[[ ele[[i]] ]] [[1]] )}] }],
{i,1,nElements}];

(*高速化のために重み係数を近似計算*)

kA=Rationalize[ N[k,7],10^-5 ];

ResultV=Table[
Flatten[{
1/10*kA[[i]].Flatten[Append[{1},ParameterPosition[[i]] ]]}],{i,1,nElements}];

ResultVSimp=Rationalize[ Simplify[ N[ResultV,5] ], 10^-5];
(*電位表示*)

ResultVC=Table[Hue[Flatten[ResultVSimp][[i]]/70],{i,1,nElements}]; //Simplify (*グラフ表示のために適当な係数をかける*)
[Graphics:indexgr2.gif][Graphics:indexgr36.gif]

ResultVO=Table[Flatten[{							
ParameterPosition[[i]],
ResultVSimp[[i]]*50,
ResultVC[[i]] }],
{i,1,nElements}];(*グラフ表示のために適当な係数をかける*)
[Graphics:indexgr2.gif][Graphics:indexgr37.gif]

ParametricPlot3D[Evaluate[ResultVO],{sP,0,1},{tP,0,1},Compiled->False,
PlotPoints->5,FaceGrids->All,ViewPoint->{2.446, 1.299, 1.945},
Lighting->False,Shading -> True, PlotRange->All]
Out of memory.  Exiting.

ParametricPlot3D[Evaluate[ResultVO],{sP,0,1},{tP,0,1},Compiled->False,
PlotPoints->10,FaceGrids->All,Shading -> True,
ViewPoint->{0.000,0.000,9.000},AxesLabel->{"X","Y","Z"},Lighting->False,PlotRange->All]
[Graphics:indexgr2.gif][Graphics:indexgr38.gif]
[Graphics:indexgr2.gif][Graphics:indexgr39.gif]

ParametricPlot3D[Evaluate[ResultVO],{sP,0,1},{tP,0,1},Compiled->False,
PlotPoints->10,FaceGrids->All,Lighting->False,
AxesLabel->{"X","Y","Z"},Lighting->False,PlotRange->All,ViewPoint->{-2.392, 2.081, 1.184}]
[Graphics:indexgr2.gif][Graphics:indexgr40.gif]
[Graphics:indexgr2.gif][Graphics:indexgr41.gif]
[Graphics:indexgr2.gif][Graphics:indexgr42.gif]
[Graphics:indexgr2.gif][Graphics:indexgr43.gif]
(*電界表示*)

ResultEC=Table[Hue[N[(Part[ kA[[i]],2]^2+Part[ kA[[i]],3]^2)^(1/2)]/8],{i,1,nElements}];  //Simplify
[Graphics:indexgr2.gif][Graphics:indexgr44.gif]
[Graphics:indexgr2.gif][Graphics:indexgr45.gif]
[Graphics:indexgr2.gif][Graphics:indexgr46.gif]

ResultEO=Table[Flatten[{							
ParameterPosition[[i]],
ResultVSimp[[i]]*50,
ResultEC[[i]] }],
{i,1,nElements}];(*グラフ表示のために適当な係数をかける*)
[Graphics:indexgr2.gif][Graphics:indexgr47.gif]

ParametricPlot3D[Evaluate[ResultEO],{sP,0,1},{tP,0,1},Compiled->False,
PlotPoints->6,FaceGrids->All,ViewPoint->{-1.199,2.909,1.245},
Lighting->False,PlotRange->All]
[Graphics:indexgr2.gif][Graphics:indexgr48.gif]
[Graphics:indexgr2.gif][Graphics:indexgr49.gif]

ParametricPlot3D[Evaluate[ResultEO],{sP,0,1},{tP,0,1},Compiled->False,
PlotPoints->6,FaceGrids->All,ViewPoint->{-1.312, 2.713, 1.540},
Lighting->False,PlotRange->All]
[Graphics:indexgr2.gif][Graphics:indexgr50.gif]
[Graphics:indexgr2.gif][Graphics:indexgr51.gif]

ParametricPlot3D[Evaluate[ResultEO],{sP,0,1},{tP,0,1},Compiled->False,
PlotPoints->6,FaceGrids->All,ViewPoint->{-0.017,9.000,-0.009},
Lighting->False,PlotRange->All]
[Graphics:indexgr2.gif][Graphics:indexgr52.gif]
[Graphics:indexgr2.gif][Graphics:indexgr53.gif]