function [angle]=myangle(a,b)
nda=size(a);
ndb=size(b);
if (nda(1)~=1 || ndb(1)~=1)
    disp('ERROR: array have dimension >1')
    return
end
if (length(a)~=3 ||length(b)~=3)
    disp ('ERROR: wrong length of array')
return
end
angle=acosd(dot(a,b)/(norm(a)*norm(b)));
return
end

