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))