Constrained Nonlinear Optimization with expectations functions in the constraints

I am modeling an optimal compensation contract with a CRRA utility function

sol= Minimize[{a + s*100 + o*25, Expectation[(1/(1-c))((((w + a) Exp[rT]) + (s Exp[dT] +ol)*p Exp[(r-d-(1/2)q^2)T +(xSqrt[T]q)])^(1 – c)), x [Distributed] NormalDistribution[0, 1]]=-0.1349, Expectation[((((a + w) Exp[r* T]) + (p (lo + sExp[dT])* Exp[(-d-q^2*(1/2)+ r)T + qSqrt[T]x]))^(-c))(lo+ sExp[d*T])Exp[(-d-q^2(1/2)+r)T+qSqrt[T]*x], x [Distributed] NormalDistribution[0, 1]]>=0.00075, a >= 0, 0 <= s <= 1, 0 <= o <= 1}, {a,s,o}]

I need to find the results of expectation in the form of function(a,s,o) to be used in the constraints. This program have to be solved for multiples values of the parameters; in one case we have {{d -> 0.0269, q -> 0.3315, r -> 0.0435, T -> 1, w -> 1.1465, p -> 100, l -> 1, c -> 3}} and the initial values of a, s and o are: {{a -> 0.0941, s -> 0.0029, o -> 0.0089}} It is used to calculate the value of the constraint. I need to find the results of expectation in the form of function(a,s,o) to be used in the constraints. thank you very much in advance. Sawsen

This entry passed through the Full-Text RSS service – if this is your content and you’re reading it on someone else’s site, please read the FAQ at fivefilters.org/content-only/faq.php#publishers.

Recent Questions – Mathematica Stack Exchange