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