%The function calculates the matrix for rotating the elastic stiffness matrix around y-axis.
%Inputs:
    %alpha- rotation angle.
%Outputs:
    %M- the matrix used to rotate elastic stiffness matrix around y-axis.
function M=BondY(alpha)
M=[cos(alpha)^2 0 sin(alpha)^2 0 -sin(2*alpha) 0;
   0 1 0 0 0 0;
   sin(alpha)^2 0 cos(alpha)^2 0 sin(2*alpha) 0;
   0 0 0 cos(alpha) 0 sin(alpha);
   sin(2*alpha)/2 0 -sin(2*alpha)/2 0 cos(2*alpha) 0;
   0 0 0 -sin(alpha) 0 cos(alpha)];