Me parece que has transformado los números racionales en números de tres decimales, eso no está bien y produce errores.
Ya habíamos Visto que k = 1/360
E[X] = [1/360 + 2·2/360 + 3·3/360 +...+ 20·20/360 + 21·5/360+ 22·5/360+...+ 50·5/360] =
Sacamos de factor común 1/360
(1/360)[1+4+9+...+400 + 5(20+21+22+...+50)] =
Lo primero es la suma de los 20 cuadrados primeros, normalmente tendrás que hacerla entera porque no creo que hayas dado esta fórmula:
http://es.easycalculation.com/álgebra/sum-of-consecutivesquare.php
Que dice que la suma de los primeros n cuadrados es:
Suma de los cuadrados consecutivos = n(n+1)(2n+1)/6
Los 20 primeros cuadrados, del 1 al 400, suman 20·21·41/6 = 2870
Luego tienes un 5(20+21+22+ ... + 50)
Esto se calcula usando la formula de la suma de 30 términos de una sucesión aritmética
5·30(20+50)/2 = 5·15·70 = 5250
En resumen te queda:
E[X] = (1/360) (2870+5250) = (1/360)8120 = 22,555...
Que puedes ver es bastante distinto del 22,89 que te daba. Habías metido muchos redondeos y hay que procurar operar con fracciones y no con números decimales
La Varianza se calcula como bien has hecho mediante:
V[X] = E[X^2] - (E[X])^2
Y para calcular esto de nuevo hay que usar fracciones y al final fórmulas no habituales para evitar hacer muchas sumas.
E[X^2] = (1/360)[1+ 2^3 + 3^3 + ...+ 20^3 + 5(21^2 + 21^2 + ....+ 50^2)]
La suma de los primeros cubos también tiene una fórmula:
http://es.wikipedia.org/wiki/Cubo_(aritm%C3%A9tica)
Suma de los n primeros cubos = [n(n+1)/2]^2
En otros sitios dice que es la suma de los n primeros números elevada al cuadrado. Son exactamente lo mismo.
Luego la suma de los 20 cubos primeros
Suma 20 cubos primeros = (20·21/2)^2 = 210^2 = 44100
También tendrías que hallar la suma de los cuadrados del 21 al 50. Eso se puede hacer como la suma de los 50 cuadrados y restarle la suma de los 20
Suma cuadrados de 21 al 50 = 50·51·101/6 - 20·21·41/6 =
42925 - 2870 = 40055
Y esos cuadrados tenían un 5 multiplicándolos
5 · 40055 = 200275
E[Y^2] = (1/360)(44100 + 200275) = 244375/360 = 678,8194444
V[Y] = E[Y^2] - (E[Y])^2 = 678,8194444 - 22,555...^2 =
678,8194444 - 508,7530864 = 170,066358
Y finalmente:
s = sqrt(170,066358) = 13,04094928
Como puedes ver de nuevo hay diferencia apreciables, también motivadas por el error en el calculo de la media que se arrastraba.
Y eso es todo, no sé si estos métodos son demasiado avanzados, pero sirven para comprobar la exactitud de las cuentas Y si se saben manejar nos evitan hacer muchas.