zero value in beronoulli trials !!!

Hi sir,I want to know,why some values in my program equal to zero

This is my program

N=10;

P=0.5;

for K=0:N;

 format long e
    M =exp((P)^K*((1-P)^N-K)*(factorial(N)/factorial(K)*factorial(N-K)))

display(M)

end

Thanks in advance

Answer 1

[Andrew Newell]

The correct expression is

M =P.^K.*(1-P).^(N-K).*factorial(N)./(factorial(K).*factorial(N-K))

If you have the Statistics Toolbox, you can confirm that this gives the same output as

binopdf(0:N,N,P)

[reem]

Thank you so so much Oleg and Andrew

I will not forget your help forever

Kind regards

Answer 2

[Oleg Komarov]

Some other parenthesis mistakes corrected:

N=10;
P=0.5;
K=0:N;
M =exp(P.^K.*(1-P).^(N-K).*factorial(N)./(factorial(K).*factorial(N-K)))

[reem]

Thank you so much Oleg

I hope all best thing for you

Answer 3

[bym]

I think you want N!/(K!*(N-K)!)<--note parentheses

[reem]

Yes I want that,so what is the problem here in my program