MathGroup Archive 1996

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

Search the Archive

optimal control

  • To: mathgroup at
  • Subject: [mg4839] optimal control
  • From: George Melikian <melikian at>
  • Date: Fri, 20 Sep 1996 01:12:52 -0400
  • Sender: owner-wri-mathgroup at

I am having some difficulty setting up an optimal control problem on
Mathematica. The problem is as follows:
The control variable (u) denotes a threshold setting (0<u<1) that is placed
between 2 probability densities (h & t). The objective is to maxize the
payoff (w) while considering changing risks over time (r(t)). Here is the

v = BetaDistribution[9,3];
r[t_] := PDF[v,t];
h = NormalDistribution[0,.1];
e = NormalDistribution[.3,.1];
g[u_] := CDF[h ,u];
b[u_] := CDF[e ,u];
w[u_,t_] := g[u]+(1-b[u])r[t]- b[u]r[t]- (1-g[u]);
te= Simplify[ D[w[u,t],t]];

hamiltonian = w[u,t] + lambda te;

The problem comes when I attempt to solve for the partial derivatives of
the hamiltonian with respect to: each  (u,t,lambda). 

D[Hamiltonian, lambda]
Solve[ D[Hamiltonian,lambda]==0, lambda] 

I get an error telling me either 1/0 encountered or that the function
cannot be solved. Any help here would be greatly appreciated.

Thanks in advance,
George Melikian
Melikian at


  • Prev by Date: FindRoot problem
  • Next by Date: Mathematica 3.0 Tour
  • Previous by thread: Re: FindRoot problem
  • Next by thread: Mathematica 3.0 Tour