function z=Czfee(x,d)         %計算總費用
%輸入:x 一組標定值組合,d:一組容差組合
%輸出:該組合下的損失費用
dh=0.001;
for i=1:7        %分別求y對x1...x7的偏導數
   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;%求數值導數。
end

mu=F(1,1);%???
sigma=sqrt((g.^2)*((x.*d).^2)');%????
%normcdf(n,mu,sigma)是求y小于n的概率,其中y服從正態分布
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);%e(i)=1表示等級A;e(i)=2表示等級B;e(i)=3表示等級C
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));% 'p*ones(7,1)'  has the same meaning with  'sum(p)'