function [theta]=m2visAngle(D,S)
theta=(2*atan(S/(2*D)))*(180/pi);
end