function z=Czfee(x,d)
dh=0.001;
d=[0.05,0.05,0.05,0.1,0.1,0.05,0.05]/3;
for i=1:7
   for j=1:7
      x1(:,j)=x(j)*ones(2,1);
   end
   x1(:,i)=linspace(x(i),x(i)+dh,2)';
   F(:,i)=FUN(x1(:,1),x1(:,2),x1(:,3),x1(:,4),x1(:,5),x1(:,6),x1(:,7));
   g(i)=diff(F(:,i),1)/dh;%求數(shù)值導數(shù)。
end
mu=F(1,1);
sigma=sqrt((g.^2)*((x.*d).^2)');% has the same meaning with  ' sigma=sqrt(sum(g.^2 .* (x.*d).^2)) '
y1=normcdf(1.4,mu,sigma)-normcdf(1.2,mu,sigma)+normcdf(1.8,mu,sigma)-normcdf(1.6,mu,sigma);
y2=1-normcdf(1.8,mu,sigma)+normcdf(1.2,mu,sigma);
e=ones(1,7);
for i=1:7
if d(i)==0.1/3 e(i)=1;
elseif d(i)==0.05/3 e(i)=2;
    else e(i)=3;
    end
end
f=[0,25,0;20,50,0;20,50,200;50,100,500;50,0,0;10,25,100;0,25,100];
for i=1:7
    p(i)=f(i,e(i));
end
z=1000*(1000*y1+9000*y2)+1000*(p*ones(7,1));