jm98a3fun.m

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

function [f,g]=jm98a3fun(xx)
%1998年全国大学生数学建模竞赛A题：收益与风险 非线性模型优化函数
%《数学的实践与认识》p39-42
global M r q p u lemda;
xx=xx(:);len=length(xx);x=xx(1:(len-1));
y=(x>100*eps).*(x<u/M/100).*x*100+(x>=u/M/100).*max(x,u/M);
y=y.*p;
f=lemda*xx(len)-(1-lemda)*sum(r.*x-y);
g(1)=sum(x+y)-1;
g(2:(len-1))=q(1:(len-2)).*x(1:(len-2))-xx(len);g=g(:);