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))