similar to: Derivative of a Function Expression

Hi, Could anyone tell me what is the command used in R to do 1. Differentiation 2. Newton Raphson method (Numerical Analysis in general...) Are there any packages separately for this? Thanks for your help! BR, Shubha [[alternative HTML version deleted]]

Hi (If you're wondering, this is a Project Euler question :)) If I wanted to calculate the sum of the digits in the decimal representation of 2^1000, what would be a good way to go about that? I've tried the following methods: # Calculate the sum of digits in the decimal representation of 2^n # Only works for smaller values of n bsum <- function(n) { s <- 0 e <-

Hi all, I am trying to implement the following matlab code with Ryacas : syms U x x0 C d1=diff(U/(1+exp(-(x-x0)/C)),x); pretty(d1) d2=diff(U/(1+exp(-(x-x0)/C)),x,2); pretty(d2) solx2 = solve(d2 == 0, x, 'Real', true) pretty(solx2) slope2=subs(d1,solx2) I have tried the following : library(Ryacas) x <- Sym("x");U <- Sym("U");x0 <-

Thanks for the response. Yes, I did study the vignette but did not understand it fully. Anyway, I have tried once again now. I am happy to say that I have got what I wanted. library(Ryacas) x <- Sym("x");U <- Sym("U");x0 <- Sym("x0");C <- Sym("C") my_func <- function(x,U,x0,C) { return (U/(1+exp(-(x-x0)/C)))} FirstDeriv <-

Have you studied the "Introduction to Ryacas" vignette that come with the package? Cheers, Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Tue, Sep 19, 2017 at 2:37 AM, Vivek Sutradhara <viveksutra at gmail.com> wrote:

Dear all, I want to equalize a symbolic derivative (of a function with two variables) with zero and convert it to a variable, e.g. x. I'm computing the derivative by: (found it in the archive) library(Ryacas) f <- function(x,y) (100-x-y)*x-10*x yacas(f) # register f with yacas Df <- f body(Df) <- yacas(expression(deriv(f(x,y),x)))[[1]] Df R-Output: function (x, y) 100 - x - y -

Hello, I am trying to construct two matrices, F and V, composed of partial derivatives and then find the eigenvalues of F*Inverse(V). I have the following equations in ryacas notation: > library(Ryacas) > FIh <- Expr("betah*Sh*Iv") > FIv <- Expr("betav*Sv*Ih") > VIh <- Expr("(muh + gamma)*Ih") > VIv <- Expr("muv*Iv") I

I added "plogis" to the derivative table in the development version of R; the patch against yesterday's R-devel src/deriv/main.c is available at http://www.zoology.ufl.edu/bolker/deriv_patch.txt . I pretty much followed the framework of the other symbols; here was my incantation - } else if (CAR(expr) == PlogisSymbol) { - ans = simplify(TimesSymbol, -

Ryacas is an R interface to the free yacas computer algebra system. Ryacas allows one to send R expressions, unprocessed yacas strings and certain other R objects to a separate yacas process from R and get back the result. It also has facilities for manipulating yacas strings and R expressions destined for yacas processing. It can be used for exact arithmetic, symbolic math, ASCII pretty

Ryacas is an R interface to the free yacas computer algebra system. Ryacas allows one to send R expressions, unprocessed yacas strings and certain other R objects to a separate yacas process from R and get back the result. It also has facilities for manipulating yacas strings and R expressions destined for yacas processing. It can be used for exact arithmetic, symbolic math, ASCII pretty

Ryacas is an R interface to the free yacas computer algebra system. Ryacas allows one to send R expressions, unprocessed yacas strings and certain other R objects to a separate yacas process from R and get back the result. It also has facilities for manipulating yacas strings and R expressions destined for yacas processing. It can be used for exact arithmetic, symbolic math, ASCII pretty

Hello, I am wondering if someone can help me. I have the following function that I derived using nls() SSlogis. I would like to find its derivative. I thought I had done this using deriv(), but for some reason this isn't working out for me. Here is the function: asym <- 84.951 xmid <- 66.90742 scal <- -6.3 x.seq <- seq(1, 153,, 153) nls.fn <- asym/((1+exp((xmid-x.seq)/scal)))

Dear All I am doing the following: > x <- yacas("3/2") > for (i in 2:400) + x <- yacas(paste(x,"*",x)) > x expression(Inf^1.260864167e+117/Inf^6.304320836e+116) > Eval(x) [1] NaN No luck this way. However, I am successful with > y <- yacas("(3/2)^400") > y expression(7.05507910865533e+190/2.58224987808691e+120) > Eval(y) [1]

Hello yet again, R People: I was working with Ryacas and yacas last night and all was well. Now this morning, I keep getting the following: > a <- Sym("a") > a Error in summary.connection(x) : invalid connection > When I go to yacas from the command line, it works fine. Any suggestions, please? I'm thinking that a port might be open, but here I have: >

I'm having trouble loading Ryacas. I've downloaded and extracted Ryacas 0.2-9 (also tried Ryacas ) and yacas 1.0.63, have the latest version of R and have tried the following (this works for installing other packages): install.packages("Ryacas") library(Ryacas) install.packages("yacas") library(yacas) and get this: Warning message: package ‘yacas’ is not available

I am having issues with Ryacas errors. I searched and found this error was reported in 2006, but nothing since. > library(Ryacas) Loading required package: XML > yacas(expression(Factor(x^2-1))) [1] "Starting Yacas!" CommandLine(1) : Expecting ) closing bracket for sub-expression, but got ^ instead > yacas("D(x)Sin(x)") CommandLine(1) : Expecting ) closing bracket

There is either something wrong with either me or is Yacas/Ryacas doing odd things. See below I ask to simpify an expression and the there's output! If this is a simplification.. well.. Do you think there is something set wrong somewhere. Thanks. > library(Ryacas) Loading required package: XML > a=Sym("a") > N=Sym("N") > b=Sym("b") >

I'm tinkering around in Ryacas trying to find an easy way to get the first and second partial derivatives of mu and sigma from a normal distribution (actually a bitterly ugly likelihood but this example works for now). I've done all of the work in Mathematica, but I then need to manually write R code whereas I think Ryacas will give me the expression such that I can export for Tex and

Hello again. Now I have the following: > xx [1] "Solve(1 - R/100==(2*y)/10000,R)" > yacas(xx) expression(list(R == 100 * (1 - 2 * y/10000))) > I would like to put in a value for y and obtain R. I've tried more stuff with eval and Eval, but no luck yet. Any suggestions would be much appreciated. Thanks, Erin -- Erin Hodgess Associate Professor Department of Computer

Let's say I had a matrix like this: library(Ryacas) x<-Sym("x") m<-matrix(c(cos (x), sin(x), -sin(x), cos(x)), ncol=2) How can I use R to obtain the eigenvalues and eigenvectors? Thanks, John [[alternative HTML version deleted]]