Ask Your Question

# algebraic substitution

How can I imitate the algsubs function of Maple? A little example:

p = x^18

algsubs(x^2=x+1,p)

The reult in Maple is:

x^9+9x^8+36x^7+84x^6+126x^5+126x^4+84x^3+36x^2+9x+1.

edit retag close merge delete

## Comments

What is the mathematical meaning of this function ? If you want x^2 to be equal to x+1, then the optimal result would be a polynomial of degree 1, namely 2584*x + 1597.

## 2 Answers

Sort by » oldest newest most voted

It seems that there are not the same function in sagemath but we can use " ratsubst " function of maxima

sage : maxima('p:x^18')
sage : maxima('ratsubst(a+1,x^2,p)')


how to use certain maxima functions (e.g. ratsubst)?

more var('y')
P = x^18
eqn = x^2 == y + 1
soln = solve(eqn,x)
P = P.subs(x=soln.rhs())
P.expand().subs(y=x)

more

## Comments

For example, if P = x^17 = x^16*x=(x^2)^4*x=(x+1)^4*x. But it is impossible with this method.

Out of curiosity, what do you want when you have x^{17}? Something with a square root? And if so, which square root?

## Your Answer

Please start posting anonymously - your entry will be published after you log in or create a new account.

Add Answer

## Stats

Asked: 2014-02-06 11:38:29 +0200

Seen: 626 times

Last updated: Feb 07 '14