wp=8.88*
13
10
;​​w=1.776*
14
10
;​​υ=1*
13
10
;​​wpv=
wp
υ
;​​wpw=
wp
w
;​​​​er=2.5;​​​​ao=1;​​wv=
w
υ
;​​​​mu=0.3;​​om=4*
12
10
;​​​​to=2.5;​​z=2*
-6
10
;​​v=3*
7
10
;​​t1=
τ
υ
-
z
v
;​​​​td=2*
-12
10
;​​​​u=
mu*om
υ
;​​R=30*
-7
10
;​​a=10*
-7
10
​​b=0.2745;​​b1=
2
a
2
R
*
2*1.098
1-
2
(a/R)
;​​l=Cos[om*t1];​​l1=0.86​​​​s=DSolveT'[τ]-ao^2*wv^2*
-1.5
T[τ]
*(1+l1)^2+u*
1-
2
l1
*(1+l1)*1-
2
wpw
*(1+b)(2*er)((wv*wpv-wpw*(1+b)/(2*er))^2+
-3
T[τ]
)0,T[0]0.1,T[τ],τ​​t1=Plot[Abs[T[τ]]/.s,{τ,1,4},PlotStyle{Thickness[0.006]},AxesLabel->{τ,T},PlotRangeAll]​​Show[t1]
1
1000000
0.86
Solve::tdep:The equations appear to involve the variables to be solved for in an essentially non-algebraic way.

DSolve::bvnul:For some branches of the general solution, the given boundary conditions lead to an empty solution.

{}
cos4000000000000-
1
10000000000000
+
τ
10000000000000

cos4000000000000-
1
10000000000000
+
τ
10000000000000
