From: Robert D. <rob...@gm...> - 2023-01-23 02:16:25
|
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 |