07/01/11 18:27:01
ｄσｖ/ｄｘ +４(ｋσｖtanφ+Ｃｗ)/Ｄ＝γω＾２*（Ｒo+x）
dσｖ/ｄ（ｘ－Ｌ） +２（Ｋσｖtanα+２Ｃｗ）/{Ｄ-２（ｘ－Ｌ）tanα}＝γω＾２*（Ｒo+x）
解析的には解けないのでルンゲクッタ法で解きます。
初期条件は（１）がｘ＝０、σｖ＝σ₀、γ＝γ₀
（２）のσｖ、γはｘ＝Ｌにおける（１）の解
これのＥＸＣＥＬ ＶＢＡプログラムをお願いします。
正直時間がないのでどうか今週中にお願いできませんでしょうか。自分ではにっちもさっちも行きません。
この先に進めないと本当に危ないんです。