From: Andrii S. <and...@gm...> - 2023-01-23 09:30:26
|
Thanks for feedback and some testing. Created a bug report https://sourceforge.net/p/maxima/bugs/4084/ Andrii On Sun, 2023-01-22 at 18:16 -0800, Robert Dodier wrote: > On Sun, Jan 22, 2023 at 9:40 AM Andrii Sokolov <and...@gm...> wrote: > > > Recently I was not able to use 'solve' for an elementary linear equation--- > > if 'orderless' was called before. The commands are invoked are the following > > (sorry for a lengthy input): > > > > orderless(w1, w2, W1, W2, k1, k2, g, f1, f2, K)$ > > > > eq : -%i*W1*B10 = K * ((%i*g*f1^3)/((%i*w1-%i*W1+k1)*((-%i*w1)+%i*W1+k1) > [...] > > *k2^2)$ > > > > solve(eq, B10); > > > > which spits out: > > > > PQUOTIENT: Quotient by a polynomial of higher degree (case 1) > > Hmm, I can't reproduce this error. I tried it with Maxima 5.46.0 and > current Git version (post-5.46) with SBCL on Linux, and tried it with > orderless and without, and in all cases, solve(eq, B10) concludes > without error, and the result appears to be correct -- I see that > rhs(eq) is free of B10, so the result should be just rhs(eq)/(-%i*W1), > which I tested with ratsimp(rhs(eq)/(-%i*W1) - rhs(%o3[1])) where %o3 > is the output from solve, and got 0 as expected. > > Not sure what's going on here. Can you create a bug report? Maybe > someone can try it with GCL. > > best, > > Robert |