最大の微分方程式

Question

私はmaximaで新しく、簡単な質問をすると本当にすみません。私は微分方程式を持って、

(%i1) -(x-x/2*sinh(x/2)+'diff(y,x))*(1/y+'diff(y,x)*x/y^2)+(x-x^2/sinh(x/2)+x^2*cosh(x/2)/(4*(sinh(x/2))^2)+'diff(y,x)*x+'diff(y,x,2)*x^2)/y+y^2-1-0.9*(x-x^2/(2*sinh(x/2)))=0; 
          2  x 
     2    2  x cosh(-) 
    2 d y  dy  x    2 
    x --- + x -- - ------- + ---------- + x 
     2  dx  x   2 x 
    dx   sinh(-) 4 sinh (-) 
         2   2 
(%o1) ---------------------------------------- 
        y 
         x   dy 
       x sinh(-)  x --       2 
      dy   2   dx 1  2    x 
    + (- -- + --------- - x) (---- + -) + y - 0.9 (x - ---------) - 1 = 0 
      dx  2    2 y       x 
           y      2 sinh(-) 
                   2 

(%i2) ode2(%,y,x); 
rat: replaced -0.9 by -9/10 = -0.9 
(%o2)        false 

私は何をすべきでしょうか？

Answer 1

あなたの持つ式は非線形です。 Maximaのode2は、限られた種類の微分方程式しか解くことができないため、方程式は扱えるカテゴリには含まれていないようです。

試してみることができるMaximaに別のシンボリック差分eqソルバーがあるかどうかわかりません。数値的な解法で十分なら、rk（Runge-Kuttaの実装）を見てください。