dxz2_1.m

来自「数值类综合算法 常用数值计算工具包（龙贝格算法、改进欧拉法、龙格库塔方法、」· M 代码 · 共 10 行

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