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;