Kink-Kink-Kink

quasi-6Pi-Kink - the velocities are 0.6 with difference equal 0.00001

> time_0:=time():

> Title:=`quasi-6Pi-Kink`:

> m:=3: # the total number of basic elements

> Digits:=32:

> v[1,1]:=1: v[1,2]:=0.59999: v[1,3]:=0: v[1,4]:=-10: # Kink:

> v[2,1]:=-1: v[2,2]:=0.6: v[2,3]:=0: v[2,4]:=-5: # Kink:

> v[3,1]:=1: v[3,2]:=0.60001: v[3,3]:=0: v[3,4]:=0: # Kink:

>

> B_p:=SGe_parameters(m,v):

> N:=B_p[0,0]; # the number of solitons

> Soliton:=(x,t)->SGe_backlund(x,t,N,B_p);

> Nt:=25: Nx:=25:

> L0:=-25: L1:=55:

> T0:=0: T1:=1.1*(L1-L0)/abs(v[2,2]):

> st:=time():

> A:=SGe_matrix(Soliton,L0,L1,T0,T1,Nx,Nt):

> matrixplot(A,axes=frame,labels=[t,x,``]);

> CPU_running_time:=time()-st;

> st:=time():

> SGe_animate(Soliton,Title,L0,L1,T0,T1,Nx,Nt);

> CPU_running_time:=time()-st;

> st:=time():

> Nt:=50: Na:=11:

> SGe_ribbon(Soliton,Title,L0,L1,T0,T1,Nx,Nt,Na);

> CPU_running_time:=time()-st;

> st:=time():

> Np:=10:

> SGe_pendulum(Soliton,Title,L0,L1,T0,T1,Nx,Nt,Np);

> CPU_running_time:=time()-st;

> Total_CPU_running_time_(in_seconds):=time()-time_0;

>