問題タブ [symbolic-computation]

For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.

0 投票する
2 に答える
2462 参照

r - Rの行列の逆行列

Rの行列のシンボリック逆行列を見つけるにはどうすればよいですか。例えば:

「yacas」へのインターフェイスである「Ryacas」というパッケージがあることは知っていますが、このような計算を行うためにそれを適用することはできませんでした。「yacas」は、数式の記号操作のためのプログラムです。詳細については、リンクを参照してください。

ありがとうございました

0 投票する
2 に答える
494 参照

f# - C# ライブラリは ^ 演算子をオーバーロードします。代わりに ** を使用するには?

Symbolismライブラリは、算術演算子をオーバーロードします。C# で書かれていますが、F# から使用できます。

出力:

しかし、それ^は力の過負荷にもなります。もちろん、これは F# ではうまく機能しません。

回避策として、powers のメソッド グループをエクスポートしました。

出力:

代わりにメソッド グループ**を使用するにはどうすればオーバーロードできますか?Aux.Pow

私はこのようなことができます:

そして、それはMathObject値に対して機能します:

しかし、同様Aux.Powにオーバーロードさintれています:

どんな提案でも大歓迎です!

0 投票する
1 に答える
4699 参照

python - sympy integrate が返す区分方程式の扱い

sympy には、Piecewise オブジェクトを返す積分があります。

この区分方程式の最初の枝だけを抽出したい、つまりitest.parts(0)単純に を抽出するようなことができるようにしたいのですsqrt(pi)*z/2。これを行う方法が見つからないようですが、ドキュメントで間違った検索用語を使用している可能性があります。何か案は?

編集

itest.args[0][0]少し調べてみると、この式を抽出できることがわかりました。ただし、これはちょっとしたハックのようです。より良いアプローチはありますか?

0 投票する
0 に答える
809 参照

python - 大きな三角関数の統合

SymPy で、以下に示すような大きな三角関数を統合するための推奨される方法は何ですか? 数値積分により が得られ-27151.225ます。

このタイプの関数は、ここで詳述されている問題で生成されます。

0 投票する
1 に答える
159 参照

python - Sympy、クラスを強制的に式の左側または右側に留まらせることは可能ですか?

この回答では、などargsの SymPy クラスでソートされないようにする方法が説明されています。MulAdd

ここで説明されているような新しく作成されたクラスの場合、 sympy.core.numbers.Floatsympy.core.numbers.Piまたはを掛けるsympy.core.numbers.Integerと、たとえば、次のようにすると右側に移動できます。

元の式は0.評価されると を返しますが、新しい式は を返しますD(x)

この微分演算子が適切に機能するためには、左側にとどまる必要があります。

_op_priority優先度を変更する場合のように__mul__()、SymPy に左または右に留まるタイプを伝える隠しパラメーターはありますか?

0 投票する
4 に答える
686 参照

python - 前の用語 sympy を使用してシーケンスを生成する

前に生成された要素が次の要素に含まれるようなシーケンスを生成したいのですが、これを行う方法がわかりません。

つまり、その項目が次のようになるようなリストを生成します。

ここで、x は単なる Sympy シンボルです

[x,(x)*(x+1),(x)*(x+1)*(x+2)]

それよりも[x,x+1,x+2]

みたいなことを考えている

k.append(k*(K+o))

しかし、型エラーが発生し続けます

どんな助けでも大歓迎です!

0 投票する
2 に答える
601 参照

python - cos(x)*cos(2x)*...*cos(mx) を SAGE 経由で統合する

$I_m=\int_0^{2\pi} \prod_{k=1}^m cos(kx){}dx$ を検索します。ここで、$m=1,2,3\ldots$ です。

シンプルな SAGE コード:

出力:

ご覧のとおり、数値の答えは正しいのですが、integrate(...) の結果は $m=1,2,\ldots,7$ に対して正しく、バグがあります。不定積分を出力できます。

そして出力:

したがって、$m=7$ の場合、WolframAlpha を介した Indef_I_7と比較して正しい回答が得られ、WolframAlpha を介したIndef_I_8の場合 、回答が正しく ありません Indef_I_8=$\frac{7x}{128}+\ldots$ が必要で、$\ はありませんsin(x)$、合計で $\sin(3x)$、$k=1,2,3,\ldots 18$ の場合は $\sin(2k)$ のみ

体積計算でごめんなさい!質問は、シンボリック統合のバグであるというのは正しいですか?