MATLAB: Simulation for a triangle

monte carlotriangle

i'm trying to creat a simple script to find the area of a triangle using monte carlo simulation. however, i'm basically stuck. could someone please help me.
any help is appreciated,
thannks

Best Answer

The theory of the Monte Carlo simulation which you describe is, I presume, that if you enclose a triangle completely in a square of known area, the area of the triangle can be estimated by multiplying the area of the square by the fraction of the random points within the square which happen to land inside the triangle. Of course this requires a very large number of random points to obtain a reasonably accurate estimate.
With the triangle you have specified with vertices (0,0), (0,1), (1,0), it is very easy to do this, and in fact you have done just that in your code. Writing x = rand and y = rand will fill the unit square with vertices at (0,0), (1,0), (1,1), and (0,1) uniformly. Since x and y are non-negative, the only test you need for being inside your triangle is for the hypotenuse, whether or not x+y <=1, and that is what your 'count' is recording. That means your code as you have shown it is almost done. Multiply the known area of the square by the fraction count/tries and you have your estimate for the area. You were almost there!
The simulation is a little more challenging with an arbitrary triangle. You have to define a square which is guaranteed to completely enclose it and then generate a large number of random points that are distributed uniformly in the square. Each of these points is tested to see if it does or does not lie within the triangle. The fraction that do lie within it is to be multiplied by the area of the square to give an estimate of the triangle's area.
Related Question