variation in Um for household members

KNITRO 9.0.0: Locally optimal solution.
objective -364296.7893; feasibility error 1.27e-09
223 iterations; 444 function evaluations
_solve_time = 8.280000

beta = 0.950000
theta = 0.200000
VoT = 60.000000

:         Um             b            c         :=
1 1   2099.669581      0.000000   319.949983
1 2   1772.003815    731.299728   177.286912
1 3    415.751772   1104.578470    67.722766
2 1   2149.214930      0.000000   386.260732
2 2   1460.631835    669.370193   242.959757
2 3    190.966813   1109.107776    72.178958
;

rho [*] :=
1   0.218584
2  -0.017779
3   0.156557
;
tweak solver options for estimating joint MDP

KNITRO 9.0.0: Locally optimal solution.
objective -369866.8647; feasibility error 6.21e-08
242 iterations; 500 function evaluations
_solve_time = 8.890000

beta = 0.950000
theta = 0.200000
VoT = 60.000000

:         Um             b            c         :=
1 1   2243.027989      0.000000   334.493799
1 2   1895.784633    734.336305   181.913438
1 3    448.532516   1104.170709    70.699313
2 1   2191.269612      0.000000   390.403832
2 2   1509.580667    669.152599   239.391992
2 3    221.667755   1108.995724    71.590607
;

rho [*] :=
1   0.208258
2  -0.034983
3   0.052484
;

:         Um0           b0            c0        :=
1 1   3600.000000      0.000000   320.000000
1 2   2500.000000    840.000000   180.000000
1 3   2000.000000   1140.000000   210.000000
2 1   3600.000000      0.000000   320.000000
2 2   2500.000000    840.000000   180.000000
2 3   1500.000000   1140.000000   210.000000
;

rho0 [*] :=
1  0.300000
2  0.000000
3  0.200000
;