getDispersion

 function [dx,dy]=getDispersion(THERING,bpmindx)
 
 determines dispersion taking two orbits at plus and minus DE=0.0001

0001 function [dx,dy]=getDispersion(THERING,bpmindx)
0002 % function [dx,dy]=getDispersion(THERING,bpmindx)
0003 %
0004 % determines dispersion taking two orbits at plus and minus DE=0.0001
0005 %
0006 
0007 
0008 DE=0.001;
0009 
0010 Op=findorbit4(THERING,DE,bpmindx);
0011 Opx=Op(1,:);
0012 Opy=Op(3,:);
0013 
0014 Om=findorbit4(THERING,-DE,bpmindx);
0015 Omx=Om(1,:);
0016 Omy=Om(3,:);
0017 
0018 dx=(Opx-Omx)./(2*DE);
0019 dy=(Opy-Omy)./(2*DE);
0020 
0021 return