exc10f2.m

来自「薛定宇高等应用数学问题的MATLAB求解例子源码」· M 代码 · 共 10 行

   function [sol,y]=exc10f2(sol,opts)
   x=sol(1:6);
   y=-1/(2*cos(x(6)))*(x(1)*x(2)*(1+x(5))+x(3)*x(4)*(1+31.5/x(5)));
   c=[-0.003079*x(1)^3*x(2)^3*x(5)+cos(x(6))^3;
      -0.1017*x(3)^3*x(4)^3+x(5)^2*cos(x(6))^3;
      -0.09939*(1+x(5))*x(1)^3*x(2)^2+cos(x(6))^2;
      -0.1076*(31.5+x(5))*x(3)^3*x(4)^2+x(5)^2*cos(x(6))^2;
      -x(3)*x(4)*(x(5)+31.5)+x(5)*(2*(x(1)+5)*cos(x(6))+x(1)*x(2)*x(5))];
   if any(c>=0), y=-1e10; end