% SH-Modell mit Crank Nicholson periodische Randbedingungen
nwave=20
xmax=nwave*2*pi
alpha=.6
R=.5
fig1=figure(1);

nx=200
dx=xmax/(nx+1)
%tmax=120
tmax=100
dt=0.1
nt=0;
ntplot=1;

x=[1:nx]*dx;

% 2 phase slips creating wavelengths
eps=4
q=.8
% no phase slips
eps=4
q=1
% 2 phase slips deleting wavelengths
eps=4
q=1.23

% random initial conditions
amp=.000001
y=amp*randn(1,nx);

x0=10
% insert the appropriate value for q
y=.4*sin(1.23*x+eps*sin((x-x0)/(2*nwave)));

 

