

function myangle( a,b )
ang=acosd(dot(a,b)/(norm(a)*norm(b)));
disp(['angle=',num2str(ang)]);
end

