function S=rotation(S,delta_theta,M_rot,sigma_rot)

M=size(S(1,:));

for(j=1:M)
   E=rand_N(M_rot,(sigma_rot*delta_theta)^2,1,1);
   S(3,j)=S(3,j)+delta_theta+E; %theta
end