変数係数の連立微分方程式

質問

下記の連立微分方程式を数値解析しようとしています。変数係数のため、どう解いていいのかわかりません。ルンゲクッタギル法などでも変数係数に対応したプログラムは見あたらないのですが、どなたかご存じでしょうか？

dA1(x)/dx=i*a*A1(x)+i*f(x)*A2(x)
dA2(x)/dx=i*q(x)*A1(x)+i*b*A2(x)

iは複素数、a,bは定数、f(x),q(x)が変数係数
A1(0)=1、A2(0)=0です。

回答No.1

Rossana · ベストアンサー

#define f(x) ■■■■
#define q(x) ■■■■
としておいて普通にルンゲクッタを適用してできるのでは？