function [d1 d2] =drift_control(cycle,a)
        phi1 =  2*sin(0.1*cycle*a)*a;  
        phi2 = -2*sin(0.1*cycle*a)*a; 
        d1 = (rand.*phi1 + rand.*phi2);
        d2 = (rand.*phi1 + rand.*phi2);
end

