olor
drawPlane3d(plane0, 'g');


%% Compute intersections

% compute intersection between 2 planes
line = intersectPlanes(plane0, plane);
drawLine3d(line, 'lineWidth', 2);

% compute intersection between a plane and a line
inter = intersectLinePlane(line12, plane0);
drawPoint3d(inter, 'marker', '+', 'markerSize', 10, 'linewidth', 3);

% create a line perpendicular to a plane
normal = planeNormal(plane);
perpLine = [80 60 40 normal];
drawLine3d(perpLine);

% compute intersection of line wi