[angle]=myangle(a,b)
nda=length(size(a));
ndb=length(size(b));
if (nda>1 || ndb>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*(norma()*norm(b))
return
end

