ある数値計算プログラムの中に次のような式が出てきました。

lambda = (re * re * re * re / 12.0 - re * ra * ra * ra / 3.0 + ra * ra * ra * ra / 4.0) / (re * re / 2.0 - re * ra + ra * ra / 2.0);

恐らく数値積分に関する何らかの公式から導かれた計算式だと思うのですが、
その公式が何なのか見当がつきません。台形公式とかルンゲクッタ法とか
そういうものを考えましたが、どうも違うようです。お気づきの点がありましたら
どんなことでも結構ですのでぜひご指導下さい。

なお、ra は以下のように定義されています。

ra = disa / 1.7724539; /* M_PI * ra^2 = disa^2 より */

また、re と disa には数値計算のパラメタとして定数が与えられています。

よろしくお願いいたします。