function x=Q(a)

% returns the matrix for rotaion around the y-axis
%
% Initial version: Magnus Billvik, Swedish Institute of Space Physics,
%Uppsala, June 2005

x=[[cos(a) 0 -sin(a)];[0 1 0];[sin(a) 0 cos(a)]];
