+- HP Forums (https://www.hpmuseum.org/forum)
+-- Forum: HP Calculators (and very old HP Computers) (/forum-3.html)
+--- Forum: HP Prime (/forum-5.html)
+--- Thread: Polynomial derivative causes Prime to Crash? (/thread-6549.html)
Polynomial derivative causes Prime to Crash? - Spybot - 07-13-2016 05:05 PM
Hello!
Perhaps just trying to save some time while finding the derivative of this polynomial function:
d/dx(5x^(3/5)+4x^(3/4)) ... using the Prime Virtual Calculator, it failed causing the Emulator to crash, then I tried to solve it separating each term to find out the symbolic results are no the ones I expected.
[attachment=3753]
[attachment=3754]
I noticed that evaluating the derivative numerically it returns correct values... but symbolically I'm not sure, what do you think?
Thank you.
RE: Polynomial derivative causes Prime to Crash? - Arno K - 07-13-2016 05:16 PM
Well, the answer is correct, you got to simplify yourself.
Arno
RE: Polynomial derivative causes Prime to Crash? - parisse - 07-13-2016 06:19 PM
I don't think it crashes (at least it does not inside Xcas), but you would have to wait much too long, I will change some parameters to avoid that (maximum degree for common algebraic extension built for simplifications).
RE: Polynomial derivative causes Prime to Crash? - roadrunner - 07-13-2016 10:56 PM
In xcas diff(5*x^(3/5)+4*x^(3/4),x) returns:
3/(x^(1/5))^2+3*1/(x^(1/4))^3*sqrt(x)
However, simplify(diff(5*x^(3/5)+4*x^(3/4),x)) returns:
Code:
rootof([[1849290596989009920*x^3-144653156352000000*x^2+21432560400000000000*x+76904296875000,3178468213574860800*x^3+7094037381120000000*x^2+2440993500000000000*x-116729736328125,9635845027921920000*x^3-386546688000000000*x^2+229471875000000000*x,7879278565785600000*x^3-1874283840000000000*x^2-22113281250000000*x,(-9423782220272762880)*x^4+1713490624512000000*x^3-107266723200000000000*x^2+1757812500000000*x,(-23106632681898639360)*x^4-34563695837184000000*x^3-98018458200000000000*x^2+71411132812500*x,(-62818287878563430400)*x^4-25623209410560000000*x^3-10937965500000000000*x^2+487518310546875*x,(-72469455631810560000)*x^4+9139958784000000000*x^3-886802343750000000*x^2,19445417291078959104*x^5-40962633184051200000*x^4+223432480800000000000*x^3+86163574218750000*x^2,(-449515361885057187840)*x^5+101669760663552000000*x^4-5546714328000000000000*x^3-29462585449218750*x^2,(-708857389430907863040)*x^5-1855984556261376000000*x^4-502953623550000000000*x^3+31589813232421875*x^2,(-2431921742405015961600)*x^5+130928497781760000000*x^4-44452603406250000000*x^3-796508789062500*x^2,(-20610723694661074944)*x^6-1971674608293642240000*x^5+265932066816000000000*x^4+7454288671875000000*x^3,(-1129025475352763301888)*x^6-157730966234726400000*x^5-13308929464680000000000*x^4-767644409179687500*x^3,(-2769378177893825249280)*x^6-4396264306606080000000*x^5-8788632003900000000000*x^4+115511627197265625*x^3,(-6561699281264726507520)*x^6-2377594621145088000000*x^5-1056820966931250000000*x^4+32634887695312500*x^3,12198562540686409728*x^7-9550672312568433868800*x^6+1825856668567680000000*x^5-57416284212890625000*x^4+711822509765625*x^3,(-156421643493552685056)*x^7-2674904881032069120000*x^6-2496402776412000000000*x^5+2680500750732421875*x^4,(-138243318795485577216)*x^7-1439651046187008000000*x^6+2043538964595000000000*x^5-236938568115234375*x^4,(-298436070998016000000)*x^7+493404453398016000000*x^6-194979280837500000000*x^5+10898437500000000*x^4],[1,0,0,0,-5*x,-4*x,0,0,10*x^2,-260*x^2,6*x^2,0,-10*x^3,-620*x^3,-340*x^3,-4*x^3,5*x^4,-140*x^4,110*x^4,-20*x^4,-x^5+x^4]])with these warnings:
Warning, choice of an algebraic branch for root of a polynomial with parameters might be wrong. The choice is done for parameters value=0 if 0 is regular, otherwise randomly. Actual choice is Vector [92]
Warning, replacing 92 by 92.0, a substitution variable should perhaps be purged.
Warning, replacing 92 by 92.0, a substitution variable should perhaps be purged.
Warning, replacing 92 by 92.0, a substitution variable should perhaps be purged.
Warning, replacing 92 by 92.0, a substitution variable should perhaps be purged.
Warning, replacing 92 by 92.0, a substitution variable should perhaps be purged.
Warning, replacing 92 by 92.0, a substitution variable should perhaps be purged.
Warning, replacing 92 by 92.0, a substitution variable should perhaps be purged.
Warning, replacing 92 by 92.0, a substitution variable should perhaps be purged.
Evaluation time: 66.129
What am I doing wrong?
-road
RE: Polynomial derivative causes Prime to Crash? - Spybot - 07-14-2016 12:21 AM
Thank you Parisse, I appreciate you will do that for us.
RE: Polynomial derivative causes Prime to Crash? - compsystems - 07-14-2016 01:57 AM
3*x^(-2/5) == diff(5*x^(3/5)) [enter] returns 1 (true)
[up] [up] [copy] => 3*x^(-2/5) == 5*x^(3/5)' [enter] returns 1 (true) ok, but shifted the single quotation mark out of the entire expression, in this case not fails, because 5 is taken as constant
should be
[up] [up] [copy] =>3*x^(-2/5) == (5*x^(3/5))'
in the following case fails
diff(x/y) [enter] returns 1/y ok
but
[up] [up] [copy] => x/y' [enter] returns ±∞. BUG CONFIRMED
should be
[up] [up] [copy] => (x/y)' [enter] 1/y ok
FOR HP-PRIME TEAM
solution to bug
diff(expr) => parser as (expr)' and not as expr'
---------------------------------
for Bernard, we would see the solution as we make to naturally
diff(5*x^(3/5) => 3*x^(-2/5) or 3/x^(2/5) and not 3*(x^(1/5))^3/x (is not so simplified superficially)
as you can improve the simplifier?
RE: Polynomial derivative causes Prime to Crash? - parisse - 07-14-2016 05:15 AM
roadrunner: nothing is wrong. You get the result as rootof(P,Q) i.e. P(alpha) where Q(alpha)=0. This is because the simplifier rewrite all algebraic extensions in a unique one, in order to certify that it will not miss a simplification, the analog as when you add fractions like 1/2+1/3+1/6 you can get a simplification to an integer. However here it requires heavy computations.
I would recommend not to run simplify on expressions with fractional powers, by the way the title of the topic is misleading, this is not a polynomial derivative.
Regarding fractional powers, the CAS rewrites them as polynomials over simpler intermediate variables, here x^(3/5) as x^(1/5)^3.