Is it possible to solve an algebraic equation in R?

M. Beausoleil picture M. Beausoleil · Sep 16, 2015 · Viewed 8.7k times · Source

I want to find the solution of:

-x^3+6*x^2+51*x+44=0

but with R. Is it possible?

I found the package Ryacas, but nobody seems to be able to make it work.

May sound trivial, but I'm not able to find an easy way to do this...

Do you have an alternative?

Thanks guys!

Answer

user3710546 picture user3710546 · Sep 16, 2015

You can use polynom package:

library(polynom)
p <- polynomial(c(44,51,6,-1))
# 44 + 51*x + 6*x^2 - x^3 
solve(p)
# [1] -4 -1 11

But you simply can use the function polyroot from base package:

polyroot(c(44,51,6,-1))
# [1] -1+0i -4+0i 11+0i

If you keep the real part with Re:

Re(polyroot(c(44,51,6,-1)))
# [1] -1 -4 11