check if the position(s) are inside a hexagon

pos: coordinates of the positions, size (#positions, 2);
radius: radius or edge length of the hexagon
center: of the hexagon, by default is [0 0]
rotation: of the hexagon, by default is 0 degree
output: a bool vector of size (#positions, 1); true if inside.

by Yuanye Wang, Powerwave Technologies
yuanye.wang@pwav.com

An example is given below, which generates positions inside a hexagon with edge length 10, center position [3 10], and rotated 40 degree:

clear all; close all; clc
center = [3 10];
radius = 10;
rotation = 40;
pos = randn(50000,2) * radius + repmat(center, 50000,1);
x = inside_hex(pos, radius, center, rotation);
plot(pos(x,1),pos(x,2),'o')