Calculus in Finite Element Method
Chain rule

■簡単な関数をxで微分■
微分の作業をやっているとき、なにげなく使っているのが、Chain rule です。例えば、sin(αx) を x で微分することを考えてみましょう。結果は、αcos(αx) ですね。信じられない人は、この結果を下の差分で確認してみて下さい。

α=2, x=1, Δx=0.001 を上の式に代入すると、-0.83229が得られます。そして、αcos(αx)を計算しても-0.83229が得られます。

では、sin(αx) からどの様な過程を経て、αcos(αx) が得られたかを考えてみましょう。 まず、sin(αx) は、直接 x で微分出来ませんので、u=αx と置きます。すると、sin(u) となりますから、u で簡単に微分できますね。答えは、cos(u) です。また、u は、x で簡単に微分できます。言うまでも有りませんが、答えは、α ですね。すると、sin(αx) の x での微分は、αとcos(u) の積で得られます。まー、こんな風に教わったと思います。
上の事柄を、f=sin(αx), u=αx とおいて、まとめると、次の様になります。

■複雑な関数をxで微分■
上で説明した手順で、微分を行う方法を Chain rule と言います。もーちょと、複雑な関数を取り上げてみましょう。例えば、f=sin(eαx) を x で微分してみます。この場合、f=sin(u), u=ev, v=αx の様におきます。すると、df/dx は、次の様になる。

上の式の右辺を、よーく見ると、分子のduと分母のdu、そして、分子のdvと分母のdv をキャンセルできますよね。結果的に、左辺と同じになります。

上で説明した事柄を延長すると、Chain rule として、以下の様にまとめることができますね。ここでやっていることは、簡単な微分ですから、いつでも、差分で Chain rule の結果をチェックできます。

有限要素法では、この Chain rule を座標変換のところで、頻繁に使います。ですから、Chain rule の使い方をしっかり勉強しておいて下さいね。

Coordinate Transform L'hospital's rule
Menu Fnc Drv Int(1) Dif Crd Chain Lhsptl Vec Green IndxN Int(2)