問題タブ [symbolic-computation]

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

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

ありがとうございました

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

出力：

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

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

出力：

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

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

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

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

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

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

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

編集

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

$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)$ のみ

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