MathGroup Archive 2006

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

Search the Archive

problems with Nmaximization

  • To: mathgroup at smc.vnet.net
  • Subject: [mg67815] problems with Nmaximization
  • From: "narsi" <narasimhan.sowmyanarayanan at gmail.com>
  • Date: Sat, 8 Jul 2006 04:57:02 -0400 (EDT)
  • Sender: owner-wri-mathgroup at wolfram.com

I am trying to use nmax to evaluate a function using Nmaximize. The
function is as follows.

\!\(m[1] = 75\[IndentingNewLine]
  m[2] = 36\[IndentingNewLine]
  m[3] = 12\[IndentingNewLine]
  m[4] = 4\[IndentingNewLine]
  m[5] = 2\[IndentingNewLine]
  m[6] = 0\[IndentingNewLine]
  m[7] = 0\[IndentingNewLine]
  m[8] = 0\[IndentingNewLine]
  m[9] = 0\[IndentingNewLine]
  n[1] = 571 - 75\[IndentingNewLine]
  n[2] = 338 - 36\[IndentingNewLine]
  n[3] = 173 - 12\[IndentingNewLine]
  n[4] = 93 - 4\[IndentingNewLine]
  n[5] = 41 - 2\[IndentingNewLine]
  n[6] = 21\[IndentingNewLine]
  n[7] = 10\[IndentingNewLine]
  n[8] = 6\[IndentingNewLine]
  n[9] = 4\[IndentingNewLine]\[IndentingNewLine]
  NMaximize[{â??\+\(t = 1\)\%9\((m[
  t]*Log[\((\(h\ *p*\((1 -
      p)\)\^\((t - 1)\)\)\/\(1 - h + h*\((1 - p)\)\^\((t - 1)\)\))\)] -
     n[t]*Log[1 - \((\(h\ *p*\((1 - p)\)\^\((t - 1)\)\)\/\(1 - h +
   h*\((1 - p)\)\^\((t - 1)\)\))\)])\), 1 â?¥ h â?¥ 0, 1 â?¥ p â?¥ 0},
{h, p}]\)

I cannot get this to converge. It gives me the following message

\!\(ower::"infy" \(\(:\)\(\ \)\) "
    Infinite expression \!\(1\/0.`\) encountered"\)


Can someone please help with this error in mathematica. I dont know
what the problem is. 

Thanks

narsi


  • Prev by Date: Re: 3D graphics headache
  • Next by Date: ReplacePart in an If[] construct?
  • Previous by thread: Re: failure to Integrate in orthogonal polynomials
  • Next by thread: Re: problems with Nmaximization