Volume of n-sphere? What is wrong with my version?

4 views (last 30 days)
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.

Accepted Answer

the cyclist
the cyclist on 11 Nov 2013
Edited: the cyclist on 11 Nov 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.)

More Answers (0)

Categories

Find more on Image Processing Toolbox in Help Center and File Exchange

Products

Community Treasure Hunt

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

Start Hunting!