sympy

0熱

2答えて

例えば、私は値aを持っています。 import sympy.physics.units as u a = 133.05*u.kg/u.m 私はfloat型の値に、この物理量を変換したい、と私は知っている唯一の方法は、float(a/(u.kg/u.m))です。

0熱

1答えて

SympifyError解析多項式

私は多項式関数を書いたPythonプログラムを持っており、Sympyのsolve()関数を使ってその逆を見つける。私は次のエラーを取得するプログラムの実行時に x = Symbol('x') y = Symbol('y') for i in range(10): p = Poly.create_random() print("p orig " + str(p))

1熱

1答えて

任意の3次元関数をグラフ化しようとするとエラーが発生する

私の多変数微積分クラスのプロジェクトに取り組んでいます。私の目的は、任意の関数f（x、y）をグラフ化し、等高線プロットを使って、/dx、df/dy）と関数の勾配のための震度プロットがありますが、より複雑な関数をグラフにしながら問題があります。 f（x、y）=（x + y）** 2のような関数入力の場合、プログラムはうまく動作しますが、より複雑な数学的概念を必要とする入力。X、Y）=罪（のx * y

2熱

1答えて

複素数を有するドット積

次のコード例は、複素数を持つ2つのベクトルの内積を実行します。 import sympy as sp from sympy import I sp.init_printing() b1 = sp.Rational(1, 2) * sp.Matrix([I, 1, 0, -I, 1]) v2 = sp.Matrix([2 * I, 1 + I, 0, 1 - I, 2 * I]) pr

2熱

1答えて

sympy plot_implicitを使用すると、どのように1つのウィンドウを表示できますか？

私はplot_implicit窓に関する質問があります。私はこれを入力し： from sympy import * from sympy.plotting import * r1, r2 = symbols('r1 r2') f = Function('f') f = 1/r1 - 1/r2 p0=plot_implicit(Eq(f,1),(r1,-1,1),(r2,-1,1))

1熱

2答えて

PythonのODEソルバーでエラーが発生しました

私はSpark Combustion Engineモデルで作業しています。私は燃焼をモデル化するためにPythonを使用している理由もあります。私はODEのソルバーを使用しようとしていますが、歩留まりは完全に現実のものです。私はシリンダーの体積の統合が間違っていることを発見した。私はすでに "odeint"と "ode"ソルバーを使ってみましたが、結果は同じです。このコードでは、Volumeをt

0熱

1答えて

テンソルフィールドのsympyでの発散

ベクトルフィールドの相違をsympyに入れることは可能です。たとえば、次のように >>> from sympy.physics.vector import ReferenceFrame >>> from sympy.physics.vector import divergence >>> R = ReferenceFrame('R') >>> vfield = 3*R[0]**2*R.x+2

0熱

1答えて

リポジトリ内のモジュールがPythonパッケージのローカルバージョンに見つかりません

physics.quantum.TensorProduct' method in sympy`を使用します。 the repoを見ると、この方法は間違いありません。私のpythonのセッションにインポートしようとすると、しかし、私は、次を得る： >>> from sympy import * >>> physics.quantum.TensorProduct(v1,v2) Traceback

3熱

2答えて

行列を要素ごとに掛ける

2つ以上の行列を持つsympyライブラリでmatrix_multiply_elementwiseから行うことはできますか？または、行列の2乗を要素ごとに乗算する他の方法はありますか？ P.S. numpyのでは、それは簡単ですが、私は、私はsympy

0熱

1答えて

シンボルをインデックス変数に置き換えるには？ XとH1におけるシンボルH0及びGを交換する方法

import sympy as s H0,z,g=s.symbols('H0 z g') def f(H0,g,z): H1=(6/H0)/(g+(6-g)*(1+z)**3) return H1 [0]、X [1]？（未置換することにより、私の実際のH1が大きい表現であるため）