Problem for doing Monte Carlo Integration

2 vues (au cours des 30 derniers jours)
Deck Zhan Sim
Deck Zhan Sim le 23 Nov 2021
Modifié(e) : Deck Zhan Sim le 23 Nov 2021
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 :)

Réponse acceptée

KSSV
KSSV le 23 Nov 2021
Modifié(e) : KSSV le 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 commentaires
Deck Zhan Sim
Deck Zhan Sim le 23 Nov 2021
Thanks for helping, now I can get the answer that I want already.
KSSV
KSSV le 23 Nov 2021
But you need to improve the code. If h = 1, it will throw error.
Thanks is accepting/ voting the answer. :)

Connectez-vous pour commenter.

Plus de réponses (0)

Catégories

En savoir plus sur Programming dans Help Center et File Exchange

Produits


Version

R2021a

Community Treasure Hunt

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

Start Hunting!

Translated by