>> x=[179,178,175,173,172,171,177,179];
C1=150*ones(1,length(x));
for j=1:x(1)-C1(1)
y(j,:)=j*binopdf(j+C1,x,0.85);
end
m=sum(y)
m =
3.1839 2.6272 1.3069 0.7264 0.5186 0.3585 2.1270 3.1839
>>
clear
x=[179;178;175;173;172;171;177;179];
C1=150*ones(length(x(:,1)),1);
n=x(:,1);
for i=1:8
for j=1:n(i)-C1(i)
y=j*binopdf(j+C1,x(i,1),0.85);
end
m(i)=sum(y);
end
m