function slipsize(slip)
% print cycle-slip size (slip : L1,L2 slip cycle)
if nargin<1, slip=1:10; end

C=299792458; f1=1.57542E9; f2=1.22760E9;
lam1=C/f1; lam2=C/f2; lamw=1/(1/lam1-1/lam2); lamn=1/(1/lam1+1/lam2);
c1=f1^2/(f1^2-f2^2); c2=f2^2/(f1^2-f2^2);

ns=length(slip);
for n=1:ns, for m=1:ns, s1(n,m)=lam1*slip(n)-lam2*slip(m); end, end
for n=1:ns, for m=1:ns, s2(n,m)=lamw*(slip(n)-slip(m)); end, end
for n=1:ns, for m=1:ns, s3(n,m)=-(2*c2+1)*lam1*slip(n)+2*c2*lam2*slip(m); end, end
for n=1:ns, for m=1:ns, s4(n,m)=-2*c1*lam1*slip(n)+(2*c1-1)*lam2*slip(m); end, end

disp('---'),disp(' Cycle-Slip Size LG (cm)'),disp([' L1/L2 ',sprintf('%6.1f ',slip)])
for n=1:ns,disp([sprintf('%6.1f ',slip(n)),sprintf('%6.1f ',s1(n,:)*100)]); end

disp('---'),disp(' Cycle-Slip Size MW (cm)'),disp([' L1/L2 ',sprintf('%6.1f ',slip)])
for n=1:ns,disp([sprintf('%6.1f ',slip(n)),sprintf('%6.0f ',s2(n,:)*100)]); end

disp('---'),disp(' Cycle-Slip Size MP1 (cm)'),disp([' L1/L2 ',sprintf('%6.1f ',slip)])
for n=1:ns,disp([sprintf('%6.1f ',slip(n)),sprintf('%6.0f ',s3(n,:)*100)]); end

disp('---'),disp(' Cycle-Slip Size MP2 (cm)'),disp([' L1/L2 ',sprintf('%6.1f ',slip)])
for n=1:ns,disp([sprintf('%6.1f ',slip(n)),sprintf('%6.0f ',s4(n,:)*100)]); end
