ニホンズイセン（日本水仙）

今日は仕事でした。

シンビジュームの花芽

一昨年末に買ったシンビジューム（品種名：マリーローランサン）の花芽がどんどん伸びています。
開花が楽しみです。

　　　　　　　

夕 日

日曜日，仕事決定。

体を動かさないと体がナマってしまうのですが・・・。

間もなくこの場所での夕日が撮れなくなります。

ユキワリイチゲ育つ

ユキワリイチゲが育っています。

むつこさんに頂いた白い花のユキワリイチゲ。


発泡スチロールに植えている水色の花が咲くユキワリイチゲ。
葉っぱの色が微妙に違うように見える。

無線LANカード

寒くなりましたね。

雪花が散ってました。

光通信無線ＬＡＮカードを取り付けました。

夕 景

二上山の見える風景

今朝は本当に久し振りに絵になる霧が出ました。

今日は田舎で農作業でした。

人参といちごの手入れ，一週間分の野菜を持って帰りました。

どこかでパーン，パーンと狩猟の音がしてました。

午前中は暖かく11時過ぎても霧がありました。
カメラは撮る緊張感が無い為かメモリーを忘れました。
なので携帯で。


午後は段々雲り空になって寒くなりました。

ナンテンの実

今日の天気予報は曇でしたが，雨でした。

今使ってるパソコン（中古）が３年程経ったので（新品時から７年）そろそろ次を考えなければと思い

今日は外付けハードディスクにデータを移動しました。

壊れてからでは遅いのでいつ壊れても良いように準備しておかなくてはいけません。

ＷＩＮ７では動かない古いソフトも使いたいのでこのマシンを温存したいのです。


写真は今年は野鳥がいないのでいつまでも赤い実をつけてるナンテンです。

ヒヨの鳴き声を一度も聞いて無いなー・・・。




　　　　　　　　

技術計算関数-ベッセル関数・Γ関数

日記として残しておきたくて書きます。

年末から正月にかけて勉強していたベッセル関数・Γ（ガンマ）関数。

地震時,水の入った水槽の側面・底面に作用する動水圧を求める

速度ポテンシャル法による式に用います。

まだ，奥が深くてこの式による答えが出ていません。

ネット検索でＶＢによるソースプログラムがありました。

これをＥＸＣＥＬのＶＢＡに書き換えたかったのですが途中になっています。

式とＶＢの抜粋です。

詳しくは ”ホームページ　ＷＡＮｔａｒｏHP”をご覧下さい。



　　　　　Pow:POWER関数　　　CDbl:CDBL関数　　FIはfloor関数でなく上記式Fi値　rbh=r/Hの記述は必要ないか？　

　　　　　ここまでの変数宣言等は省略
　　　　　　NN = 30
　　　　　　dz = 0.5
　　　　　　 z = -dz
　　　　　　Do
　　　　　　　z = z + dz 'サージタンク底面からの距離
　　　　　　　　sum = 0.0
　　　　　　　　　For i = 0 To NN
　　　　　　　　　　lambda = (2.0 * CDbl(i) + 1.0) / 2.0 * Math.PI
　　　　　　　　　　　sum = sum + Math.Pow(-1.0, CDbl(i)) / lambda * FI(r / H, lambda) * Math.Cos(lambda * z / H)
　　　　　　　　　Next i
　　　　　　　 '地震時動水圧強度(速度ポテンシャル)
　　　　　　　　pp = gamma * kh * r * Math.Cos(theta) * sum


　　　　　　　Private Function FI(ByVal rbh As Double, ByVal lambda As Double) As Double
　　　　　　　　Dim frbh As Double
　　　　　　　　frbh = 2.0 / rbh * FMBF(1, lambda * rbh) / (lambda * FMBF(0, lambda * rbh) - FMBF(1, lambda * rbh) / rbh)
　　　　　　　　Return frbh
　　　　　　　　End Function

　　　　　　　Private Function FMBF(ByVal nu As Integer, ByVal z As Double) As Double
　　　　　　　'第1種変形Bessel関数
　　　　　　　　Dim sum As Double
　　　　　　　　Dim i As Integer
　　　　　　　　sum = 0.0
　　　　　　　　　For i = 0 To NN
　　　　　　　　　　sum = sum + 1.0 / GAMMAF(CDbl(i + 1)) / GAMMAF(CDbl(nu + i + 1)) * Math.Pow(0.5 * z, CDbl(2 * i))
　　　　　　　　　Next i
　　　　　　　　Return Math.Pow(0.5 * z, CDbl(nu)) * sum
　　　　　　　　End Function

　　　　　　　　Private Function GAMMAF(ByVal z As Double) As Double
　　　 'ガンマ関数
　　　　　　　　Dim a As Double
　　　　　　　　　a = z / Math.E * Math.Sqrt(z * Math.Sinh(1.0 / z) + 1.0 / 810.0 / Math.Pow(z, 6.0))
　　　　　　　　Return Math.Sqrt(2.0 * Math.PI / z) * Math.Pow(a, z)
　　　　　　　　End Function

風 景

昼食を早めに済ませて昼一から夕方まで打合せ。
三日寝室で寝てませんが今日は普通に寝れます。
目を閉じて次の朝がすぐ来ると,又・・・。

今日は久し振りの雨。
携帯で撮ったタイミングの悪い窓越しの風景です。