(Monte Carlo Integration) The average value of a function f over an interval [a, b] is given by Savg = -a f(2) dr. This formula can be used to approximate an integral by approximating the average value of ƒ by sampling the value of f(x) at a large number of randomly chosen points in [a, b] and computing the sample mean of these values. That is, f(r) dx = N i=1

C++ for Engineers and Scientists
4th Edition
ISBN:9781133187844
Author:Bronson, Gary J.
Publisher:Bronson, Gary J.
Chapter5: Repetition Statements
Section5.7: Do While Loops
Problem 6E: (Numerical analysis) Here’s a challenging problem for those who know a little calculus. The...
icon
Related questions
Question

make c+++ code 

(Monte Carlo Integration) The average value of a function f over an interval [a, b] is
given by
1
favg = f(2) dz.
b - a
This formula can be used to approximate an integral by approximating the average
value of ƒ by sampling the value of f(x) at a large number of randomly chosen points
in [a, b] and computing the sample mean of these values. That is,
N
[ f(x) dz =
N
i=1
where the r; are chosen randomly according to a uniform distribution in the interval
[a, b].
Note: While this idea is actually inefficient in 1 dimension, it is a common strategy
for very high dimensional integrals and there are numerous improvements that can be
made to make it even mroe efficient.
mc_int Function:
Input variables:
• an anonymous function representing f
• a scalar representing a
• a scalar representing b
• a scalar representing N
Output variables:
a scalar representing the value of the integral
A possible sample case is:
>> approx_int
approx_int = 0.99821
mc_int(@(x) sin(x), 0, pi/2, 1000)
%3D
Transcribed Image Text:(Monte Carlo Integration) The average value of a function f over an interval [a, b] is given by 1 favg = f(2) dz. b - a This formula can be used to approximate an integral by approximating the average value of ƒ by sampling the value of f(x) at a large number of randomly chosen points in [a, b] and computing the sample mean of these values. That is, N [ f(x) dz = N i=1 where the r; are chosen randomly according to a uniform distribution in the interval [a, b]. Note: While this idea is actually inefficient in 1 dimension, it is a common strategy for very high dimensional integrals and there are numerous improvements that can be made to make it even mroe efficient. mc_int Function: Input variables: • an anonymous function representing f • a scalar representing a • a scalar representing b • a scalar representing N Output variables: a scalar representing the value of the integral A possible sample case is: >> approx_int approx_int = 0.99821 mc_int(@(x) sin(x), 0, pi/2, 1000) %3D
Expert Solution
trending now

Trending now

This is a popular solution!

steps

Step by step

Solved in 2 steps

Blurred answer
Knowledge Booster
Fundamentals of Boolean Algebra and Digital Logics
Learn more about
Need a deep-dive on the concept behind this application? Look no further. Learn more about this topic, computer-science and related others by exploring similar questions and additional content below.
Similar questions
  • SEE MORE QUESTIONS
Recommended textbooks for you
C++ for Engineers and Scientists
C++ for Engineers and Scientists
Computer Science
ISBN:
9781133187844
Author:
Bronson, Gary J.
Publisher:
Course Technology Ptr