9年前ある振動の方程式をC言語を用いて、ルンゲクッタで解析するプログラムについての質問です。式で表すと、x"+rx'+(k^2)sinx-((a^2)sin2(x-j))/4となる振動モデルをルンゲクッタで解析します。xはある一定の値に落ち着きます。いう条件でプログラムを実行すると、xは徐々に増加し1.489に落ち着きます。値を入れると出力xにマイナスの値が出力されます。j=10(度)のときは、xは-1.318...j=50(度)のときは、xは-0.650...が出力されてしまいます。ミクプラは自由にリンクできます。