z=0.55;s=0
for i=1:6
   if z<=(6-i+1)*0.1&z>(6-i)*0.1
      f=1-i+6*(1-z*1000/600)
else if z<=(6-i)*0.1&z>(6-i-1)*0.1
      f=1+i-6*(1-z*1000/600)
else
      f=0;
    end
   end
    s=s+a(i)*f
end
    pn=a(1)*f(1)+a(2)*f(2)+a(3)*f(3)+a(4)*f(4)+a(5)*f(5)+a(6)*f(6)
%deduced pv
    for i=7:12
   if z<=(6-(i-6)+1)*0.1&z>(6-(i-6))*0.1
      f=1-(i-6)+6*(1-z*1000/600)
else if z<=(6-(i-6))*0.1&z>(6-(i-6)-1)*0.1
      f=1+(i-6)-6*(1-z*1000/600)
else
      f=0
    end
   end
    s=s+a(i)*f
end
    pv=a(7)*f(7)+a(8)*f(8)+a(9)*f(9)+a(10)*f(10)+a(11)*f(11)+a(12)*f(12)