b=0.1;
a=0.1;
γ=0.2;
β=0.4;
g=9.81;
α=0.3;
m=0.1;
θ=Pi;
ω=100;
κ=1;
v=0.5;
A=(((m+1)/(b(v+2a*κ)))(α*γ^2+ω(v+2a*κ)^2+κ^2(a+α)(v+2a*κ)^2))^(1/2m);
B=A^m*Sqrt[b/((m+1)(β-α))];
Animate[Plot[(A/Cosh[B(x-v*t)]^(1/m))Re[Exp[I(-κ*x+ω*t+θ)]],{x,0,4},PlotRange->{2,-2}],{t,0,5}]