function bearing = getbearingbetween(v1,v2)
        bearing = acos(dot(v1,v2)/(norm(v1)*norm(v2)))*180/pi;
        crossproduct = cross(v1,v2);
        dotproduct = dot(crossproduct, [0,1,0]);
        if dotproduct > 0
            bearing = -bearing +360;
        end