Next: Number Theory, Previous: atensor [Contents][Index]
Next: Introduction to Series, Previous: Sums, Products, and Series, Up: Sums, Products, and Series [Contents][Index]
それぞれの和や積に独自のインデックスを与えることで、
式exprを変換します。
これは、
和や積と一緒に機能する時
changevarによりよい精度を与えます。
独自のインデックスの形式はjnumberです。
量numberは、gensumnumに参照することで決定されます。
これを、ユーザーは変更することができます。
例えば、gensumnum:0$はそれを再設定します。
Lの中のそれぞれの要素xに関するexprの和を表します。
もし引数Lがリストに評価されなければ、名詞形'lsumが返されます。
例:
(%i1) lsum (x^i, i, [1, 2, 7]);
7 2
(%o1) x + x + x
(%i2) lsum (i^2, i, rootsof (x^3 - 1, x));
====
\ 2
(%o2) > i
/
====
3
i in rootsof(x - 1, x)
和の外側の掛け算因子を内側に移動します。
もし外側の式でインデックスが使われているなら、
関数は、sumcontractに関してするのと同じように、合理的なインデックスを見つけようとします。
これは、本質的に、和のoutativeプロパティの逆の考えですが、
このプロパティを取り除かず、ただ無視するだけであることに注意してください。
いくつかの場合、
intosumの前に、scanmap (multthru, expr)が必要になるかもしれません。
インデックスiがi_0からi_1まで変えたexprの値の積を返します。
名詞形'productは、大文字Πとして表示されます。
productは、exprと下限上限i_0、i_1を評価し、
productは、インデックスiをクォートします（評価しません）。
もし上限と下限が整数差だけ違うなら、 exprは、インデックスiのそれぞれの値に関して評価され、 結果は陽な積です。
そうでなければ、インデックスの範囲は不定です。
積を整理するためにいくつかの規則が適用されます。
グローバル変数simpproductがtrueの時、
更なる規則が適用されます。
いくつかの場合、式整理は、積でない結果を出力します;
そうでなければ、結果は名詞形'productです。
nounsとevflagも参照してください。
例:
(%i1) product (x + i*(i+1)/2, i, 1, 4);
(%o1) (x + 1) (x + 3) (x + 6) (x + 10)
(%i2) product (i^2, i, 1, 7);
(%o2) 25401600
(%i3) product (a[i], i, 1, 7);
(%o3) a a a a a a a
1 2 3 4 5 6 7
(%i4) product (a(i), i, 1, 7);
(%o4) a(1) a(2) a(3) a(4) a(5) a(6) a(7)
(%i5) product (a(i), i, 1, n);
n
/===\
! !
(%o5) ! ! a(i)
! !
i = 1
(%i6) product (k, k, 1, n);
n
/===\
! !
(%o6) ! ! k
! !
k = 1
(%i7) product (k, k, 1, n), simpproduct;
(%o7) n!
(%i8) product (integrate (x^k, x, 0, 1), k, 1, n);
n
/===\
! ! 1
(%o8) ! ! -----
! ! k + 1
k = 1
(%i9) product (if k <= 5 then a^k else b^k, k, 1, 10);
15 40
(%o9) a b
デフォルト値: false
simpsumがtrueの時、
sumの結果は、整理されます。
この整理は、時々、閉形式を生成することができるかもしれません。
もしsimpsumがfalseもしくは、もしクォートされた形'sumが使われたなら、
値は、数学で使われるΣ表示の表現である和の名詞形です。
インデックスiがi_0からi_1まで変えたexprの値の和を返します。
名詞形'sumは、大文字Σとして表示されます。
sumは、被和exprと下限上限i_0、i_1を評価し、
sumは、インデックスiをクォートします（評価しません）。
もし上限と下限が整数差だけ違うなら、 被和exprは、インデックスiのそれぞれの値に関して評価され、 結果は陽な和です。
そうでなければ、インデックスの範囲は不定です。
積を整理するためにいくつかの規則が適用されます。
グローバル変数simpsumがtrueの時、
更なる規則が適用されます。
いくつかの場合、式整理は、和でない結果を出力します;
そうでなければ、結果は名詞形'sumです。
evflag(評価フラグ) cauchysumがtrueの時、
和の積は、コーシー積として表現されます。
コーシー積では、内側の和のインデックスは、独立に変化するのではなく、外側の和のインデックスの関数になります。
グローバル変数genindexは、
和の次のインデックスを生成するのに使われるアルファベット前置です。
gensumnumは、
自動生成されるインデックスが必要な時、
和の次のインデックスを生成するのに使われる数値接尾です。
gensumnumがfalseの時,
自動生成されるインデックスは、接尾なしのgenindexのみです。
sumcontract, intosum, bashindices, niceindices,
nouns, evflag, zeilbergerも参照してください。
例:
(%i1) sum (i^2, i, 1, 7);
(%o1) 140
(%i2) sum (a[i], i, 1, 7);
(%o2) a + a + a + a + a + a + a
7 6 5 4 3 2 1
(%i3) sum (a(i), i, 1, 7);
(%o3) a(7) + a(6) + a(5) + a(4) + a(3) + a(2) + a(1)
(%i4) sum (a(i), i, 1, n);
n
====
\
(%o4) > a(i)
/
====
i = 1
(%i5) sum (2^i + i^2, i, 0, n);
n
====
\ i 2
(%o5) > (2 + i )
/
====
i = 0
(%i6) sum (2^i + i^2, i, 0, n), simpsum;
3 2
n + 1 2 n + 3 n + n
(%o6) 2 + --------------- - 1
6
(%i7) sum (1/3^i, i, 1, inf);
inf
====
\ 1
(%o7) > --
/ i
==== 3
i = 1
(%i8) sum (1/3^i, i, 1, inf), simpsum;
1
(%o8) -
2
(%i9) sum (i^2, i, 1, 4) * sum (1/i^2, i, 1, inf);
inf
====
\ 1
(%o9) 30 > --
/ 2
==== i
i = 1
(%i10) sum (i^2, i, 1, 4) * sum (1/i^2, i, 1, inf), simpsum;
2
(%o10) 5 %pi
(%i11) sum (integrate (x^k, x, 0, 1), k, 1, n);
n
====
\ 1
(%o11) > -----
/ k + 1
====
k = 1
(%i12) sum (if k <= 5 then a^k else b^k, k, 1, 10);
10 9 8 7 6 5 4 3 2
(%o12) b + b + b + b + b + a + a + a + a + a
定数だけ異なる上限と下限を持つ足し算の和すべてを結合します。
結果は、
そんな和のそれぞれの集合が、この和を形成するよう抽出されなければならないすべての適切な余分の項に加えられた１つの和を含む式です。
sumcontractは、互換性のある和すべてを結合し、
可能なら和の１つからインデックスの１つを使い、もし供給されたどれもが使えないなら、合理的なインデックスを形成するよう試みます。
sumcontractの前に、intosum (expr)を実行する必要があるかもしれません。
デフォルト値: false
sumexpandがtrueの時、
和の積と、指数和は、入れ子の和に整理されます。
cauchysumも参照してください。
例:
(%i1) sumexpand: true$ (%i2) sum (f (i), i, 0, m) * sum (g (j), j, 0, n);
m n
==== ====
\ \
(%o2) > > f(i1) g(i2)
/ /
==== ====
i1 = 0 i2 = 0
(%i3) sum (f (i), i, 0, m)^2;
m m
==== ====
\ \
(%o3) > > f(i3) f(i4)
/ /
==== ====
i3 = 0 i4 = 0
Next: Functions and Variables for Series, Previous: Sums and Products, Up: Sums, Products, and Series [Contents][Index]
Maximaは、微分可能な関数の級数を見つけるために、
関数taylorとpowerseriesを含みます。
ある級数の閉形式を見つける性能があるnusumのようなツールも持ちます。
足し算や掛け算のような演算は、級数上で普通に機能します。
この節は、展開を制御するグローバル変数を提供します。
Next: Poisson series, Previous: Introduction to Series, Up: Sums, Products, and Series [Contents][Index]
デフォルト値: false
上限としてinfを持つ和同士を掛ける時、
もしsumexpandがtrue、かつ、cauchysumがtrueなら、
通常の積ではなくCauchy積が使われます。
Cauchy積では、
内側の和のインデックスは、独立に変化するのではなく、外側のインデックスの関数です。
例:
(%i1) sumexpand: false$
(%i2) cauchysum: false$
(%i3) s: sum (f(i), i, 0, inf) * sum (g(j), j, 0, inf);
inf inf
==== ====
\ \
(%o3) ( > f(i)) > g(j)
/ /
==== ====
i = 0 j = 0
(%i4) sumexpand: true$
(%i5) cauchysum: true$
(%i6) ''s;
inf i1
==== ====
\ \
(%o6) > > g(i1 - i2) f(i2)
/ /
==== ====
i1 = 0 i2 = 0
deftaylorは、
ある変数x_iの関数f_iそれぞれに関して、
expr_iをゼロの回りのTaylor級数と定義します。
expr_iは、典型的には、x_iの多項式か和です;
deftaylorは、もっと一般的な式も問題なく受け付けます。
powerseries (f_i(x_i), x_i, 0)は、
deftaylorで定義された級数を返します。
deftaylorは、
関数f_1, ..., f_nのリストを返します。
deftaylorは、引数を評価します。
例:
(%i1) deftaylor (f(x), x^2 + sum(x^i/(2^i*i!^2), i, 4, inf));
(%o1) [f]
(%i2) powerseries (f(x), x, 0);
inf
==== i1
\ x 2
(%o2) > -------- + x
/ i1 2
==== 2 i1!
i1 = 4
(%i3) taylor (exp (sqrt (f(x))), x, 0, 4);
2 3 4
x 3073 x 12817 x
(%o3)/T/ 1 + x + -- + ------- + -------- + . . .
2 18432 307200
デフォルト値: true
maxtayorderがtrueの時、
(切り詰められた)Taylor級数の代数操作の間、
taylorは、厳密とわかっているできるだけ多くの項を保とうとします。
exprの中の和や積のインデックスを改名します。
niceindicesは、
その名前が被加数や非積数の中に現れないなら、
インデックスそれぞれをniceindicespref[1]の値に改名しようとします。
現れた場合、
niceindicesは、
未使用の変数が見つかるまでniceindicesprefの次の要素を順に試します。
もしリスト全部が使い果たされたら、
例えば, i0, i1, i2, ....というように、
niceindicespref[1]の値に整数を追加することで、
追加のインデックスが構成されます。
niceindicesは式を返します。
niceindicesは引数を評価します。
例:
(%i1) niceindicespref;
(%o1) [i, j, k, l, m, n]
(%i2) product (sum (f (foo + i*j*bar), foo, 1, inf), bar, 1, inf);
inf inf
/===\ ====
! ! \
(%o2) ! ! > f(bar i j + foo)
! ! /
bar = 1 ====
foo = 1
(%i3) niceindices (%);
inf inf
/===\ ====
! ! \
(%o3) ! ! > f(i j l + k)
! ! /
l = 1 ====
k = 1
デフォルト値: [i, j, k, l, m, n]
niceindicesprefは、
niceindicesが和や積のインデックスの名前を取ってくる
リストです。
The elements of
niceindicesprefの要素は、
niceindicesによって強制されませんが、
通常、変数名です。
例:
(%i1) niceindicespref: [p, q, r, s, t, u]$
(%i2) product (sum (f (foo + i*j*bar), foo, 1, inf), bar, 1, inf);
inf inf
/===\ ====
! ! \
(%o2) ! ! > f(bar i j + foo)
! ! /
bar = 1 ====
foo = 1
(%i3) niceindices (%);
inf inf
/===\ ====
! ! \
(%o3) ! ! > f(i j q + p)
! ! /
q = 1 ====
p = 1
R.W. Gosperによる決定手続きを使って、 xに関するexprの不定超幾何総和を実行します。 exprと結果は、 整数べき、階乗、二項式、有理関数の積として表現可能でなければいけません。
用語「定」と「不定和」は、
「定」と「不定積分」へ類似して使われています。
不定に和を取ることは、
ただ、例えば0からinfまででなく、
変数の長さの区間上の和に関して、シンボリックな結果を与えることを意味します。
例えば、二項級数の一般的な部分和に関する公式はないので、
nusumはそれができません。
nusumとunsumは、有限積の和と差について少し知っています。
unsumも参照してください。
例:
(%i1) nusum (n*n!, n, 0, n);
Dependent equations eliminated: (1)
(%o1) (n + 1)! - 1
(%i2) nusum (n^4*4^n/binomial(2*n,n), n, 0, n);
4 3 2 n
2 (n + 1) (63 n + 112 n + 18 n - 22 n + 3) 4 2
(%o2) ------------------------------------------------ - ------
693 binomial(2 n, n) 3 11 7
(%i3) unsum (%, n);
4 n
n 4
(%o3) ----------------
binomial(2 n, n)
(%i4) unsum (prod (i^2, i, 1, n), n);
n - 1
/===\
! ! 2
(%o4) ( ! ! i ) (n - 1) (n + 1)
! !
i = 1
(%i5) nusum (%, n, 1, n);
Dependent equations eliminated: (2 3)
n
/===\
! ! 2
(%o5) ! ! i - 1
! !
i = 1
分子と分母の次数の和がべき級数の切り詰めレベル以下の 与えられたTaylor級数展開、すなわち、「最良」近似を持ち 加えて指定された次数範囲を満たす、有理関数すべてのリストを返します。
taylor_seriesは１変数Taylor級数です。 numer_deg_boundとdenom_deg_boundは、 分子と分母上の次数範囲を指定する 正の整数です。
taylor_seriesは
Laurent級数も可能です。
次数範囲は、infも可能で、
総次数が、冪級数の長さ以下の有理関数すべてを返すことになります。
総次数は
numer_deg_bound + denom_deg_boundとして定義されます。
べき級数の長さは
"truncation level" + 1 - min(0, "order of series")として定義されます。
(%i1) taylor (1 + x + x^2 + x^3, x, 0, 3);
2 3
(%o1)/T/ 1 + x + x + x + . . .
(%i2) pade (%, 1, 1);
1
(%o2) [- -----]
x - 1
(%i3) t: taylor(-(83787*x^10 - 45552*x^9 - 187296*x^8
+ 387072*x^7 + 86016*x^6 - 1507328*x^5
+ 1966080*x^4 + 4194304*x^3 - 25165824*x^2
+ 67108864*x - 134217728)
/134217728, x, 0, 10);
2 3 4 5 6 7
x 3 x x 15 x 23 x 21 x 189 x
(%o3)/T/ 1 - - + ---- - -- - ----- + ----- - ----- - ------
2 16 32 1024 2048 32768 65536
8 9 10
5853 x 2847 x 83787 x
+ ------- + ------- - --------- + . . .
4194304 8388608 134217728
(%i4) pade (t, 4, 4);
(%o4) []
このべき級数展開を持つ次数4の 分子/分母の有理関数はありません。 一般的に、 解くのに十分な数の未知の係数を持つために、 その和が少なくともべき級数の次数になるまで 分子の次数と分母の次数を増やさなければいけません。
(%i5) pade (t, 5, 5);
5 4 3
(%o5) [- (520256329 x - 96719020632 x - 489651410240 x
2
- 1619100813312 x - 2176885157888 x - 2386516803584)
5 4 3
/(47041365435 x + 381702613848 x + 1360678489152 x
2
+ 2856700692480 x + 3370143559680 x + 2386516803584)]
変数xに関する点a
(無限大のためにはinfかもしれません)
の回りのexprのべき級数展開の一般形式を返します:
inf
====
\ n
> b (x - a)
/ n
====
n = 0
もしpowerseriesがexprを展開することができないなら、
taylorが、級数の最初のいくつかの項を与えることができます。
verboseがtrueの時、
powerseriesは進捗メッセージを印字します。
(%i1) verbose: true$
(%i2) powerseries (log(sin(x)/x), x, 0);
can't expand
log(sin(x))
so we'll try again after applying the rule:
d
/ -- (sin(x))
[ dx
log(sin(x)) = i ----------- dx
] sin(x)
/
in the first simplification we have returned:
/
[
i cot(x) dx - log(x)
]
/
inf
==== i1 2 i1 2 i1
\ (- 1) 2 bern(2 i1) x
> ------------------------------
/ i1 (2 i1)!
====
i1 = 1
(%o2) -------------------------------------
2
デフォルト値: false
psexpandがtrueの時、
拡張有理関数展開が完全に展開されて表示されます。
スイッチratexpandは同じ効果を持ちます。
psexpandがfalseの時、
多変数式がちょうど有理関数パッケージにあるかのように表示されます。
psexpandがmultiの時、
変数に関する同じ総次数の項は一緒にまとめられます。
これらの関数は、
変数xに関するゼロの回りのTaylor級数exprの反転を返します。
revertは、
exprの最高次数と等しい次数の多項式を返します。
revert2は、次数nの多項式を返します。
nは、exprの次数よりも大きい値も小さい値も同じ値も取り得ます。
load ("revert")はこれらの関数をロードします。
例:
(%i1) load ("revert")$
(%i2) t: taylor (exp(x) - 1, x, 0, 6);
2 3 4 5 6
x x x x x
(%o2)/T/ x + -- + -- + -- + --- + --- + . . .
2 6 24 120 720
(%i3) revert (t, x);
6 5 4 3 2
10 x - 12 x + 15 x - 20 x + 30 x - 60 x
(%o3)/R/ - --------------------------------------------
60
(%i4) ratexpand (%);
6 5 4 3 2
x x x x x
(%o4) - -- + -- - -- + -- - -- + x
6 5 4 3 2
(%i5) taylor (log(x+1), x, 0, 6);
2 3 4 5 6
x x x x x
(%o5)/T/ x - -- + -- - -- + -- - -- + . . .
2 3 4 5 6
(%i6) ratsimp (revert (t, x) - taylor (log(x+1), x, 0, 6));
(%o6) 0
(%i7) revert2 (t, x, 4);
4 3 2
x x x
(%o7) - -- + -- - -- + x
4 3 2
taylor (expr, x, a, n)は、
式exprを、変数xのaの周りのTaylorもしくはLaurent級数を
(x - a)^nまで展開します。
もしexprが形式f(x)/g(x) の形であり、
g(x)がn次まで項を持たないなら、
taylorはg(x)を2 n次まで展開しようとします。
もしまだ0でない項がないなら、
taylorは、展開の次数がn 2^taylordepth以下である限り
g(x)の展開の次数を倍にしていきます。
taylor (expr, [x_1, x_2, ...], a, n)
は、すべての変数x_1, x_2, ...について
点(a, a, , ...)の周りでn次までのべき級数を返します。
taylor (expr, [x_1, a_1, n_1], [x_2, a_2, n_2], ...)は、変数x_1, x_2, ...について
点(a_1, a_2, ...)の回りで
n_1次, n_2次, ....まで展開したべき級数を返します。
taylor (expr, [x_1, x_2, ...], [a_1, a_2, ...], [n_1, n_2, ...])は、
変数x_1, x_2, ...について、
点(a_1, a_2, ...)の回りで
n_1次, n_2次, ....まで展開したべき級数を返します。
taylor (expr, [x, a, n, 'asymp])は、
exprのx - aの負のべき乗展開を返します。
最高次の項は(x - a)^-nです。
maxtaylorderがtrueの時、
（丸められた）Taylor級数の代数操作の間、
talyorは正確とわかっている限り多くの項を保とうとします。
psexpandがtrueの時、
拡張有理関数式は、フルに展開されて表示されます。
スイッチratexpandは同じ効果を持ちます。
psexpandがfalseの時、
有理関数パッケージかのように多変数式が表示されます。
psexpandがmultiなら、同じ総次数の項が一緒にグループ化されます。
展開を制御するには、taylor_logexpandスイッチも参照してください。
例:
(%i1) taylor (sqrt (sin(x) + a*x + 1), x, 0, 3);
2 2
(a + 1) x (a + 2 a + 1) x
(%o1)/T/ 1 + --------- - -----------------
2 8
3 2 3
(3 a + 9 a + 9 a - 1) x
+ -------------------------- + . . .
48
(%i2) %^2;
3
x
(%o2)/T/ 1 + (a + 1) x - -- + . . .
6
(%i3) taylor (sqrt (x + 1), x, 0, 5);
2 3 4 5
x x x 5 x 7 x
(%o3)/T/ 1 + - - -- + -- - ---- + ---- + . . .
2 8 16 128 256
(%i4) %^2;
(%o4)/T/ 1 + x + . . .
(%i5) product ((1 + x^i)^2.5, i, 1, inf)/(1 + x^2);
inf
/===\
! ! i 2.5
! ! (x + 1)
! !
i = 1
(%o5) -----------------
2
x + 1
(%i6) ev (taylor(%, x, 0, 3), keepfloat);
2 3
(%o6)/T/ 1 + 2.5 x + 3.375 x + 6.5625 x + . . .
(%i7) taylor (1/log (x + 1), x, 0, 3);
2 3
1 1 x x 19 x
(%o7)/T/ - + - - -- + -- - ----- + . . .
x 2 12 24 720
(%i8) taylor (cos(x) - sec(x), x, 0, 5);
4
2 x
(%o8)/T/ - x - -- + . . .
6
(%i9) taylor ((cos(x) - sec(x))^3, x, 0, 5);
(%o9)/T/ 0 + . . .
(%i10) taylor (1/(cos(x) - sec(x))^3, x, 0, 5);
2 4
1 1 11 347 6767 x 15377 x
(%o10)/T/ - -- + ---- + ------ - ----- - ------- - --------
6 4 2 15120 604800 7983360
x 2 x 120 x
+ . . .
(%i11) taylor (sqrt (1 - k^2*sin(x)^2), x, 0, 6);
2 2 4 2 4
k x (3 k - 4 k ) x
(%o11)/T/ 1 - ----- - ----------------
2 24
6 4 2 6
(45 k - 60 k + 16 k ) x
- -------------------------- + . . .
720
(%i12) taylor ((x + 1)^n, x, 0, 4);
2 2 3 2 3
(n - n) x (n - 3 n + 2 n) x
(%o12)/T/ 1 + n x + ----------- + --------------------
2 6
4 3 2 4
(n - 6 n + 11 n - 6 n) x
+ ---------------------------- + . . .
24
(%i13) taylor (sin (y + x), x, 0, 3, y, 0, 3);
3 2
y y
(%o13)/T/ y - -- + . . . + (1 - -- + . . .) x
6 2
3 2
y y 2 1 y 3
+ (- - + -- + . . .) x + (- - + -- + . . .) x + . . .
2 12 6 12
(%i14) taylor (sin (y + x), [x, y], 0, 3);
3 2 2 3
x + 3 y x + 3 y x + y
(%o14)/T/ y + x - ------------------------- + . . .
6
(%i15) taylor (1/sin (y + x), x, 0, 3, y, 0, 3);
1 y 1 1 1 2
(%o15)/T/ - + - + . . . + (- -- + - + . . .) x + (-- + . . .) x
y 6 2 6 3
y y
1 3
+ (- -- + . . .) x + . . .
4
y
(%i16) taylor (1/sin (y + x), [x, y], 0, 3);
3 2 2 3
1 x + y 7 x + 21 y x + 21 y x + 7 y
(%o16)/T/ ----- + ----- + ------------------------------- + . . .
x + y 6 360
デフォルト値: 3
もしまだ非ゼロ項がないなら、
展開の次数がn 2^taylordepth以下である限り、
taylorは、
g(x)の展開の次数を倍にします。
Taylor級数exprについての情報を返します。 戻り値はリストのリストです。 リストそれぞれは、変数名、展開点、展開次数から成ります。
もしexprがTaylor級数でないなら、
taylorinfoはfalseを返します。
例:
(%i1) taylor ((1 - y^2)/(1 - x), x, 0, 3, [y, a, inf]);
2 2
(%o1)/T/ - (y - a) - 2 a (y - a) + (1 - a )
2 2
+ (1 - a - 2 a (y - a) - (y - a) ) x
2 2 2
+ (1 - a - 2 a (y - a) - (y - a) ) x
2 2 3
+ (1 - a - 2 a (y - a) - (y - a) ) x + . . .
(%i2) taylorinfo(%);
(%o2) [[y, a, inf], [x, 0, 3]]
もしexprがTaylor級数なら、trueを、
そうでないなら、falseを返します。
デフォルト値: true
taylor_logexpandは、
taylor級数の中の対数の展開を制御します。
taylor_logexpandがtrueの時、
対数すべては完全に展開されるので、対数的恒等式を含むゼロ認識問題は
展開プロセスを邪魔しません。
しかしながら、分岐情報を無視するので、この方法はいつも数学的にただしいわけではありません。
taylor_logexpandがfalseに設定されている時、
生じる対数の唯一の展開は、
形式的なべき級数を得るのに必要なものです。
デフォルト値: true
taylor_order_coefficientsは、
Taylor級数の中の係数の順序付けを制御します。
taylor_order_coefficientsがtrueの時、
Taylor級数の係数は標準に順序付けられます。
べき級数exprの係数を整理します。
taylorはこの関数をコールします。
デフォルト値: true
taylor_truncate_polynomialsがtrueの時、
多項式は入力切り詰めレベルを基礎に切り詰められます。
そうでないなら、
taylorへの多項式入力は、
不定の精度を持つと考えられます。
taylor形式から標準有理式(CRE)形式に
exprを変換します。
効果はrat (ratdisrep (expr))と同じですが、より速いです。
一般式exprの内部表現をアノテートするので、まるでその和が切り詰められたTaylor級数かのように表示されます。 exprは別に変更されません。
例:
(%i1) expr: x^2 + x + 1;
2
(%o1) x + x + 1
(%i2) trunc (expr);
2
(%o2) 1 + x + x + . . .
(%i3) is (expr = trunc (expr));
(%o3) true
最初の後方差f(n) - f(n - 1)を返します。
従って、
unsumは、ある意味、sumの逆です。
nusumも参照してください。
例:
(%i1) g(p) := p*4^n/binomial(2*n,n);
n
p 4
(%o1) g(p) := ----------------
binomial(2 n, n)
(%i2) g(n^4);
4 n
n 4
(%o2) ----------------
binomial(2 n, n)
(%i3) nusum (%, n, 0, n);
4 3 2 n
2 (n + 1) (63 n + 112 n + 18 n - 22 n + 3) 4 2
(%o3) ------------------------------------------------ - ------
693 binomial(2 n, n) 3 11 7
(%i4) unsum (%, n);
4 n
n 4
(%o4) ----------------
binomial(2 n, n)
デフォルト値: false
verboseがtrueの時、
powerseriesは進捗メッセージを印字します。
Previous: Functions and Variables for Series, Up: Sums, Products, and Series [Contents][Index]
aをPoisson符号に変換します。
aをPoisson符号から一般表現に変換します。
もしaがPoisson形式でないなら、
outofpoisは変換を実行し、
すなわち、その戻り値は、outofpois (intopois (a))です。
例えば、この関数は、
特定のタイプのサインやコサイン項のべきの和に関する
標準整理器です。
aをbに関して微分します。 bは三角関数の引数の中だけ、または係数の中だけにいなければいけません。
関数的にintopois (a^b)と同一です。
bは正の整数でなければいけません。
(poisdiffと)似て制限された意味で積分します。
もしbが三角関数の引数の中にあるなら、
bの中の非周期的項は落とされます。
デフォルト値: 5
poislimは、三角関数の引数の中の係数の領域を決定します。
初期値5は
区間[-2^(5-1)+1,2^(5-1)]、すなわち[-15,16]に対応しますが、
[-2^(n-1)+1, 2^(n-1)]に設定することができます。
関数sinfnを与えられたPoisson級数のサイン項に、 cosfnをコサイン項に マップします。 sinfnと cosfnは、２引数関数です。 引数それぞれは、級数の中の項の係数と三角関数部です。
関数的にintopois (a + b)と同一です。
aを、一般表現のaに関するPoisson級数に変換します。
シンボル/P/は、Poisson級数式の行ラベルに続きます。
aをcの中のbに代入します。 cはPoisson級数です。
(1) bが変数u, v, w, x, y, zのいずれかの場合、
aはそれらの変数に関して線形の式(例えば、6*u + 4*v)でなければいけません。
(2) bはそれらの変数以外の場合、 aもまたそれらの変数を含んではいけなく、さらに、サインもコサインも含んではいけません。
poissubst (a, b, c, d, n)は、
上のタイプ(1)のようにaとbに関して演算しますが、
dがPoisson級数の場合、
cの中でbにa + dを代入した結果を供給するために、
cos(d)とsin(d)を次数nに展開する
特殊なタイプの代入です。
アイデアは、
dが小さなパラメータの項に関する展開だということです。
例えば、
poissubst (u, v, cos(v), %e, 3)はcos(u)*(1 - %e^2/2) - sin(u)*(%e - %e^3/6)をもたらします。
intopois (a*b)と同じ機能です。
(もしユーザーがそれを定義したら)
Poisson乗算の間、適用する予約関数です。
項の中のu, v, ..., zの係数を引数とする６引数の述語論理関数です。
(この項の係数に関して)poistrimがtrueの項は、
乗算の間に消去されます。
可読フォーマットでPoisson級数を印字します。
outofpoisと共通で、
もし必要なら、aを最初にPoisson符号に変換します。
Next: Number Theory, Previous: atensor [Contents][Index]