Problem for doing Monte Carlo Integration

2 visualizaciones (últimos 30 días)
Deck Zhan Sim
Deck Zhan Sim el 23 de Nov. de 2021
Editada: Deck Zhan Sim el 23 de Nov. de 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 :)

Respuesta aceptada

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

Iniciar sesión para comentar.

Más respuestas (0)

Categorías

Más información sobre Programming en Help Center y File Exchange.

Productos


Versión

R2021a

Community Treasure Hunt

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

Start Hunting!

Translated by