1

これらの指示に従って、テイラー シリーズを生成します。

f(x) := ''(ratdisrep(taylor(qExct('x),'x,0,5)));

qExct は定義されていない関数です: 滑らかな関数である qExct に対して特定の計算を実行したいです。

xこれを知って、変数を特定の値 (例: 1)に設定するにはどうすればよいですか?

私がこれを行う場合：

f(1);

次に、マキシマは次のエラーを返します:

diff: variable must not be a number; found: 1

そして、私がそれを行うと:

f(D);

次に、変数を考慮Dし、すべての variablexを variableに置き換えますD。特に、d/dx の代わりに d/dD を使用して区別します。しかし、私が望むのは、変数をx^n 項のみ xで数値に置き換え、導関数をそのままにしておくことです…</p>1

これを行うにはどうすればよいですか?

4

2 に答える 2

1

atの代わりに使用する別のソリューションを次に示しますpdiff。

(%i1) f(x) := ''(ratdisrep(taylor(qExct('x),'x,0,2)));
                                !
                   2            !
               2  d             !
              x  (--- (qExct(x))!     )
                    2           !
                  dx            !                         !
                                !x = 0       d            !
(%o1) f(x) := ------------------------- + x (-- (qExct(x))!     ) + qExct(0)
                          2                  dx           !
                                                          !x = 0
(%i2) at(f(x), x=1);
                           !
                     !     !
        2            !     !
       d             !     !
       --- (qExct(x))!     !
         2           !     !
       dx            !     !                           !
                     !x = 0!                     !     !
                           !x = 1   d            !     !
(%o2)  -------------------------- + -- (qExct(x))!     !      + qExct(0)
                   2                dx           !     !
                                                 !x = 0!
                                                       !x = 1
(%i3) %, qExct=sin;
                                !
                          !     !
               2          !     !
              d           !     !
              --- (sin(x))!     !
                2         !     !
              dx          !     !                         !
                          !x = 0!                   !     !
                                !x = 1   d          !     !
(%o3)         ------------------------ + -- (sin(x))!     !
                         2               dx         !     !
                                                    !x = 0!
                                                          !x = 1
(%i4) %, nouns;
(%o4)                                  1

f(1)経由で評価されることに注意してくださいat(f(x), x=1)。

ネストされたat式は厄介です。もう起こらないように (Maxima のソース コードで) 修正しました。

于 2015-04-01T20:02:54.573 に答える
1

式の変数はMaximadiffのどこでもダミー (仮) 変数として認識されないため、 を評価しようとするf(1)と Maxima はdiff式に 1 を代入し、エラーを引き起こします。これはバグだと思います。それについてバグレポートを作成します。

回避策として、Maxima に含まれているアドオン パッケージpdiff(位置導関数) を使用できます。この表記は、Maxima でデフォルトで使用される dy/dx 表記とは少し異なります。

(%i1) load (pdiff) $
(%i2) f(x) := ''(ratdisrep(taylor(qExct('x),'x,0,2)));
                                    2
                     qExct     (0) x
                          ""(2)
(%o2)        f(x) := ---------------- + qExct     (0) x + qExct(0)
                            2                ""(1)
(%i3) f(h);
                                2
                   qExct   (0) h
                        (2)
(%o3)              -------------- + qExct   (0) h + qExct(0)
                         2               (1)
(%i4) ev (%, qExct=sin);
(%o4)                                  h
(%i5) ev (%o3, h=1);
                     qExct   (0)
                          (2)
(%o5)                ----------- + qExct   (0) + qExct(0)
                          2             (1)

""の表示のスプリアスf(x) := ...はマイナーな表示バグだと思います。無視してもいいと思います。

Maxima のインストールpdiffには、ドキュメントがあります。share/pdiff/pdiff-doc.pdf

于 2015-04-01T17:56:33.053 に答える