Problem for doing Monte Carlo Integration

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 :)

 採用された回答

KSSV
KSSV 2021 年 11 月 23 日
編集済み: KSSV 2021 年 11 月 23 日

1 投票

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 件のコメント

Deck Zhan Sim
Deck Zhan Sim 2021 年 11 月 23 日
Thanks for helping, now I can get the answer that I want already.
KSSV
KSSV 2021 年 11 月 23 日
But you need to improve the code. If h = 1, it will throw error.
Thanks is accepting/ voting the answer. :)

サインインしてコメントする。

その他の回答 (0 件)

カテゴリ

ヘルプ センター および File ExchangeProgramming についてさらに検索

製品

リリース

R2021a

質問済み:

2021 年 11 月 23 日

編集済み:

2021 年 11 月 23 日

Community Treasure Hunt

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

Start Hunting!

Translated by