reset;
model;
param r;
param m;
param M;
var x >=1 ;
var y >=1 ;
for{i in 1...34}
{
maximize U: (x[i]^0.27)*(y[i]^0.73);
subject to cons:(1+r[i])*x[i]+y[i]=m[i]*(1+r[i])+M[i];
data;
param m :=
199
233
248
226
214
222
219
196
176
175
181
169
140
140
155
167
122
134
144
154
165
176
192
210
230
255
287
328
386
472
574
635
721
;
param M :=
233
248
226
214
222
219
196
176
175
181
169
140
140
155
167
122
134
144
154
165
176
192
210
230
255
287
328
386
472
574
635
721
842
;
param r :=
0.0324
0.0378
0.0532
0.0576
0.0657
0.0684
0.0684
0.0699
0.0720
0.0730
0.0767
0.1111
0.0965
0.0789
0.0756
0.0937
0.1098
0.1098
0.0873
0.0715
0.0503
0.0294
0.0225
0.0225
0.0202
0.0198
0.0203
0.0225
0.0235
0.0320
0.0394
0.0225
0.0230
;
option solver minos;
solve;
}
option display_1col 0;
display x,m,M,r,y,u;
display x >> output7.txt;