Hi all members,
I am trying to solve a problem based on Monte Carlo simulation and I am a little bit confused.
The problem concerns the computation of the surface area of F(x) = (x^4)*exp(-x) when x belongs to (1,5) through the geometric approach of Monte Carlo simulation. My code for handling this problem is:
>res1=0
>curve((x^4)*exp(-x),xlim=c(1,5),ylim=c(1,5))
>x=runif(1000,1,5)
>y=runif(1000,1,5)
>points(x,y,col=ifelse(yfor (i in 1:100000) {
> x=runif(100000,1,5)
> y=runif(100000,1,5)
> res1
see this
https://www.projectcubicle.com/monte-carlo-simulation-example/
Xuelei Dai, as Anuraj Nayarisseri suggested, try to lower "1" in
>y=runif(1000,1,5)
to something like
>y=runif(1000,1,5) or, to
>y=runif(1000,0.36787944,5)
or even to
>y=runif(1000,0.36787944,4.21121687)
because
when x belongs to the domain (1,5) then the range for y is slightly different from (1,5)
Xuelei Dai Steftcho P. Dokov, pls check the ceiling function and its parameters. i am not sure abt it.
Anuraj Nayarisseri I am very grateful that you took the time to code for my query. My results are in. Thanks a lot
- Badges
- Science topic
- Similar topics
- Physical Sciences
- Materials Science
- Materials
- Inorganic Materials
- Ceramics
- Glass