有约束的最小最大值matlab程序

function []=fminemaxc（ab）
%调用的格式为fminemax1（ab），a=[0.30.69]b=[9.9111.17]
%由于全局变量只有两个xy所以这里就引用两个二维矩阵作为参变量
if （a（1）<0.3）|（a（2）>0.69）|（b（1）<9.91）|（b（2）>11.17）
    disp（‘wrong input‘）;
else
    runfun（ab）;
end
%以下为子函数
function runfun（ab）
 y1=9.91;
 y2=10.54;
 y3=11.17;
 Rg=10;
r4=0.0033;
Dis42=0.69;
xtem=size（a（1）:0.01:a（2））;
ytem=size（b（1）:0.01:b（2））;
Utem=zeros（6xtem（2）*ytem（2））;
i=1;
for xi=a（1）:0.01:a（2）
    for yj=b（1）:0.01:b（2）
        
        D41=sqrt（abs（Dis42-xi）^2+（yj+y1）^2）;
        D42=sqrt（abs（Dis42-xi）^2+（yj+y2）^2）;
        D43=sqrt（abs（Dis42-xi）^2+（yj+y3）^2）;
        
        d41=sqrt（abs（Dis42-xi）^2+abs（yj-y1）^2）;
        d42=sqrt（abs（Dis42-xi）^2+abs（yj-y2）^2）;
        d43=sqrt（abs（Dis42-xi）^2+abs（yj-y3）^2）;
 
        Z44=60*log（2*yj/r4）;
        Z41=60*log（D41/d41）;