%1998年全國大學(xué)生數(shù)學(xué)建模競賽A題：收益與風(fēng)險(xiǎn)非線性模型
%參看《數(shù)學(xué)的實(shí)踐與認(rèn)識(shí)》1999年第一期,p39-42
%另見 jm98a1, jm98a1,jm98a3fun

%問題一
clear;global M r q p u g lemda;
M=100;lemda=0.9;
r=[28 21 23 25 5]'/100;
q=[2.5 1.5 5.5 2.6 0]'/100;
p=[1 2 4.5 6.5 0]'/100;
u=[103 198 52 40 1]'/100;
x0=0.5*ones(length(r)+1,1)/6;
vlb=zeros(6,1);
[x options]=constr('jm98a3fun',x0,[],vlb,[])

%問題二
clear;global M r q p u g lemda;
M=2000;
lemda=0.3;
r=[9.6 18.5 49.4 23.9 8.1 14 40.7 31.2 33.6 36.8 11.8 9 35 9.4 15 5]'/100;
q=[42 54 60 42 1.2 39 68 33.4 53.3 40 31 5.5 46 5.3 23 0]'/100;
p=[2.1 3.2 6 1.5 7.6 3.4 5.6 3.1 2.7 2.9 5.1 5.7 2.7 4.5 7.6 0]'/100;
u=[181 407 428 549 270 397 178 220 475 248 195 320 267 328 131 0]';
x0=0.5*ones(length(r)+1,1)/(length(r)+1);
vlb=zeros(length(r)+1,1);
[x options]=constr('jm98a3fun',x0,[],vlb,[])