This function calculates the view factors of nRect rectangular areas.
function PhiRect
input Integer nRect "Number of rectangle areas";
input Real a1[nRect];
input Real a2[nRect];
input Real a3[nRect];
input Real b1[nRect];
input Real b2[nRect];
input Real b3[nRect];
input Real c1[nRect];
input Real c2[nRect];
input Real c3[nRect];
input Real d1[nRect];
input Real d2[nRect];
input Real d3[nRect];
input Integer r[nRect] "Normal direction of the rectangular areas";
input Integer z[nRect] "Number of segments";
input Boolean viewFacComp "true: numerical compensation of the view factors (only for closed hulls),
false: no numerical compensation of the view factors";
output Real PhiRect[nRect, nRect] "View factors for rectangular areas";
end PhiRect;