Re: NIntegrate::inum
- To: mathgroup at smc.vnet.net
- Subject: [mg59128] Re: NIntegrate::inum
- From: Peter Pein <petsie at dordos.net>
- Date: Sat, 30 Jul 2005 01:25:02 -0400 (EDT)
- References: <dca01n$cul$1@smc.vnet.net> <dcccag$3ch$1@smc.vnet.net>
- Sender: owner-wri-mathgroup at wolfram.com
Valeri Astanoff schrieb:
> If you replace f[x_, om_, w_] by f[x_?NumberQ, om_, w_]
> it settles the trouble.
>
> v.a.
>
If I try this with Ruth's code, I get a kernel crash in version 5.1
(Mathematica's kernel, of course)!?. If all NumberQs are removed, except
in rr[zz_?NumericQ, om_?NumericQ, w_?NumericQ]:=..., the results are
dependent of the starting value of om:
In[9]:=
Timing[NMinimize[{chi2f2[om, w], 0 <= om <= 1}, {om, w}]]
Out[9]=
{12.891*Second, {0.2633701433313954,
{om -> 0.9999982304774572, w -> 0.07381086751638463}}}
If NMinimize starts with an om between 0 and 1/2, we get the same
minimum as in my previous posting:
Timing[NMinimize[{chi2f2[om, w], 0 <= om <= 1}, {{om, 0, 1/2}, w},
MaxIterations -> 1000]]
Out[13]=
{10.328*Second, {0.14773793935404245,
{om -> 0., w -> -2.5414386386395766}}}
--
Peter Pein
Berlin