function [A] = aa(alpha)
%
%  This function computes the A-matrix for any value of alpha
%
radalpha = alpha*pi/180;
x = cos(radalpha);
A = [  0     1
      -1   -2*x ];
%
return
