function [x,y,z] = rotation(x,y,z,az,el);
% function [x,y,z] = rotation(x,y,z,az,el)
% modifie les coodronnées en x,y et z d'un angle d'élévation et d'azimuth (en degrés)

% rotation selon axe y d'angle el (élévation, vertical)

a = el/180*pi;
roty = [cos(a) 0 -sin(a);0 1 0;sin(a) 0 cos(a)];

% rotation selon axe z d'angle az (azimut, horizontal)

a = az/180*pi;
rotz = [cos(a) -sin(a) 0;sin(a) cos(a) 0;0 0 1];

for i = 1:size(x,1)
    M = [x(i,:);y(i,:);z(i,:)];
    M = rotz * roty * M;
    x(i,:) = M(1,:);
    y(i,:) = M(2,:);
    z(i,:) = M(3,:);
end