例によって Rhino の API を C++ のプログラムで動かしてます．Excelを使って求めていた断面2次モーメントですが，移植して，CAD上だけで求めるようにしました．これでいろいろ実験できます．
というわけでまずは基本的な四角形から．
縦40mm，横40mmの正方形です．断面2次モーメントは
です．なので，
I=40*40*40*40/12=213333.3333....
となるはず．さぁ，どうでしょう？
若干誤差がありますね．被積分関数がそれなりに複雑なので，区分求積法で積分する際に誤差が出ています．まぁ，4桁信用できれば上出来なのかなぁ．．．．．
