function myangle( a,b )
%UNTITLED Summary of this function goes here
%   Detailed explanation goes here
ang=acosd(dot(a,b)/(norm(a)*norm(b)));
disp(['angle=',num2str(ang)]);
end

