Integral CAS approximation

05012021, 08:25 PM
Post: #1




Integral CAS approximation
Curious concerning why the disagreement in approximation of an integral in CAS vs the same integral in HOME. When an exact integration is impossible for a function, should HOME be trusted to bring more accuracy? I was expecting an approximation in CAS would force a HOME type numerical approximation.
CAS approximation of integral from 1 to 1 of SQRT(1+x^3) produces 1.95275719206 In home, I get 1.95275723373 which I get on the WP 34S on a DM42 platform and 1.952757234 on a CASIO 9750GIII. 

05012021, 08:26 PM
Post: #2




RE: Integral CAS approximation
This on the beta April 28...


05022021, 05:51 AM
Post: #3




RE: Integral CAS approximation
You get identical results in both "HOME" and "CAS" if you use the extremes 1.0 and 1.0, instead of 1 and 1 (note the absence of the dot after 1 and 1). The same curious behavior can be found in "xCas" version 1.7


05022021, 06:30 AM
Post: #4




RE: Integral CAS approximation
HOME is using CAS for integral evaluation, but without trying an exact computation first (in other words it will call integrate(sqrt(1+x^3),x,1.0,1.0)).
If you run integrate(sqrt(1+x^3),x,1,1) in CAS, you get the initial integral expressed with other integrals because the CAS tried to get an exact answer, and this explains the small numeric differences, the general rule is that the precision is better without initial symbolic step (if the integral can not be expressed in closed form). 

05022021, 02:35 PM
Post: #5




RE: Integral CAS approximation
On HP emulator, integration seems work harder if we use pattern ∫(f(x), x= a .. b)
We integrate ∫(√(1+x^3), x= 1 .. 1) backwards (*), letting t = x Cas> expand(∫(√(1t^3), t, 1, 1)) √2*2*1/5 + integrate(3/5/√(t^3+1),t,1,1) Cas> float(Ans) → 1.95275719206 Cas> expand(∫(√(1t^3), t = 1 .. 1)) √2*2*1/5 + 1/5*√π*Gamma(1/3)/Gamma(5/6) + integrate(3/5/√(t^3+1),t,1,0) Cas> float(Ans) → 1.95275723373 (*) Cas knows ∫(√(1t^3), t,0,1), but not its equivalent, ∫(√(1+x^3), x,1,0) 

« Next Oldest  Next Newest »

User(s) browsing this thread: 1 Guest(s)