Need help with SymEngine.jl?
Click the “chat” button below for chat support from the developer who created it, or find similar developers for support.

About the developer

154 Stars 40 Forks MIT License 384 Commits 35 Opened issues


Julia wrappers of SymEngine

Services available


Need anything else?

Contributors list


Build Status Build status Codecov Coveralls

SymEngine SymEngine

Julia Wrappers for SymEngine, a fast symbolic manipulation library, written in C++.


You can install

by giving the following command.
julia> Pkg.add("SymEngine")

Quick Start

Working with scalar variables

Defining variables

One can define variables in a few ways. The following three examples are equivalent.

Defining two symbolic variables with the names

, and assigning them to julia variables with the same name.
julia> a=symbols(:a); b=symbols(:b)

julia> a,b = symbols("a b") (a, b)

julia> @vars a b (a, b)

Simple expressions

We are going to define an expression using the variables from earlier:

julia> ex1 = a + 2(b+2)^2 + 2a + 3(a+1)
3*a + 3*(1 + a) + 2*(2 + b)^2

One can see that values are grouped, but no expansion is done.

Working with vector and matrix variables

Defining vectors of variables

A vector of variables can be defined using list comprehension and string interpolation.

julia> [symbols("α_$i") for i in 1:3]
3-element Array{SymEngine.Basic,1}:

Defining matrices of variables

Some times one might want to define a matrix of variables. One can use a matrix comprehension, and string interpolation to create a matrix of variables.

julia> W = [symbols("W_$i$j") for i in 1:3, j in 1:4]
3×4 Array{Basic,2}:
 W_11  W_12  W_13  W_14
 W_21  W_22  W_23  W_24
 W_31  W_32  W_33  W_34

Matrix-vector multiplication

Now using the matrix we can perform matrix operations:

julia> W*[1.0, 2.0, 3.0, 4.0]
3-element Array{Basic,1}:
 1.0*W_11 + 2.0*W_12 + 3.0*W_13 + 4.0*W_14
 1.0*W_21 + 2.0*W_22 + 3.0*W_23 + 4.0*W_24
 1.0*W_31 + 2.0*W_32 + 3.0*W_33 + 4.0*W_34



julia> expand(a + 2(b+2)^2 + 2a + 3(a+1))
11 + 6*a + 8*b + 2*b^2


Performs substitution.

julia> subs(a^2+(b-2)^2, b=>a)
a^2 + (-2 + a)^2

julia> subs(a^2+(b-2)^2, b=>2) a^2

julia> subs(a^2+(b-2)^2, a=>2) 4 + (-2 + b)^2

julia> subs(a^2+(b-2)^2, a^2=>2) 2 + (-2 + b)^2

julia> subs(a^2+(b-2)^2, a=>2, b=>3) 5


Peforms differentiation

julia> diff(a + 2(b+2)^2 + 2a + 3(a+1), b)
4*(2 + b)


is licensed under MIT open source license.

We use cookies. If you continue to browse the site, you agree to the use of cookies. For more information on our use of cookies please see our Privacy Policy.