function f2=CtypeGreyIncidence(f0,f1)  
%%%%%%%%%%%%%C TYPE degree of grey incidence%%%%%%%%%%%%%%%%%%%%%%
%f0=[1000000;1000000;1000000;1000000;1000000;1000000;1000000]/1000;
%f1=[30622.0927;1477.74233;1214.19412;59811.2107;211.272102;604.756576;1485.3764];
global GreyEntropy;
    [n,n1]=size(f0);
     for i=1:n
       d0(i)=(f0(i,1)/f1(i,1)); 
    end
     for i=1:n-1
       d1(i)=((f0(i+1,1)-f0(i,1))/(f1(i+1,1)-f1(i,1))); 
    end
    for i=2:n-1
       d2(i)=(((f0(i+1,1)-2*f0(i,1))+f0(i-1,1))/((f1(i+1,1)-2*f1(i,1))+f1(i-1,1))); 
    end
    D0(1)=(d1(1)+d0(1))/2;D0(n)=d0(n);
    for i=2:n-1
        D0(i)=(d0(i)+d1(i)+d2(i))/3;
    end
    D=0;
    for i=1:n
        D=D+D0(i);
    end
    %f2=1-D/n;
    f2=D/n;
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%為什么f2不行？？？%%%%%%%%%%%%%%%%%%%%%%
     %%%%%%%%%%%%%%%%%%灰熵計算!!(改進并取負）%%%%%%%%%%%%%%%
    ff2=(2+f2)*log(2+f2);
    %f2,ff2
  %%%%%%%%%%%%%%%%%%優度計算  %%%%%%%%%%%%%%%%%%%%%
    ff3=1/(1+((1-f2)/f2))^2;
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if GreyEntropy==1;
    f2=ff2;% 使用灰熵計算;
elseif GreyEntropy==2;
    f2=ff3;
else
    f2=f2;%灰關聯計算。
end
    %%%%%%2007.8.20%%%%%%%%%%%%%%%%%%%%%