Calculaing the table of values for function defined as an infinie series
1 vue (au cours des 30 derniers jours)
Afficher commentaires plus anciens
Chuzymatics Chuzymatics
le 19 Juil 2014
Commenté : Chuzymatics Chuzymatics
le 21 Juil 2014
I need this code below to caluate the sum of the first 50 terms of the series :t to the power n for t value. That is, a kind of table of values for the series for t=1:10 It gives the error message: ??? Error using ==> power Matrix dimensions must agree. Please someone help % Calculating the sum of an infinite series
t =1:10;
n=0:50;
E(1,:)= sum(t.^n)
0 commentaires
Réponse acceptée
Jos
le 19 Juil 2014
Modifié(e) : Jos
le 19 Juil 2014
you can use a for loop to do this
t = 1:10;
n = 0:50;
for i=1:length(t)
E(1,i)=sum(t(i).^n);
end
Cheers
Joe
4 commentaires
Jos
le 20 Juil 2014
Both the method above and the one below by Roger (a lot more elegant by the way) produce exactly what you would expect
E:
51
2.25e+15
1.08e+24
1.69e+30
1.11e+35
9.70e+38
2.10e+42
1.63e+45
5.80e+47
1.11e+50
plotting E in the way you try to do will basically only show the last point as it is a lot bigger than the other points, try plotting the y axis on a logarithmic scale using semilogy(t,E), you'll get this graph
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/179751/image.jpeg)
Plus de réponses (1)
Roger Stafford
le 19 Juil 2014
Also you can use this:
E = sum(bsxfun(@power,1:10,(0:50)'),1);
4 commentaires
Roger Stafford
le 20 Juil 2014
Modifié(e) : Roger Stafford
le 20 Juil 2014
A Side Note: The value of your E(10) would be 1111...11, a string of fifty-one 1's in decimal, if it were computed exactly. This reminds me of one of the Fifty-Eighth Putnam Exam questions which posed the intriguing problem: "Let N be the positive integer with 1998 decimal digits, all of them 1; that is, N = 1111...11. Find the thousandth digit after the decimal point of the square root of N." (Surprisingly, it is solvable without using a computer.)
Voir également
Catégories
En savoir plus sur 2-D and 3-D Plots 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!