function [c, ceq, K1,s] = mycon(x, s)
% 初始化取樣區(qū)間
if isnan(s(1, 1)),
   s = [2 2];
end
% 樣本集
w1x = 1:s(1,1):100;
w1y = 1: s(1, 2):100;
[wx,wy]=meshgrid(w1x,w1y);
%
% 半無限約束
K1=sin(wx*x(1)).*cos(wy*x(2))-1/1000*(wx-50).^2 -...
    sin(wx*x(3))-x(3)+sin(wy*x(2)).*cos(wx*x(1))-...
    1/1000*(wy-50).^2-sin(wy*x(3))-x(3)-1.5;
% 無有限非線性約束
c=[]; ceq=[];
% 網(wǎng)格圖
mesh(K1)
title('Semi-infinite constraint')
drawnow