~~~~~~~~~~~~~~
DR_24=atdrift('DR_24',2.4);
QD=atquadrupole('QD',0.2,-3.2243,'StrMPoleSymplectic4Pass');
DR_03=atdrift('DR_03',0.3);
QF = atquadrupole('QF',0.2,4.601,'StrMPoleSymplectic4Pass');
DR_031309 = atdrift('DR_031309',0.3131);
B = atrbend('B',1,0.7854,0,'BndMPoleSymplectic4E2Pass');
QDM = atquadrupole('QDM',0.2,-1.9590,'StrMPoleSymplectic4Pass');
SD = atsextupole('SD',0.1,-35.4586,'StrMPoleSymplectic4Pass');
DR_04 = atdrift('DR_04',0.4);
SF = atsextupole('SF',0.1,13.5023,'StrMPoleSymplectic4Pass');
QFM = atquadrupole('QFM',0.2,3.5793,'StrMPoleSymplectic4Pass');
M = atmarker('M','IdentityPass');
SF = atsextupole('SF',0.1,13.5023,'StrMPoleSymplectic4Pass');