In[]:=
F[L_?NumericQ]:=NIntegrate-
L
2
+x
3/2
2

L
2
+x
+
2
(y-y1)
+
2
(z-z1)
,y1,-
19
10000
,
19
10000
,z1,-
115
10^7
,
115
10^7
,x,-
6
1000
,
6
1000
,y,-
19
10000
,
19
10000
,z,-
125
10^6
,
125
10^6

In[]:=
{F[-20],F[1],F[10]}
Out[]=
{9.9636×
-18
10
,-3.98596×
-15
10
,-3.98545×
-17
10
}
In[]:=
Show[{ListLinePlot[Table[{L,F[L]},{L,-10,-10^-1,1/5}],PlotLabels{"F[L]"}],ListLinePlot[Table[{L,F[L]},{L,10^-1,10,1/5}],PlotLabels{"F[L]"}]},PlotRange->All,ImageSize->600]
Out[]=