function vcmbound(orig, r, thb, the, x, y, lm, M, t)
%VCMbound(orig, r, thb, the, x, y, lm, M, t)

xb = zeros((the-thb+1), 2); 
for theta = thb:the
	xt = orig + r * [cos(theta/180*pi), sin(theta/180*pi)];
	xb(theta-thb+1,:) = vcmline(orig, xt, x, y, lm, M, t);
end;
plot(xb(:,1), xb(:,2), 'k-');