ジュリアでのシンボリック演算の方法
概要
Juliaでの記号演算は、SymEngine.jl
1パッケージを通じて使うことができる。
コード
シンボルの定義
シンボルは、以下の方法で定義することができる。
julia> using SymEngine
julia> x = symbols(:x)
x
julia> x, y = symbols("x y")
(x, y)
julia> @vars x, y
(x, y)
julia> x = symbols(:x)
x
julia> f = 2x + x^2 + cos(x)
2*x + x^2 + cos(x)
ベクトルと行列
julia> v = [symbols("v_▷eq1◁i$j") for i in 1:2, j in 1:3]
2×3 Matrix{Basic}:
a_11 a_12 a_13
a_21 a_22 a_23
julia> A*v
2-element Vector{Basic}:
v_1*a_11 + v_2*a_12 + v_3*a_13
v_1*a_21 + v_2*a_22 + v_3*a_23
julia> @vars a, b, c, d, x, y
(a, b, c, d, x, y)
julia> [a b; c d] * [a x; b y]
2×2 Matrix{Basic}:
a^2 + b^2 a*x + b*y
a*c + b*d c*x + d*y
微分
記号微分は、Calculus.jl
パッケージでも使用することができる。
julia> f = 2x + x^2 + cos(x)
2*x + x^2 + cos(x)
julia> diff(f, x)
2 + 2*x - sin(x)
環境
- OS: Windows10
- バージョン: Julia v1.7.1, SymEngine v0.8.7