a=[1,2,3];
b=[0,-3,2];
myangle(a,b)

function angle=myangle( a,b )
angle=acosd(dot(a,b)/(norm(a)*norm(b)));
end

