dharr

Dr. David Harrington

6399 Reputation

21 Badges

20 years, 17 days
University of Victoria
Professor or university staff
Victoria, British Columbia, Canada

Social Networks and Content at Maplesoft.com

I am a retired professor of chemistry at the University of Victoria, BC, Canada. My research areas are electrochemistry and surface science. I have been a user of Maple since about 1990.

MaplePrimes Activity


These are replies submitted by dharr

Glad it worked for you. Just a tip: diff(y(x), x, x, x, x); is simpler as diff(y(x), x, x, x, x); which can be further abbreviated as diff(y(x), x$4);
That's interesting, and very mysterious. Do you get it OK with my example? I am using Maple 10 under Windows. The only other thing I can think of would be that using a list might force the order, and I was just lucky skipping it. display([points, curve]); display([curve,points]);
That's interesting, and very mysterious. Do you get it OK with my example? I am using Maple 10 under Windows. The only other thing I can think of would be that using a list might force the order, and I was just lucky skipping it. display([points, curve]); display([curve,points]);
Thanks, this works. I look forward to Maple 11 then.
Thanks, this works. I look forward to Maple 11 then.
Yes, sorry about that.
Yes, sorry about that.
Thanks for pointing me there, I would never have found. I can now do what I want to. BTW, you can save at least simple *.mw files as help pages with makehelp (which just reads the lines and calls INTERFACE_HELP), but not sure which features you might lose.
Thanks for pointing me there, I would never have found. I can now do what I want to. BTW, you can save at least simple *.mw files as help pages with makehelp (which just reads the lines and calls INTERFACE_HELP), but not sure which features you might lose.
I agree about the difficulty and annoyance in a long code of finding such problems. But qualified responses could quickly get complicated. Consider solve(a*x^2+b*x+c,x); This needs a qualification that a<>0, which is simple enough. But should it also be mentioned that there is one answer, not two, if b^2=4*a*c? And that if a=0, then the solution is -c/b unless b=0, in which case the equation doesn't have x in it any more ... Cheers, David.
I think Maple frequently gives results that have restrictions; a simpler example is simplify(c/c), which is also not true for c=0. Likewise solve(1=(x+1)/(x+1),x); returns x, but x=-1 is not acceptable. If Maple is too fussy, then one never gets to a useful result, so I don't consider this a bug as such.
Try with(StringTools): sentence:="This sentence, isn't interesting"; Split(sentence);remove(type,%,"");Join( %, " " ); Words(sentence);Join( %, " " ); Split keeps punctuation, but Words throws it away, so it depends on what you want.
Sorry, I was thinking all roots real. If the signs of the coefficients are known, then Descartes' rule of signs can tell max number of real roots, but there can be less, so finding when there is an exact number of real roots is a pretty difficult problem.
Try: R:=alpha->3*alpha^6-6*alpha^5*delta+(6*delta*a-c_p+mu^2-4*a^2)*alpha^4+mu^2*(6*delta*a-c_p+mu^2-4*a^2)*alpha^2-6*delta*alpha*mu^4+3*mu^6: with(PolynomialTools): Hurwitz(R(alpha),alpha,'s','g'); s; which gives the output [0, -1/2*alpha/delta, -6*alpha*delta/(6*delta*a-c_p+mu^2-4*a^2), 1/6*(6*delta*a-c_p+mu^2-4*a^2)^2*alpha/(mu^2*(-2*mu^2-4*a^2+6*delta*a-c_p)*delta), (-12*mu^2-24*a^2+36*delta*a-6*c_p)*delta*alpha/((6*delta*a-c_p+mu^2-4*a^2)*(-2*c_p-mu^2-8*a^2+12*delta*a)), -1/12*(-2*c_p-mu^2-8*a^2+12*delta*a)^2*alpha/(mu^2*(-2*mu^2-4*a^2+6*delta*a-c_p)*delta), -4*(-2*mu^2-4*a^2+6*delta*a-c_p)*delta*alpha/(mu^2*(-2*c_p-mu^2-8*a^2+12*delta*a))] I haven't used this much, but understand that for real roots you need the coefficients of alpha in the (Stieltjes) sequence (after the first entry) all to be positive. So looking at the second entry: -1/2*alpha/delta delta must be negative. So now looking at the next one: -6*alpha*delta/(6*delta*a-c_p+mu^2-4*a^2) Since delta is negative the denomimator must be negative to get this positive, etc, etc. So you can build up conditions on the parameters to get real roots. More information of this method is given in the reference at the end of the Hurwitz help page: N. Levinson and R.M. Redhoffer, Complex Variable, 1970.
Hi, Thomas, I think I owe you one. Hard to say without seeing it, but sounds like you just need to force Maple to re-order it, perhaps to get terms close enough. I'd try collect (in z1 or some combination of the variables; probably not important which) and then simplify again. Sort might work as well. Cheers, David.
First 62 63 64 65 66 Page 64 of 66