The hardware and bandwidth for this mirror is donated by METANET, the Webhosting and Full Service-Cloud Provider.
If you wish to report a bug, or if you are interested in having us mirror your free-software or open-source project, please feel free to contact us at mirror[@]metanet.ch.

The ‘giacR’ package

R interface to ‘Giac’

R-CMD-check

Giac is a general purpose symbolic algebra software. It powers the graphical interface Xcas. This package allows to execute Giac commands in R. You can find the documentation of Giac here.

Installation

remotes::install_github("rstudio/chromote")
remotes::install_github("stla/giacR")

Initialisation of a Giac session

The ‘chromote’ package is used to create a Giac session. If the find_chrome() function of ‘chromote’ returns NULL, you can set the path to the Chrome executable (or Chromium, Brave, etc) to the environment variable CHROMOTE_CHROME. Or you can pass it to the the Giac$new function. Since the Chrome executable is in my system path, I can use Sys.which("chrome").

library(giacR)
giac <- Giac$new(Sys.which("chrome"))

Examples

Elementary calculus

giac$execute("2 + 3/7")
## [1] "17/7=2.42857142857"

Gröbner basis

giac$execute("gbasis([x^3 - 2*x*y, x^2*y - 2*y^2 + x], [x, y])")
## [1] "[x^2,x*y,2*y^2-x]"

Antiderivative

giac$execute("integrate(ln(x))")
## [1] "x*ln(x)-x"

Infinite sum

giac$execute("sum(1/(n^2), n, 1, +(infinity))")
## [1] "1/6*pi^2=1.64493406685"

Exact rational roots of a polynomial

giac$execute("crationalroot(2*x^3 - 3*x^2 + 8*x - 12)")
## [1] "[2*i,3/2,-2*i]"

Solve a system of equations (and simplify the solutions)

giac$execute(
  "apply(simplify, solve([x^2+y+z=1, x+y^2+z=1, x+y+z^2=1], [x, y, z]))"
)
## [1] "list[[0,1,0],[1,0,0],[0,0,1],[sqrt(2)-1,sqrt(2)-1,sqrt(2)-1],[-sqrt(2)-1,-sqrt(2)-1,-sqrt(2)-1]]"

Determinant of a matrix with symbolic entries

giac$execute("det([[1, 2, 3], [3/4, a, b], [c, 4, 5]])")
## [1] "(-6*a*c+10*a+4*b*c-8*b+3)/2"

Check whether a variable occurs in an expression

giac$execute("has(x*y + u^2*z, u)")
## [1] "3"
giac$execute("has(x*y + u^2*z, w)")
## [1] "0"

Close session

giac$close()
## [1] TRUE

Blog posts

These binaries (installable software) and packages are in development.
They may not be fully stable and should be used with caution. We make no claims about them.