Re: Integrate this how?
- To: mathgroup at smc.vnet.net
- Subject: [mg116329] Re: Integrate this how?
- From: DrMajorBob <btreat1 at austin.rr.com>
- Date: Fri, 11 Feb 2011 04:15:00 -0500 (EST)
I get this in version 8 (with error messages for one or more methods): r[t_] = (1 + Sin[t]) (1 + 9/10 Cos[8 t]) (1 + 1/10 Cos[24 t]) (9/10 + 5/100 Cos[200 t]); arcLength[t_] = Sqrt[r[t]^2 + D[r[t], t]^2] // Simplify; NIntegrate[arcLength@t, {t, -Pi, Pi}, Method -> #] & /@ {"GlobalAdaptive", "LocalAdaptive", "MonteCarlo", "AdaptiveMonteCarlo", "QuasiMonteCarlo", "AdaptiveQuasiMonteCarlo"} {48.37, 48.0666, 48.0852, 48.3858, 47.985, 48.0663} but Integrate simply fails. The following plot suggests how messy this is, numerically: Plot[arcLength@t, {t, -Pi, Pi}] I think you're trying to approximate a fractal, and that fractal has infinite arc length... so it's not surprising you run into difficulties. Bobby On Thu, 10 Feb 2011 04:21:36 -0600, BoLe <bole79 at gmail.com> wrote: > Curve is closed (t from -Pi to Pi) and in polar form > > r[t_]:=(1+Sin[t]) (1+9/10 Cos[8 t]) (1+1/10 Cos[24 t]) (9/10+5/100 > Cos[200 t]) > > Would like to find out its arc length, so the integral of > > Sqrt[r[t]^2+D[r[t],t]^2] > > Version 7 complains but returns a somewhat reasonable number > (something around 50; drew together with a circle of the same length), > version 8 says nothing but gives just too low a number (five and > something). How can this be? New Mathematica should have no problem > with such high oscillatory integrands, shouldn't she? > -- DrMajorBob at yahoo.com