Mathematica 9 is now available
Services & Resources / Wolfram Forums
-----
 /
MathGroup Archive
2005
*January
*February
*March
*April
*May
*June
*July
*August
*September
*October
*November
*December
*Archive Index
*Ask about this page
*Print this page
*Give us feedback
*Sign up for the Wolfram Insider

MathGroup Archive 2005

[Date Index] [Thread Index] [Author Index]

Search the Archive

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


  • Prev by Date: Re: Re: 3d plots in mathematica 5.0
  • Next by Date: Re: Add terms surrounded by zero together in matrix
  • Previous by thread: Re: NIntegrate::inum
  • Next by thread: Colorfunction along a line...