Wolfram Cloud Document

In[]:=

L=10^-2;s=0.2*10^-6;m=10^-6;he=0.02;distribution=

L

s*

2*Pi

Exp

-L^2(r-m/L)^2

2*s^2

;Plot[distribution,{r,0,2*10^-4},PlotRangeAll,AxesOrigin{0,-1000}]

Out[]=

In[]:=

q=0.02*p1*r^4-

8

3

r^3+

16

3

1

p1^3

*distribution;PlotIntegrateq,r,

2

p1

,Infinity,{p1,1,10^5},AxesOrigin{0,-10^-14},PlotRangeAll

Out[]=

Forumula:(p1^6/(128he))5(

2

δ

flow)

2

δp1

+p1(

3

δ

flow)

3

δp1



d1=DIntegrateq,r,

2

p1

,Infinity,{p1,2};d2=DIntegrateq,r,

2

p1

,Infinity,{p1,3};ListPlotTable

(Range[1,10^5,100])^6

128*he

((5*Table[d1,{p1,1,10^5,100}])+(Range[1,10^5,100]*Table[d2,{p1,1,10^5,100}])),{r,0,2*10^-4},PlotRangeAll,AxesOrigin{0,-1000}

General

:Exp-4.9995×

9

10

 is too small to represent as a normalized machine number; precision may be lost.

General

:Exp-4.9995×

9

10

 is too small to represent as a normalized machine number; precision may be lost.

General

:Exp-4.9995×

9

10

 is too small to represent as a normalized machine number; precision may be lost.

General

:Further output of General::munfl will be suppressed during this calculation.

Out[]=