Hi lippie,
My post was aimed specifically at the very simple problem that vintage70 posted.
Essentially this only requires the combination of the section properties of two basic sections (The main I beam and the pair of angles), whose individual properties are known from standard section property tables, and which have a common axis of symmetry (so we know from the start that the principal axes are horizontal and vertical).
I cannot lay my hands on a text that sets out the method of calculation, so what follows is as good as I can offer. Please excuse the limitations that are inherent in these posts (I cannot, for example, set out the neat tabulation of values that I was taught to use back in the 1950's).
Let me start by defining my variables:
Datum line : I usually use a horizontal line through the bottom of the I beam as a datum, from which all vertical offsets

are taken. (Some prefer to use the NA of the main beam as the start datum - there are arguments in favour of either choice).
Variables:
Acomb = total area of combined section; Ybar = distance of NA of combined section above the assumed datum line;
AYcomb = 1st moment of area of the combined section about the assumed datum;
Icomb00 = "Moment of Inertia" (or 2nd moment of area if that is what you were taught to call it) I of the combined section
about the datum line;
Icombna = I of the combined section
about its own Neutral Axis;
ybtm = distance from combined NA to bottom fibre; ytop = distance from combined NA to top fibre.
A1=main beam area; y1=vertical distance to centroid of main beam; I1 = I of main beam about its own Neutral Axis;
A2=area of
pair of angles; y2=vertical distance to centroid of angles; I2 = I of
pair of angles about their Neutral Axis;
For easy legibility of any formal calculations I tabulate 6 columns for the following properties for every component section involved - Description of component section (with dimensions where appropriate),Area, y, Area*y, Area*y^2, I.
Particularly if you are calculating somewhat more complex sections (eg a PS concrete I beam, with triangular fillets between the web and both flanges) the tabulation keeps everything clear for anyone else looking at the calcs.
Basic equations:
Acomb = A1 + A2 = sum of my table column 2
AYcomb=A1*Y1+A2*Y2 = sum of table column 4
Icomb00 = A1*Y1^2 + I1 + A2*Y2^2 + I2 =sum of table columns 5 and 6
[for the more complex sections, these just become Acomb=A1+A2+A3+A4+..etc]
That is the hard work done

.
Now we get ybar from ybar=Aycomb/Acomb.
Icombna=Icomb00-Acomb*ybar^2
Here comes the benefit of taking the datum through the bottom of the main beam:
ybtm = ybar, Section modulus to bottom = Icombna/ybtm.
[with any other datum you have to calculate the bottom fibre distance as the distance from the datum to the bottom fibre of the main beam plus ybar]
ytop = Top fibre distance = overall depth - ybtm, Section modulus to top = Icombna/ytop.
All done. In fact it takes much longer to describe than to do on paper. Total calculations involved: 2 mental doublings (single angle area and I to give properties of the pair), 4 multiplications (Ay and Ay^2 ), 4 column additions (A,Ay,Ay^2,I) plus 4 simple expressions.