ttp://s-mac-p92.sap.hokkyodai.ac.jp/info/ex3/text/HTML/phys-ex1.html
すごい見にくいんですが，
これの４式を計算するために７式をn_Aで偏微分しているのですが
７式中にある密度dは（近似的に）一次関数的にn_Aに依存しているんです
(d = a n_A + bのように近似しています　これは確かです)
この微分がうまくいかないのか、最終的な４式で計算した値がおかしくなってしまうんです…
なんか数学のスレなのに本当に申し訳内です