自分のブログログアウトねむねこ幻想郷ネムネコ・ファミリーのために、大昔に書いたファンタジー小説を公開。数学のサイトではありません。前々回、メインの中に与えた境界積分サブルーティン(関数)は、以下でした。s)B1,B2,H1,H2の計算前に必要な積分パラメータを全て揃えます。そこには、B_Inte_?の中身として[内点方程式の積分公式(ラプラス型)]と[内点方程式から境界方程式へ(ラプラス型)]で導出した式を単純に並べたいという意図があります。積分パラメータを与える関数を検討します。最初にLength関数ですが、明らかに次でOKです。Function前々回に予告したように、関数とサブルーティンはローカル変数を持てる、外部手続きで統一します。Normal_LengthとNormal_Footの計算に疑問の余地はないのですが、ちょっと長そうなので後にします。次にBeta_kの扱いですが、前回の検討からGamm_1,Gamm_2の特殊ケースとして扱って良いとわかりました。いう訳でGamm_1とGamm_2の計算です。これらは単独に計算しては駄目で、Gamm_2-Gamm_1の形で計算する必要がありました。それをGam_12で表します。特異点(x0,y0)が要素節点(x1,y1)と(x2,y2)に一致しないケースから考え始めると、わかりやすいと思います。以下、前回を参照しつつ・・・。計算開始の前提として、図-1のr1,r2および(x0,y0),(x1,y1)と(x2,y2)は引数渡しされるとします。特異点と要素節点の一致/不一致判定のために、0距離判定値Epsも必要です。Ifの設定が基本です。Eps(r1=0)の時は、(x01,y01)=(1,0)とするのでした。