程序出现如下问题，---options.MaxFunEvals = 400 (the default value).论文时间紧张，有时间的同学可以看一下吗？谢谢

x0 = [100;100;100;100];

A = [];

b = [];

Aeq = [1,1,1,1];

beq = [2000];

VLB = [300;400;0;0];

VUB = [600;1000;500;500];

[x,fval] =fmincon(@ht,x0,A,b,Aeq,beq,VLB,VUB)

function f =ht(x);

a=5;

b=5;

f = 0.5*x(1)^2+1.25*x(1)+1+...,

0.625*x(2)^2+1*x(2)+1.25+...,

1*x(3)+1.05*x(4)+...,

a*(.717e-11*(-.37e10*exp(x(3))^10-.162e11*x(3)^2*exp(x(3))^10+.559e11*x(3)*exp(x(3))^10+.401e11+.681e13*x(3)^4+.134e13*x(3)^2+.583e13*x(3)^9-.464e11*erf(.333+.667*x(3))*exp(x(3))^10+.121e14*x(3)^6+.329e12*x(3)^12+.102e14*x(3)^5+.116e13*x(3)^11+.296e13*x(3)^10+.355e13*x(3)^3+.332e12*x(3)+.917e13*x(3)^8+.539e10*x(3)^14+.117e14*x(3)^7+.604e11*x(3)^13)/exp(x(3))^10)+...,

a*(.717e-11*(-.37e10*exp(x(4))^10-.162e11*x(4)^2*exp(x(4))^10+.559e11*x(4)*exp(x(4))^10+.401e11+.681e13*x(4)^4+.134e13*x(4)^2+.583e13*x(4)^9-.464e11*erf(.333+.667*x(4))*exp(x(4))^10+.121e14*x(4)^6+.329e12*x(4)^12+.102e14*x(4)^5+.116e13*x(4)^11+.296e13*x(4)^10+.355e13*x(4)^3+.332e12*x(4)+.917e13*x(4)^8+.539e10*x(4)^14+.117e14*x(4)^7+.604e11*x(4)^13)/exp(x(4))^10)+...,

b*(.717e-11*(-.401e11-.464e11*erf(.333+.667*x(3))*exp(x(3))^10-.917e13*x(3)^8-.539e10*x(3)^14-.117e14*x(3)^7-.604e11*x(3)^13-.121e14*x(3)^6-.329e12*x(3)^12-.102e14*x(3)^5-.116e13*x(3)^11-.681e13*x(3)^4-.296e13*x(3)^10-.355e13*x(3)^3-.583e13*x(3)^9-.134e13*x(3)^2-.332e12*x(3)-.389e11*x(3)*exp(x(3))^10+.569e11*exp(x(3))^10+.348e11*x(3)^2*exp(x(3))^10)/exp(x(3))^10)+...,

b*(.717e-11*(-.401e11-.464e11*erf(.333+.667*x(4))*exp(x(4))^10-.917e13*x(4)^8-.539e10*x(4)^14-.117e14*x(4)^7-.604e11*x(4)^13-.121e14*x(4)^6-.329e12*x(4)^12-.102e14*x(4)^5-.116e13*x(4)^11-.681e13*x(4)^4-.296e13*x(4)^10-.355e13*x(4)^3-.583e13*x(4)^9-.134e13*x(4)^2-.332e12*x(4)-.389e11*x(4)*exp(x(4))^10+.569e11*exp(x(4))^10+.348e11*x(4)^2*exp(x(4))^10)/exp(x(4))^10);