Moving Kink and Standing Breather
> time_0:=time():
> Title:=`Moving Kink and Standing Breather`:
> m:=2: # the total number of basic elements
> v[1,1]:=1: v[1,2]:=0.1: v[1,3]:=0: v[1,4]:=-10: # Kink:
> v[2,1]:=1: v[2,2]:=0: v[2,3]:=0.1: v[2,4]:=0: # Breather:
> 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:=35: Nx:=35:
> L0:=-22: L1:=12:
> T0:=0: T1:=0.85*(L1-v[1,4])/(v[1,2]):
> st1:=time():
> A:=SGe_matrix(Soliton,L0,L1,T0,T1,Nx,Nt):
> matrixplot(A,axes=frame,labels=[t,x,``]);
> time()-st1;
> st:=time():
> SGe_animate(Soliton,Title,L0,L1,T0,T1,Nx,Nt);
> CPU_running_time:=time()-st:
> st:=time():
> Na:=10:
> 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);
> time()-st;
> Total_CPU_running_time_(in_seconds):=time()-time_0;