function x=R(a)

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

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