dxz2f.m

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

function [f,g]=fun(xx)
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(:);