Problem for doing Monte Carlo Integration

1 view (last 30 days)
Hi all, I have a problem for doing the monte carlo integration.
Below here are all the codes, the %% Monte Carlo is the section for monte carlo integration:
%% Monte Carlo
monte_carlo=0; % Initialize Input
for i=1:h
monte_carlo=c*(sum(y)/h(i));
end
toc % Time Execution Stop
hold on
The outputs is looks something like this
I would be grateful that someone can help me to solve this matter. Thanks :)

Accepted Answer

KSSV
KSSV on 23 Nov 2021
Edited: KSSV on 23 Nov 2021
Your y in the ocde is a anonymous function:
function_handle with value:
@(x)sin(x)
You need to input some value into it....Like
y(1), y(2), y(pi) etc...
You cannot use sum(y) straight away.
May be you have to repalce this line:
monte_carlo=c*(sum(y)/h(i));
with
monte_carlo=c*(sum(y(x(1:i)))/h);
  2 Comments
Deck Zhan Sim
Deck Zhan Sim on 23 Nov 2021
Thanks for helping, now I can get the answer that I want already.
KSSV
KSSV on 23 Nov 2021
But you need to improve the code. If h = 1, it will throw error.
Thanks is accepting/ voting the answer. :)

Sign in to comment.

More Answers (0)

Products


Release

R2021a

Community Treasure Hunt

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

Start Hunting!