HP Forums

Full Version: Integral CAS approximation
You're currently viewing a stripped down version of our content. View the full version with proper formatting.
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.
This on the beta April 28...
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
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).
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(∫(√(1-t^3), t, -1, 1))

√2*2*1/5 + integrate(3/5/√(-t^3+1),t,-1,1)

Cas> float(Ans)       → 1.95275719206

Cas> expand(∫(√(1-t^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 ∫(√(1-t^3), t,0,1), but not its equivalent, ∫(√(1+x^3), x,-1,0)
Reference URL's