Other option : rSymPy tutorial

Ryacas Vignettes

\[\lim_{n\rightarrow\infty} (1-(1/n))^n\]

# install.packages('Ryacas')
library(Ryacas)


test <- yacas("Limit(n,Infinity)  (1+(1/n))^n")  ## n -> Inf
Eval(test)
## [1] 2.718282
print(test)
## expression(exp(1))

\[\lim_{h\rightarrow 0} \frac{\sin(x+h)-\sin(x)}{h}\]

test1 <- yacas("Limit(h,0) (Sin(x+h)-Sin(x))/h")
Eval(test1, list(x = pi))
## [1] -1
print(test1)
## expression(cos(x))