Re: Re: More weird integration issues...
- To: mathgroup at smc.vnet.net
- Subject: [mg35765] Re: [mg35748] Re: [mg35728] More weird integration issues...
- From: Andrzej Kozlowski <andrzej at tuins.ac.jp>
- Date: Tue, 30 Jul 2002 07:22:15 -0400 (EDT)
- Sender: owner-wri-mathgroup at wolfram.com
It seems that you have not evaluated Binesh's defintion of F. You should be getting a numerical answer, not a formula involving F, Integrate etc. Andrzej On Monday, July 29, 2002, at 07:57 PM, DrBob wrote: > I get the same answer both ways, to as many places as you see here: > > 1.8430622009642987*Integrate[ > (-0.9025328739879791 + X)^5* > F[Subscript[p, 0], > 0.2856264419666387, > 0.04012528894153399, > 0.2662836604257217, > 0.7959934248376723, > 0.15354837805257518, X], > {X, -Infinity, Infinity}] > > I'm using Version 4.2. > > Bobby > > -----Original Message----- > From: Andrzej Kozlowski [mailto:andrzej at tuins.ac.jp] To: mathgroup at smc.vnet.net > Sent: Monday, July 29, 2002 2:13 AM > Subject: [mg35765] [mg35748] Re: [mg35728] More weird integration issues... > > There does sem to be a bug here, but it is not quite what you think. > It's the first answer that you get that is probably wrong. On the other > hand, in the second case it is very unlikely that Mathematica enters an > infinite loop, rather it is still trying to arrive at the answer and > there is no guarantee that it will reach one after, say a week or a > month. > > As for first case, the reason why the answer is probably this. Evaluate > the formula: > > > formula=Integrate[F[Subscript[p, 0], Subscript[p, 1], Subscript[p, 2], > Subscript[p, 3], Subscript[p, 4], Subscript[p, 5], X]*((X - > m)/sd)^5, > {X, -Infinity, Infinity}] > > Now set > > > Evaluate[Table[Subscript[p, i], {i, 1, 5}]] = Table[Random[], {5}] > > and also > > sd = Random[]; m = Random[]; > > Now evaluate again > > > Integrate[F[Subscript[p, 0], Subscript[p, 1], Subscript[p, 2], > Subscript[p, 3], Subscript[p, 4], Subscript[p, 5], X]*((X - > m)/sd)^5, > {X, -Infinity, Infinity}] > > and > > formula > > You will almost certainly get different answers, while they clearly > ought to be the same. It seems that it is nto the fact that the names > you are suing are different that leads to different results in both of > your integrals but the fact that the names of the parameters in the > first case are not symbols. "Officially" there is no reason why they > should be, but in practice using non-symbols in formulas makes them more > > complicated and is more likely to result in errors. > > The really bad news as far as your problem is concerned is that it is > the 6 hour fruitless computation that appears to be the correct one ... > Andrzej Kozlowski Toyama International University JAPAN http://platon.c.u-tokyo.ac.jp/andrzej/ > > On Saturday, July 27, 2002, at 07:43 PM, Binesh Bannerjee wrote: > >> Hi. I'm trying to compute the 5th moment of a tweak to the normal pdf. >> Here's what I have: >> >> F[a_, b_, c_, d_, e_, f_, X_] := (a + b X + c X^2 + d X^3 + e X^4 + f >> X^5)* >> Exp[-(((X - m)/sd)^2)/2]/(Sqrt[2Pi]sd) >> >> The strange thing, and I'd appreciate someone shedding some light on >> this, >> is that this: >> >> \!\(Integrate[ >> F[p\_0, p\_1, p\_2, p\_3, p\_4, p\_5, >> X]*\((\((X - m)\)/sd)\)^5, {X, \(-Infinity\), Infinity}, >> Assumptions -> {sd > 0}]\) >> >> (I cut and pasted that it looks like it works) >> >> Anyway, THAT gives me an answer really quickly... (within 5 minutes on > >> my box) >> >> JUST changing it from p0..5 to a,b,c,d,e,f like so: >> >> Integrate[F[a, b, c, d, e, f, X]*((X - m)/sd)^5, {X, -Infinity, >> Infinity}, >> Assumptions -> {sd > 0}] >> >> Causes mathematica to go into an infinite loop (seemingly after 6 >> hours). >> >> This ... sucks. How am I to know if a certain equation is solvable, if >> only I choose the right variables?? >> >> Binesh Bannerjee >> >> -- >> "For in much wisdom is much grief, and he that increaseth knowledge >> increaseth sorrow." -- Ecclesiastes 1:18 >> >> PGP Key: http://www.hex21.com/~binesh/binesh-public.asc >> SSH2 Key: http://www.hex21.com/~binesh/binesh-ssh2.pub >> SSH1 Key: http://www.hex21.com/~binesh/binesh-ssh1.pub >> OpenSSH Key: http://www.hex21.com/~binesh/binesh-openssh.pub >> >> >> > > > > >