程序是这样的
clear
R=normrnd(0,10,100);
d(1)=150;
for k=1:99
d(k+1)=100+0.7*d(k)+R(k+1);
end
z(1)=150;
for k=1:99
z(k+1)=d(k+1)-R(k+1);
end
for x=1:800
    zz=myfun622(x,z);
    y(x,:)=zz;
end

其中：myfun622是
function F = myfun622(x,z)
h=10;
b=40;
r=30;
w=15;
a=0.3;
jf=quad('t.*normpdf(t,333,10)',0,x);
F = (1-normcdf(x,333,10))-2*a*(h+b)*(1-normcdf(x,333,10))*x*normcdf(x,333,10)+2*a*((h+r)-(h+b)*normcdf(x,333,10))*jf-2*a*(r-b)*z-(h+w)/(h+b);


现在我想加上约束：-2*a*(h+b)*(1-normcdf(x,333,10))*x*normcdf(x,333,10)+2*a*((h+r)-2*a*(r-b)*z这个式子大于0小于1该怎么写，写在哪里呢？

求教~！谢谢~~