Volume of n-sphere? What is wrong with my version?
1 visualización (últimos 30 días)
Mostrar comentarios más antiguos
Steve
el 11 de Nov. de 2013
Editada: the cyclist
el 11 de Nov. de 2013
Hello Experts,
I have to calculate the volume of 10-dim sphere.
Here what I did:
N = 100000;
RandVars = -1 + 2*rand(10,N);
V = zeros(N,1);
for i = 1:N
if (sum(RandVars(:,i).^2) < 1)
V(i) = 1;
end
end
% Mean values using the Monte-Carlo method:
M = (1/N)*sum(V,1);
% Error estimation:
C = V - M*ones(N,1);
Error = sqrt((1/(N-1))*sum(C.^2,2));
Please have a look what is wrong as I get M = 0.002. Thanks a lot in advance.
2 comentarios
Respuesta aceptada
the cyclist
el 11 de Nov. de 2013
Editada: the cyclist
el 11 de Nov. de 2013
You forgot to take into account that your bounding hypercube has volume 2^N, not 1, where N is the number of dimensions. (N=10 in your specific case.)
0 comentarios
Más respuestas (0)
Ver también
Categorías
Más información sobre Monte-Carlo en Help Center y File Exchange.
Productos
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!