I have troubles in for loop with 2 variables.

2 ビュー (過去 30 日間)
Chris
Chris 2022 年 10 月 7 日
コメント済み: Rik 2022 年 10 月 9 日
for n = 2:30
for x=1:n;
p(x) = (factorial(n)./(factorial(x).* factorial(n-x))).*((10/230).^x).*(1-(10/230).^(n-x));
end
end
I want to make like - when n=2: x=1,2. when n=3: x=1,2,3. , ..., when n=30: x=1,2,3, ....,30.
Therefore I think it should be 29 values but it has 30 values.
Would you tell me what do I have to fix?
  2 件のコメント
Torsten
Torsten 2022 年 10 月 7 日
編集済み: Torsten 2022 年 10 月 7 日
I think you will have to include x = 0 to get the complete binomial distribution.
And (factorial(n)./(factorial(x).* factorial(n-x))) can be replaced by nchoosek(n,x).
Rik
Rik 2022 年 10 月 9 日
I recovered the removed content from the Google cache (something which anyone can do). Editing away your question is very rude. Someone spent time reading your question, understanding your issue, figuring out the solution, and writing an answer. Now you repay that kindness by ensuring that the next person with a similar question can't benefit from this answer.

サインインしてコメントする。

採用された回答

Walter Roberson
Walter Roberson 2022 年 10 月 7 日
"Therefore I think it should be 29 values but it has 30 values."
But you said that when n = 30, x=1,2,3,...,30 -- which is 30 values.
Note that you are storing into p(x), where x is the inner loop control variable. For each new n value, it is going to overwrite each p(1:n-1) entries from the previous value of n, and then will write in a new value for p(n) in a location that has not been written into before. Overwriting locations is potentially the correct thing to do, but only in cases where you use the existing value -- for example, y(K) = Y(K) + exp(-t)/factorial(K);
Perhaps you want
for n = 2:30
for x=1:n;
p(n-1, x) = (factorial(n)./(factorial(x).* factorial(n-x))).*((10/230).^x).*(1-(10/230).^(n-x));
end
end
p
p = 29×30
0.0832 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1302 0.0054 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1739 0.0113 0.0003 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.2174 0.0189 0.0008 0.0000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.2609 0.0284 0.0016 0.0001 0.0000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.3043 0.0397 0.0029 0.0001 0.0000 0.0000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.3478 0.0529 0.0046 0.0003 0.0000 0.0000 0.0000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.3913 0.0681 0.0069 0.0005 0.0000 0.0000 0.0000 0.0000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.4348 0.0851 0.0099 0.0008 0.0000 0.0000 0.0000 0.0000 0.0000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.4783 0.1040 0.0136 0.0012 0.0001 0.0000 0.0000 0.0000 0.0000 0.0000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0

その他の回答 (0 件)

カテゴリ

Help Center および File ExchangeLoops and Conditional Statements についてさらに検索

製品

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by