function x=P(a)

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

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