Aniž bychchtěl jakkoliv kritizovat práci pana Michla, je zcela zřejmé proč jsou výsledkyprá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ž bychviděl do kódu, který je zde ke stažení si dovoluji tvrdit, že pan Michl pouzeprohodil proměnné uvádějící vzdálenost těžiště k příslušnému okraji. Opravit tobude otázkou několika minut jeho práce.

Kromě tohobude celá utilitka (soudím podle délky kódu) pracovat pouze s obdélníkovýmprůřezem a tedy vzorci Ix=1/12 b h^3 => Wx=Ix/cy => pro obdélníkpak Wx = 1/6 b * h^2, respektive s údaji Autocadem odečtených hodnot Ix/Iy asouřadnic těžiště, a pak nelze použít dělení vypočteného Ix údajem o polozetěž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émto ale asi ani tak nevyřeší. Pokud se spokojíte s výpočtem pouze triviálníchpří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šittento problém obecně, pak vám nezbude než se pustit do integrování(Wx=Int(Ix/cy)dx)  a to vám paknezávidím, protože to si v LSP i s numerickými metodami asi užijete dostlegrace.

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