In[]:=
m=0.05;B=0.007;sigma=0.007;l=0.08;​​k=
6*B
m*
2
l
;h=
m*l*
E
0
2*a
;
E
0
=
2
9
*sigma;a=
m*
2
l
3
;
n
1
=
B
a
2;
n
0
=0;
n
2
=2*
n
1
;
n
3
=
n
2
*2;​​
A
infinity
=l*h;​​
A
rez
=
2
3
*sigma;​​
p
*
=
2
k
Sqrt[
2
k
-2*
2
n
]
;​​
A
max
=
A
infinity
Sqrt
2
2
k
p
-1
+4*
2
n
1
p

;​​
A
n0
=
A
infinity
Sqrt
2
2
k
p
-1
+4*
2
n
0
p

;​​
A
n2
=
A
infinity
Sqrt
2
2
k
p
-1
+4*
2
n
2
p

;​​
A
n3
=
A
infinity
Sqrt
2
2
k
p
-1
+4*
2
n
3
p

;​​​​{maxA,pmax}={#1,p/.#2}&@@FindMaximum[Subscript[A,max],p];​​Plot{Subscript[A,max],Subscript[A,n0],Subscript[A,n2],Subscript[A,n3]},{p,0,550},PlotStyle{{Thick},{Dashed},{Dashed},{Dashed}},AxesLabelp[
-1
c
],A[м],PlotRange{0,0.007},PlotLegendsLineLegend[{
n
1
,
n
0
,
n
2
,
n
3
}],GridLines{{{pmax,Red},{k,Blue}},{{maxA,Red},{
A
rez
,Blue}}},GridLinesStyleDashed
Out[]=
32.8125
0
65.625
131.25