title, 'test of rf-multipole';
option, debug, -echo;

bunchcharge := 2.e10;

twss1 : beta0, 
        betx=3.3999, bety=3.7645,
        alfx=-1.0879, alfy= 1.1907;

!!!
!!! it creates a sequence with an RF-Multipole and multipole that cancel each other
!!!

k:=sqrt(pi);
rfm1 : rfmultipole,  knl:={0, k, k, k }, ksl:={0, k, k, k, k}, pnl:={ 0, 0, 0,0,0,0,0,0,0 }, psl:={0,0,0,0,0,0,0,0}, freq=12, volt=100, lag=0.1;

test_line : line =( rfm1 );

!!! 
!!! TEST: together with test-rf-multipole-4.m (octave) it checks the symplecticity of the RF-Multipole 
!!!

!!! creeates the twiss table

select,flag=sectormap;
beam, sequence=test_line, bv=-1;
use, sequence=test_line;
twiss, betx=1,bety=1,x=1,px=1,t=1,rmatrix,sectormap,file='sectormap';