First time here? Check out the FAQ!

Ask Your Question
0

error in factor

asked 0 years ago

pfeifhns gravatar image

updated 0 years ago

Emmanuel Charpentier gravatar image

Hi, I have a problem with the factor() function: Consider the following code:

%display latex

j_00 = 1/2*(a^2 + eta)*(a^2 + zeta)/((a^2 - b^2)*(a^2 - c^2)*sqrt((a^2 + eta)*(a^2 + xi)*(a^2 + zeta)/((a^2 - b^2)*(a^2 - c^2))))

show(j_00)

j_00_simplified = j_00.canonicalize_radical().factor())

print(j_00_simplified)

As a result, I get:

1/2*sqrt(a^2 + eta)*sqrt(a^2 + zeta)/(sqrt(a^2 + xi)*sqrt(a + b)*sqrt(a + c)*sqrt(-a + b)*sqrt(-a + c))

The result is false! It should be:

1/2*sqrt(a^2 + eta)*sqrt(a^2 + zeta)/(sqrt(a^2 + xi)*sqrt(b^2-a^2)*sqrt(a^2 - c^2))

The term -a+c is wrong, it should be c-a. I think this is an error, or is there anything I don't see? If it is an error, what can I do that factor() returns the expected result?

Thank's in adcance,

Hans

Preview: (hide)

Comments

Edited for legibility.

When you need to insert a block of code, indent it by four spaces (or select it an hit Ctrl-K) (or use the 101010 button). Inline code can be surrounded by backticks (`).

Emmanuel Charpentier gravatar imageEmmanuel Charpentier ( 0 years ago )

[ Edited as an answer. Sorry for the noise... ].

Emmanuel Charpentier gravatar imageEmmanuel Charpentier ( 0 years ago )

What's the difference between -a+c and c-a? Maybe I'm misunderstanding.

John Palmieri gravatar imageJohn Palmieri ( 0 years ago )

1 Answer

Sort by » oldest newest most voted
0

answered 0 years ago

Emmanuel Charpentier gravatar image

updated 0 years ago

Tentative answer : on Sage 10.5.rc0, I get :

sage: j_00_simplified
1/2*sqrt(a^2 + eta)*sqrt(a^2 + zeta)/(sqrt(a^2 + xi)*sqrt(a + b)*sqrt(a + c)*sqrt(-a + b)*sqrt(-a + c))

and, BTW,

sage: (j_00.canonicalize_radical().factor()/j_00).canonicalize_radical()
1

What is your version of Sage ? On what platform ?

EDIT : Your expression is symmetric in b anc c, therefore so your simplified expression should be ; the "correct" expression you propose is not...

Preview: (hide)
link

Your Answer

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

Add Answer

Question Tools

1 follower

Stats

Asked: 0 years ago

Seen: 176 times

Last updated: Nov 25 '24