Skip to main content
8 events
when toggle format what by license comment
Dec 24, 2015 at 19:46 history edited Jack LaVigne CC BY-SA 3.0
added 3711 characters in body
Dec 24, 2015 at 16:00 comment added Jerry Guern No, I mean the 1 in the first argument of PolynomialRemainder in the code line for q2, not q1. Read the line of code right after the line in my question about generalizing.
Dec 24, 2015 at 15:45 comment added Jack LaVigne Your posted question addresses NSolve rather than Solve. Solve is giving the correct answer. If you numerically evaluate them some answers appear wrong. In particular solution[[5]]. {q1, q2} /. N[solution[[5]]] is bad. It turns out that this is a matter of precision. {q1, q2} /. N[solution[[5]], 30] gives the correct solution.
Dec 24, 2015 at 15:42 comment added Jack LaVigne I think you meant replace 1 in q1 with t rather than q2. I don't see a 1 in q2. I tried replacing 1 in q1 with t and Solve spit out a result but I don't think there is any hope of being able to plot it.
Dec 24, 2015 at 3:49 comment added Jerry Guern Actually, this is returning incorrect results. I posted a Question about it.
Dec 24, 2015 at 0:27 comment added Jerry Guern It also really surprised me that this returns complex roots that aren't paired with their complex conjugate. I know that if a single-variable polynomial has real coefficients, it's complex roots will be paired with their conjugate. Was I confused when I expected that here?
Dec 24, 2015 at 0:26 comment added Jerry Guern Oh, I was trying to eliminate p and didn't think of just solving for {c,p} as a pair. Any thoughts on the generalization part, so I can plot c[t] where Abs[t] in some reasonable time?
Dec 24, 2015 at 0:07 history answered Jack LaVigne CC BY-SA 3.0