問題タブ [sympy]

大きな数の上昇階乗を計算する必要があります。これまでに見つけた最高のものは、sympy パッケージsympy packageの上昇階乗関数です。これは本当に素晴らしいですが、それでももっと速いものが必要です。

私が正確に必要とするのは、これの本当に高速なバージョンです:

編集: 上昇階乗 x(n) = x (x + 1)(x + 2)...(x + n-1) が与えられた場合、その展開された式の特定の乗数を取得したいと思います。

例えば：

与えられた: x(6) = x(x + 1)(x + 2)(x + 3)(x + 5)(x + 6)

展開された形式: x(6) = x**6 + 15*x**5 + 85*x**4 + 225*x**3 + 274*x**2 + 120*x

これらの乗数の 1 つを取得する方法が必要です (この場合は 1、15、85、225、274、120)。

「my_rising_factorial()」を使用すると、うまく機能します...しかし、本当に遅いです

Pythonでsympyを使用して以下の式を計算する必要がありますか？

a=6、、この場合b=5、c=2Pythonでsympyを使用して式を計算する方法は？私を助けてください。

私は、シンボル と が任意の型の未知数になる可能性があるシンボリック操作を許可する python ライブラリに興味があります。

これは私が書きたいコードです：

またはそれの大まかな概算。実際にはそれほど賢い必要はありません。上記のようなものを取得するために、多くの追加のイントロスペクション メソッドを呼び出さなければならないことを十分に嬉しく思います (たとえば、論理トートロジーが直接単純化されていなくても)。

私の最初の本能はsympyを見ることでしたが、ライブラリはシンボリック変数が数値でなければならないという強い仮定をしているようです。そして、少なくともシーケンスとセットを操作したいと思います:

非数値変数を理解するためにsympyを取得する方法、または同様のことができる別のライブラリはありますか?

sympy 0.7.1 が Jython で動作するかどうかは誰にもわかりませんか?

私は 3x3 行列を持っていますが、その逆数を計算します。一部の部分式が新しい記号に置き換えられた場合にのみ、逆は読みやすく記述できます。これは、それらが複数回出現するためです。sympy にこれらの部分式を見つけて置き換えてもらうことはできますか? 次のことを試しましたが、成功しませんでした：

sympy 0.6.7 を使用 (必要に応じてアップグレードできます)。

編集：

0.7.1-git (正確には cf9c01f8f9b4b749a7f59891f546646e4b38e580) にアップグレードし、実行しました (@PreludeAndFugue の提案に感謝します):

入手

では、なぜ「-ν_xy - 2⋅ν_xz² + 1」を m に置き換えないのでしょうか。

たとえば、私は解決したい

_{(ソース: texify.com )}

これが私が試したことです：

1.インタラクティブで、シンボリック変数も取れるモジュールで関数を定義したいと思います。

関数が

_{(ソース: texify.com )}
.

次に、次のように機能させたい

ユーザー入力を促す理由は、多くの引数を持つ関数を考えているからです。また、「f(x)+g(y)=h(z)」のような方程式での使用も考えているので、固定ではなく変数を代入できると便利です。

2.これが私が試したもので、うまくいきませんでした。その理由はケンダルの答えで説明されています。

初めてモジュールを作りました。

2 回 .py :

それから私はそれを実行しました、

私はこのコードを実行しました：

インタプリタは私に次のように与えます：「割り当て前に参照されるローカル変数x」

私が試しnsolve(f,x,(0,1))てみると、「与えられた許容範囲内でルートを見つけることができませんでした...」

（また、私はfindroot(f,(0,1))成功せずに試しました（インポートfrom mpmath import *してからmp.dps = 30; mp.pretty = True）。

私のプロジェクトでは、またはオブジェクトevaluate=falseを作成するときに使用する必要があります。この場合、これらのオブジェクトに等価チェックを適用すると問題に直面します。問題は、引数の順序によるものです。AddMul

以下の例を検討してください。

結果はfalseas k2.argsare(x,y,2)とk1.argsare(x,2,y)です。そのため、比較でタプルの等価性がチェックされている間、 が返されますfalse。必要な結果を得る方法はありますか?

また、タプルに何らかの操作を加えると (順序を逆にしてからチェックするなど)、 と が異なるオブジェクトから形成されている場合 ( =と=の場合k1など)に失敗します。k2Mulk1.args2*x,yk2.args2*y,x

Add([x+y,z],evaluate=False)この場合のように、ここでは並べ替えを使用できず、Add([x+z,y],evaluate=False)2 つの異なる式になります。またevalaute=True、この場合はAdd([x+y],x])andAdd([2*x+y])を使用すると同じになりますが、これは望ましくありません。

私はPython2で次のプログラムを作成して、数学の問題セットに対してニュートン法の計算を行いました。これは完全に機能しますが、私にはわからない理由で、最初にipythonにロードしたときに、Python3%run -i NewtonsMethodMultivariate.pyの分割がインポートされません。x**(3/4)Pythonプログラムをロードした後、 「1」と入力すると、これがわかります。新しい部門を手動でインポートした後、期待どおりにx**(3/4)残りx**(3/4)ます。どうしてこれなの？

この問題を回避するためにPython3を使用しますが、私のLinuxディストリビューションにはPython2用のSymPyしか出荷されていません。誰もが提供できる支援に感謝します。

また、誰かが疑問に思っている場合に備えて、私はまだnxn Jacobians用にこのスクリプトを一般化しておらず、問題セットで2x2を処理するだけで済みました。さらに、私のマシンにインストールされているSymPy 0.7.1は、「zeros（）は正確に1つの引数を取る」と文句を言うので、コマンドを使用する代わりに2x2のゼロ行列をスライスしています。たぶん、このコマンドはgitバージョン専用です。zeros(2,1)（ゼロ関数の問題を修正した表記を修正してくれたeryksunに感謝します。）