Aniž bych chtěl jakkoliv kritizovat práci pana Michla, je zcela zřejmé proč jsou výsledky právě takové, jaké jsou. Jak se říká v Cimrmanovi "... a jak jste na mě spěchali, tak v tom fofru jsem to zapomněl vynásobit dvěma...". Aniž bych viděl do kódu, který je zde ke stažení si dovoluji tvrdit, že pan Michl pouze prohodil proměnné uvádějící vzdálenost těžiště k příslušnému okraji. Opravit to bude otázkou několika minut jeho práce.

Kromě toho bude celá utilitka (soudím podle délky kódu) pracovat pouze s obdélníkovým průřezem a tedy vzorci Ix=1/12 b h^3 => Wx=Ix/cy => pro obdélník pak Wx = 1/6 b * h^2, respektive s údaji Autocadem odečtených hodnot Ix/Iy a souřadnic těžiště, a pak nelze použít dělení vypočteného Ix údajem o poloze těžiště pro jiný průřez než symetrický podle vyšetřovaných os, (kvůli hlavním osám průřezu, které už nemusí být totožné s těmi WCS nebo UCS).

Váš problém to ale asi ani tak nevyřeší. Pokud se spokojíte s výpočtem pouze triviálních případů (obdélník, T-průřez, trojúhelník, kruh, elipsa atd.) pak Vám stačí otevřít statické tabulky nalistovat vzorec a je po problému. Pokud chcete řešit tento problém obecně, pak vám nezbude než se pustit do integrování (Wx=Int(Ix/cy)dx)  a to vám pak nezávidím, protože to si v LSP i s numerickými metodami asi užijete dost legrace.

PS: pokud používám zastaralé značení veličin, tak mě prosím nekamenujte, nedostudoval jsem včera, ale před pár lety... Normy se mění, fyzikální zákony (naštěstí) a zvyky (bohužel) málokdy