<optimum}} The code below extracts the two elements *) {tmpResults[[1]],(c /. tmpResults[[2]])} ,{i,Length[xGrid]}]; MaxValuesTm1 = -Transpose[MaxResultsTm1][[1]]; (* Extract the maxima *) cValuesTm1 = Transpose[MaxResultsTm1][[2]]; (* Extract the optimal c's *) (* Construct the numerical interpolating functions *) VTm1 = Interpolation[Transpose[{xGrid,MaxValuesTm1}],InterpolationOrder->1]; cTm1 = Interpolation[Transpose[{xGrid,cValuesTm1}],InterpolationOrder->1]; cTm1CEQ[x_] := .5 + x/2; (* Optimal policy if there is no uncertainty is consume half *) PlotVTm1Simple = Plot[VTm1[x],{x,xMin,xMax} ,PlotLabel->"VTm1" ,ImageSize->{72 6.,72 6./GoldenRatio} ]; Display[ParentDirectoryString<>"PlotVTm1Simple.EPS",PlotVTm1Simple,"EPS"]; PlotcTm1Simple = Plot[{cTm1[x],cTm1CEQ[x]},{x,xMin,xMax} ,PlotStyle->{Dashing[{}],GrayLevel[.5]} ,PlotLabel->"cTm1[XTm1] (solid line) vs. cTm1CEQ[XTm1] (gray)" ,ImageSize->{72 6.,72 6./GoldenRatio} ]; Display[ParentDirectoryString<>"PlotcTm1Simple.EPS",PlotcTm1Simple,"EPS"];