function BR_test_int_box
%

MIN_Y = -6.25;
MAX_Y = 6.25;
MIN_X = -6.25;
MAX_X = 6.25;
MIN_Z = -6.25;
MAX_Z = 6.25;

faces = zeros(6,3);
faces(1:3,:) = [MIN_X,MIN_Y,MIN_Z;MAX_X,MIN_Y,MIN_Z;MIN_X,MAX_Y,MIN_Z];
faces(4:6,:) = [MIN_X,MIN_Y,MIN_Z;MAX_X,MIN_Y,MIN_Z;MIN_X,MIN_Y,MAX_Z];
faces(7:9,:) = [MAX_X,MIN_Y,MIN_Z;MAX_X,MAX_Y,MIN_Z;MAX_X,MIN_Y,MAX_Z];
faces(10:12,:) = [MIN_X,MAX_Y,MIN_Z;MAX_X,MAX_Y,MIN_Z;MAX_X,MAX_Y,MAX_Z];
faces(13:15,:) = [MIN_X,MIN_Y,MIN_Z;MAX_X,MAX_Y,MIN_Z;MIN_X,MIN_Y,MAX_Z];
faces(16:18,:) = [MIN_X,MIN_Y,MAX_Z;MAX_X,MIN_Y,MAX_Z;MIN_X,MAX_Y,MAX_Z];

normals = zeros(6,3);
normals(1,:) = [0,0,1];
normals(2,:) = [0,1,0];
normals(3,:) = [-1,0,0];
normals(4,:) = [0,-1,0];
normals(5,:) = [1,0,0];
normals(6,:) = [0,0,-1];

pos = [-6,0,-6];
p_new = [-3,0,-9];
[ip,dist2ip,normal] = BR_int_box(faces,normals,pos,p_new);

tch = 0;
