dxz2_2.m

数值类综合算法 常用数值计算工具包（龙贝格算法、改进欧拉法、龙格库塔方法、复合辛普森）

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('dxz2f',x0,[],vlb,[])