Maple/孤立子

亮孤立子

f1 := sqrt(2*(alpha^2-beta))*sech(sqrt(alpha^2-beta)*(x(t)-2*alpha*t)); with(plots)

> f2 := subs(alpha = 5, beta = 1, f1);

> animate(f2, x = -100 .. 100, t = -10 .. 10, numpoints = 200, frames = 50, color = red, thickness = 3);

暗孤立子

> f9 := 1-sech(x(t)-5*t)^2;

> animate(f9, x = -100 .. 100, t = -20 .. 20);

双孤立子

> f9a := sech(x(t)-5*t)^2+sech(-x(t)+6*t-12)^2;

> animate(f9a, x = -100 .. 100, t = -20 .. 20);

扭孤立子

> g1 := sech(lambda*x(t)+.3*lambda*y(t)+4*lambda*t)^2;

> g1s := subs(lambda = .8, g1);

> > animate3d(sech(x(t)^2+y(t)+100*t), x = -100 .. 100, y = -100 .. 100, t = -50 .. 50, color = cos(x*y));

3D Maple 图

animate3d(sech(x(t)^2+y(t)^2+.1*t), x = -100 .. 100, y = -100 .. 100, t = -50 .. 50, color = cos(x*y))

animate3d(cos(x(t)^2+y(t)^2+10*t)/cosh(x(t)^2+y(t)^2+5*t), x = -100 .. 100, y = -100 .. 100, t = -50 .. 50, color = cos(x*y))