SetDirectory["Z:\\F2010-09-21"]

Length[Pd4]

13 16+4

Pd8=ReadList["scan_2010-09-22_010540",{Number,Number,Number,Number,Number,Number,Number,Number,Number,Number,Number,Number,Number,Number,Number,Number}];
(* resolucija 2'' na enoto, Fokus 10.50 *)

Pd7=ReadList["scan_2010-09-22_001302",{Number,Number,Number,Number,Number,Number,Number,Number,Number,Number,Number,Number,Number,Number,Number,Number}];
(* resolucija 2'' na enoto, Fokus 9.0 *)

Pd6=ReadList["scan_2010-09-21_235727",{Number,Number,Number,Number,Number,Number,Number,Number,Number,Number,Number,Number,Number,Number,Number,Number}];
(* resolucija 2'' na enoto, Fokus 12.10 *)

Pd5=ReadList["scan_2010-09-21_232853",{Number,Number,Number,Number,Number,Number,Number,Number,Number,Number,Number,Number,Number,Number,Number,Number}];
(* resolucija 2'' na enoto , Fokus 12.10*)

Pd4=ReadList["scan_2010-09-21_230918",{Number,Number,Number,Number,Number,Number,Number,Number,Number,Number,Number,Number,Number,Number,Number,Number}];
(* resolucija 2'' na enoto *)

Pd3=ReadList["scan_2010-09-21_225316",{Number,Number,Number,Number,Number,Number,Number,Number,Number,Number,Number,Number}];
(* resolucija 3'' na enoto *)

Pd2=ReadList["scan_2010-09-21_223934",{Number,Number,Number,Number,Number,Number,Number,Number,Number,Number,Number,Number}];
(* resolucija 3'' na enoto *)

Pd1=ReadList["scan_2010-09-21_222750",{Number,Number,Number,Number,Number,Number,Number,Number,Number,Number,Number,Number,Number,Number}];
(* resolucija 3'' na enoto *)

Pd=ReadList["scan_2010-09-21_215538",{Number,Number,Number,Number,Number,Number,Number,Number,Number,Number,Number,Number,Number,Number}];
(* resolucija 3'' na enoto *)

ListPlot3D[Pd,PlotRange->{0,6350},AxesLabel->{"x","y","z"}]

ListPlot3D[Pd2,PlotRange->{0,13000},AxesLabel->{"x","y","z"}]

ListPlot3D[Pd1,PlotRange->{0,13000},AxesLabel->{"x","y","z"}]

ListPlot3D[Pd3,PlotRange->{0,14000},AxesLabel->{"x","y","z"}]

ListPlot3D[Pd4,PlotRange->{0,5 14000},AxesLabel->{"x","y","z"}]

ListPlot3D[Pd5,PlotRange->{0,5 14000},AxesLabel->{"x","y","z"}]

Show[ListPlot3D[Pd6,PlotStyle->Opacity[0.9],PlotRange->{0,5 14000},AxesLabel->{"x","y","z"}]]

Show[ListPlot3D[Pd7,PlotStyle->Opacity[1],PlotRange->{0,5 14000},AxesLabel->{"x","y","z"}]]

ListPlot3D[Pd8,PlotRange->{0,5 14000},AxesLabel->{"x","y","z"}]

Max[Pd5]

ListContourPlot[Pd2,PlotRange->{0,Max[Pd2]},AxesLabel->{"x","y","z"}]

ListContourPlot[Pd3,PlotRange->{0,Max[Pd3]},AxesLabel->{"x","y","z"}]

ListContourPlot[Pd4,PlotRange->{0,Max[Pd4]},AxesLabel->{"x","y","z"}]

ListContourPlot[Pd5,PlotRange->{0,Max[Pd5]},AxesLabel->{"x","y","z"}]

ListContourPlot[Pd6,PlotRange->{0,Max[Pd6]},AxesLabel->{"x","y","z"}]

ListContourPlot[Pd7,PlotRange->{0,Max[Pd6]},AxesLabel->{"x","y","z"}]

Pq1=Table[{i,j,Pd1[[i,j]]},{i,Length[Pd1]},{j,Length[Pd1]}];

Kolaps[A_]:=Table[Sum[A[[i,j]],{j,Length[A]}],{i,Length[A]}]

KolapsY[A_]:=Table[Sum[A[[j,i]],{j,Length[A]}],{i,Length[A]}]

Show[ListPlot[Kolaps[Pd8]],PlotRange->All]

Fct[x_]=Fit[Kolaps[Pd3],{(E^(-(((x-5))^2)/2)),(E^(-(((x-6))^2)/2)),(E^(-(((x-7))^2)/2)),(E^(-(((x-8))^2)/2)),(E^(-(((x-9))^2)/2)),(E^(-(((x-10))^2)/2)),(E^(-(((x-11))^2)/2)),(E^(-(((x-12))^2)/2)),(E^(-(((x-13))^2)/2))},x]

Mx=FindMaximum[Fct[x],{x,10}]

Mx[[1]]

Solve[Fct[x]==Mx[[1]]/2,x]