”ロバートソン反応”の方程式を陰的ルンゲクッタで(Excel VBA)

ロバートソン反応を表す式

y1' = -0.04*y1 + 1.e4 * y2*y3

y2' = 0.04*y1 - 1.e4*y2*y3 - 3.e7*y2^2

y3' = 3.e7*y2^2

を陰的ルンゲクッタ（ラダウIIA)で計算。これも最初ステップ細かくして、あと粗くしないと大変。