# How could I sums up 1/n for n=1 to 10000?

3 vues (au cours des 30 derniers jours)
Wenchen Liu le 23 Sep 2022
Commenté : Walter Roberson le 23 Sep 2022
I have a code like this, but I can't run it.
n = 10000; % whatever you want
sum_harm = 0;
for i = 1:n
sum_harm = sum_harm + 1/i;
end
sum_harm
sum_harm = 9.7876
##### 0 commentairesAfficher -2 commentaires plus anciensMasquer -2 commentaires plus anciens

Connectez-vous pour commenter.

### Réponses (2)

Chunru le 23 Sep 2022
It can be run as shown below. What problem have you encountered?
n = 10000; % whatever you want
sum_harm = 0;
for i = 1:n
sum_harm = sum_harm + 1/i;
end
sum_harm
sum_harm = 9.7876
##### 3 commentairesAfficher 1 commentaire plus ancienMasquer 1 commentaire plus ancien
Chunru le 23 Sep 2022
There is no variables called "result". You should check "sum_harm" after running your code.
To see the intermediate partial sum, simply remove semicolon at the end of sum statement as shown below:
n = 10; % whatever you want
sum_harm = 0;
for i = 1:n
sum_harm = sum_harm + 1/i % ; suppress display; remove it for display
end
sum_harm = 1
sum_harm = 1.5000
sum_harm = 1.8333
sum_harm = 2.0833
sum_harm = 2.2833
sum_harm = 2.4500
sum_harm = 2.5929
sum_harm = 2.7179
sum_harm = 2.8290
sum_harm = 2.9290
sum_harm
sum_harm = 2.9290
Walter Roberson le 23 Sep 2022
n = 1:10000;
result = cumsum(1./n);
result(end)
ans = 9.7876
plot(n, result)

Connectez-vous pour commenter.

Image Analyst le 23 Sep 2022
Why not try it vectorized:
n = 1 : 10000; % whatever you want
sum_harm = sum(1 ./ n)
sum_harm = 9.7876
What problem did you have in running your code? It ran in MATLAB online. Did you type the m-file name into the command window or click the green run triangle on the tool ribbon? Either should work, but what did you do, if anything? Why do you say you can't run it? Why do you not know how to run the code? Try this link if you don't know how to run your MATLAB program:
##### 0 commentairesAfficher -2 commentaires plus anciensMasquer -2 commentaires plus anciens

Connectez-vous pour commenter.

### Catégories

En savoir plus sur Logical dans Help Center et File Exchange

### Community Treasure Hunt

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

Start Hunting!

Translated by