In[]:=
ClearAll[a,β,ξ,h,ν,k,c,γ];​​​​w[ξ_]:=-
6
hγ
-2k-ν
-
2
6
c
β[ξ]
Y
γ
-2k-ν
/.
β[ξ]
Y
-
h
2c
-
-4cg+
2
h
Tanh
1
2
-4cg+
2
h
(ξ+ξ0)
2c
;​​​​"B1(x,t)="​​w[x_,t_]:=w[ξ]*Exp^(Iθ)//Simplify​​we=w[x,t]/.{ξμx+ωt,θτx+λt}​​c=1;g=-2.5;γ=4;μ=1;ω=5/3;ξ0=-1;h=2;λ=1;τ=1;​​​​Plot3D[Re[Sqrt[6]Exp[I(tλ+xτ)]Sqrt[-4cg+h^2]γTanh[1/2Sqrt[-4cg+h^2](xμ+ξ0+tω)]],{t,-2,2},{x,-2,2},PlotStyle{Red,Blue,Green},TicksStyleDirective[Blue,Thick,Italic],Mesh4,ColorFunction"Rainbow",AxesLabel{"n","t",None}]
Out[]=
B1(x,t)=
Out[]=
6
(t+x)
Exp
10.c+
2
h
γTanh
1
2
10.c+
2
h
-1+
5t
3
+x
-2k-ν
Out[]=