Re: Help with NMinimize of NIntegrate function

• To: mathgroup at smc.vnet.net
• Subject: [mg86762] Re: Help with NMinimize of NIntegrate function
• From: antononcube <antononcube at gmail.com>
• Date: Thu, 20 Mar 2008 02:55:49 -0500 (EST)
• References: <fouedm\$m7f\$1@smc.vnet.net>

```The definition of the function EAa is wrong. Instead of

EAa[t, J, m, Q3, Q2]:= ...

it should be

EAa[t_, J_, m_, Q3_, Q2_]:=...

Then you won't get the NIntegrate::inumr messages.

Anton Antonov
Wolfram Research, Inc.

On Feb 13, 5:48 am, vasil michev <mic... at gmail.com> wrote:
> What's wrong with this code:
>
> test[t_?NumericQ, J_?NumericQ,
>   m_?NumericQ, (Q3_)?NumericQ, (Q2_)?NumericQ] :=
>  NIntegrate[
>     EAa[t, J, m, Q3, Q2]*
>      HeavisideTheta[EAa[t, J, m, Q3, Q2]], {x, -Pi, Pi}, {y, -Pi,
>      Pi}, {z, -Pi, Pi}, MaxPoints -> 1000000,
>     Method -> "LocalAdaptive"]/8/(Pi^3)
>
> In[75]:= test[10, 1, 1, 1, 1] // Timing
>
> Out[75]= {5.819, 2.22142}
>
> where
> EAa[t, J, m, Q3, Q2]:= 1/2 (-3 J - 2 m - 2 t (Cos[x] + Cos[y]) -
>     2 t (Cos[x]/3 + Cos[y]/3 + (4 Cos[z])/3)) -
>  1/2 Sqrt[16 (-(Q2/2) - (t Cos[x])/Sqrt[3] + (t Cos[y])/Sqrt[
>       3])^2 + (-2 Q3 + 2 t (Cos[x] + Cos[y]) -
>      2 t (Cos[x]/3 + Cos[y]/3 + (4 Cos[z])/3))^2]
>
> with
> NMinimize[test[10, 1, \[Mu], q, Q], {\[Mu], q, Q}]
> I get those inumr errors, have i set anything wrong? Everything should
> be numerical, x,y,z are integrated out, so what is it?

```

• Prev by Date: Re: choose elements from list based on elements in different list
• Next by Date: Re: choose elements from list based on elements in different list
• Previous by thread: Re: Creating a function that is a derivative of another
• Next by thread: Saving Packages (was: Re: importing nb files)