function z=Cfee(d)
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*(p*ones(7,1));