Problem 44694. Monte Carlo integration: area of a polygon

Created by David Verrelli

Solve Later

{"relationships":[{"relationshipType":"http://schemas.mathworks.com/matlab/code/2013/relationships/document","targetMode":"","relationshipId":"rId1","target":"/matlab/document.xml"},{"relationshipType":"http://schemas.mathworks.com/matlab/code/2013/relationships/output","targetMode":"","relationshipId":"rId2","target":"/matlab/output.xml"}],"parts":[{"partUri":"/matlab/document.xml","relationship":[{"relationshipType":"http://schemas.mathworks.com/matlab/code/2013/relationships/image","targetMode":"","relationshipId":"rId1","target":"/media/image1.png"},{"relationshipType":"http://schemas.mathworks.com/matlab/code/2013/relationships/image","targetMode":"","relationshipId":"rId2","target":"/media/image2.png"},{"relationshipType":"http://schemas.mathworks.com/matlab/code/2013/relationships/image","targetMode":"","relationshipId":"rId3","target":"/media/image3.png"}],"contentType":"application/vnd.mathworks.matlab.code.document+xml","content":"<?xml version=\"1.0\" encoding=\"UTF-8\"?>\n<w:document xmlns:w=\"http://schemas.openxmlformats.org/wordprocessingml/2006/main\"><w:body><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:r><w:t>The area of a polygon (or any plane shape) can be evaluated by</w:t></w:r><w:r><w:t> </w:t></w:r><w:hyperlink w:docLocation=\"https://www.mathworks.com/matlabcentral/cody/problems/179\"><w:r><w:t>Monte Carlo integration</w:t></w:r></w:hyperlink><w:r><w:t>. The process involves 'random' sampling of (x,y) points in the xy plane, and testing whether they fall inside or outside the shape. That is illustrated schematically below for a circle.</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:customXml w:element=\"image\"><w:customXmlPr><w:attr w:name=\"height\" w:val=\"-1\"/><w:attr w:name=\"width\" w:val=\"-1\"/><w:attr w:name=\"relationshipId\" w:val=\"rId1\"/></w:customXmlPr></w:customXml></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:r><w:t>^</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:r><w:rPr><w:b/></w:rPr><w:t>Steps:</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"ListParagraph\"/><w:numPr><w:numId w:val=\"2\"/></w:numPr></w:pPr><w:r><w:t>Define a 2D region within which to sample points. In the present problem you must choose the rectangular box (aligned to the axes) that bounds the specified polygon.</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"ListParagraph\"/><w:numPr><w:numId w:val=\"2\"/></w:numPr></w:pPr><w:r><w:t>By a 'random' process generate coordinates of a point within the bounding box. In the present problem a basic scheme using the uniform pseudorandom distribution available in MATLAB must be employed. </w:t></w:r><w:r><w:t> </w:t></w:r><w:r><w:rPr><w:i/></w:rPr><w:t>Other schemes in use include quasi-random sampling (for greater uniformity of coverage) and stratified sampling (with more attention near the edges of the polygon).</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"ListParagraph\"/><w:numPr><w:numId w:val=\"2\"/></w:numPr></w:pPr><w:r><w:t>Determine</w:t></w:r><w:r><w:t> </w:t></w:r><w:hyperlink w:docLocation=\"https://www.mathworks.com/matlabcentral/cody/problems/198\"><w:r><w:t>whether the sampled point is inside or outside the polygon</w:t></w:r></w:hyperlink><w:r><w:t>. </w:t></w:r><w:r><w:t> </w:t></w:r><w:r><w:rPr><w:i/></w:rPr><w:t>Due to working in double precision, it is extremely unlikely to sample a point falling exactly on the edge of the polygon, and consequently if it does happen it doesn't really matter whether it's counted as inside or outside or null.</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"ListParagraph\"/><w:numPr><w:numId w:val=\"2\"/></w:numPr></w:pPr><w:r><w:t>Repeat steps 2–3</w:t></w:r><w:r><w:t> </w:t></w:r><w:r><w:rPr><w:rFonts w:cs=\"monospace\"/></w:rPr><w:t>N</w:t></w:r><w:r><w:t> times.</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"ListParagraph\"/><w:numPr><w:numId w:val=\"2\"/></w:numPr></w:pPr><w:r><w:t>Based upon the proportion of sampled points falling within the polygon, report the approximate area of the polygon.</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:r><w:t>Inputs to your function will be</w:t></w:r><w:r><w:t> </w:t></w:r><w:r><w:rPr><w:rFonts w:cs=\"monospace\"/></w:rPr><w:t>N</w:t></w:r><w:r><w:t>, the number of points to sample, and</w:t></w:r><w:r><w:t> </w:t></w:r><w:r><w:rPr><w:rFonts w:cs=\"monospace\"/></w:rPr><w:t>polygonX</w:t></w:r><w:r><w:t> &amp;</w:t></w:r><w:r><w:t> </w:t></w:r><w:r><w:rPr><w:rFonts w:cs=\"monospace\"/></w:rPr><w:t>polygonY</w:t></w:r><w:r><w:t>, which together constitute</w:t></w:r><w:r><w:t> </w:t></w:r><w:hyperlink w:docLocation=\"https://www.mathworks.com/matlabcentral/cody/problems/194\"><w:r><w:t>an ordered list of polygon vertices</w:t></w:r></w:hyperlink><w:r><w:t>. </w:t></w:r><w:r><w:t> </w:t></w:r><w:r><w:rPr><w:rFonts w:cs=\"monospace\"/></w:rPr><w:t>N</w:t></w:r><w:r><w:t> will always be a positive integer. </w:t></w:r><w:r><w:t> </w:t></w:r><w:r><w:rPr><w:rFonts w:cs=\"monospace\"/></w:rPr><w:t>polygonX</w:t></w:r><w:r><w:t> &amp;</w:t></w:r><w:r><w:t> </w:t></w:r><w:r><w:rPr><w:rFonts w:cs=\"monospace\"/></w:rPr><w:t>polygonY</w:t></w:r><w:r><w:t> will always describe a single, continuous outline; however, the polygon may be self-intersecting.</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:r><w:t>For polygons that are self-intersecting, you must find the area of the enclosed point set. In the case of the cross-quadrilateral, it is treated as two simple triangles (cf. three triangles in the first figure below), and the clockwise/counterclockwise ordering of points is irrelevant. A self-intersecting pentagram (left &amp; middle of the second figure below) will therefore have the same area as a concave decagon (right).</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:customXml w:element=\"image\"><w:customXmlPr><w:attr w:name=\"height\" w:val=\"-1\"/><w:attr w:name=\"width\" w:val=\"-1\"/><w:attr w:name=\"relationshipId\" w:val=\"rId2\"/></w:customXmlPr></w:customXml></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:customXml w:element=\"image\"><w:customXmlPr><w:attr w:name=\"height\" w:val=\"-1\"/><w:attr w:name=\"width\" w:val=\"-1\"/><w:attr w:name=\"relationshipId\" w:val=\"rId3\"/></w:customXmlPr></w:customXml></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:r><w:t>^</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:r><w:t>EXAMPLE</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"code\"/></w:pPr><w:r><w:t><![CDATA[ % Input\n N = 1000\n polygonX = [1 0 -1 0]\n polygonY = [0 1 0 -1]\n % Output\n A = 2.036]]></w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:r><w:t>Explanation: the above polygon is a square of side-length √2 centred on the origin, but rotated by 45°. The exact area is hence 2. As the value of</w:t></w:r><w:r><w:t> </w:t></w:r><w:r><w:rPr><w:rFonts w:cs=\"monospace\"/></w:rPr><w:t>N</w:t></w:r><w:r><w:t> is moderate, the estimate by Monte Carlo integration is moderately accurate (an error of 0.036 in this example, corresponding to 509 of the 1000 sampled points falling within the polygon). Of course, if the code were run again then a slightly different value of</w:t></w:r><w:r><w:t> </w:t></w:r><w:r><w:rPr><w:rFonts w:cs=\"monospace\"/></w:rPr><w:t>A</w:t></w:r><w:r><w:t> would probably be returned, such as</w:t></w:r><w:r><w:t> </w:t></w:r><w:r><w:rPr><w:rFonts w:cs=\"monospace\"/></w:rPr><w:t>A</w:t></w:r><w:r><w:t>=1.992 (corresponding to 498 of the 1000 sampled points falling within the polygon), due to the random qualities of the technique.</w:t></w:r></w:p></w:body></w:document>"},{"partUri":"/matlab/output.xml","contentType":"text/xml","content":"<?xml version=\"1.0\" encoding=\"UTF-8\" standalone=\"no\" ?><embeddedOutputs><metaData><evaluationState>manual</evaluationState><layoutState>code</layoutState><outputStatus>ready</outputStatus></metaData><outputArray type=\"array\"/><regionArray type=\"array\"/></embeddedOutputs>"},{"partUri":"/media/image1.png","contentType":"image/png","content":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAPcAAADwCAYAAADcifLrAAAjU0lEQVR42u1dB5gV1dn+7tBdK01AsFIEESKIERULttgNiAZ77GBAorFrRESj2DU2NJYoEDF2fzuCigYj2FEjIghSpHdEYOf/3nvO6t3Zmd25e6edme99nvdhmb17p5zzzjnnO18hEggEAoFAIBAIBAKBQCAQCAQCgUAQPjozP2b+II9CIEgPTmLOYY4RcQsE6cJpzFbM00XcAkE6UZO4t2F2F8bOXswDNI924cH6d/vLs0oEO5kg7huZS5jThaHwO+Zs5jzmIuYy5irmGuY65npmOdMu4EYXOj+zgfkzc63+vuXMxcz5ur1nyrMPjXi+n5gi7mtkglMy6jA7MI9n3sB8SYv6J+bXzFeY9zIvZv6BeSCzK7M1s1GR56rLbMHchbkv8/fMPzFvZT6jjah4iSxlTmDexTxTjzgNpalKxm9F3OnG5sxDmEOZrzNXML9nPs8cxjyOuRMzF+ZF8BC+m53LPc2cZFvWffz/pgW/bqGv8RLmKOYXeqT/gHkHs59elglE3JkfmfdiDmdO0dPq9/Vo2UcLKVKwkNuxqFcz7QJ+wsfrVfNnm+pZw1/1bGKZXjrcp9f2ZdLUyRd3Sz3l+7Ne77XWbCji9o3NmCeS2k5cpBt0BPOgWkylgxe3Zd3iELYi0aFFfI3F7Ma8jPmWXrtjJnKB7i+CBIp7vsMAU8FjRdzVAiNbf712RUd/kXlqHCOzD3E/5iHuU0r4WozchzMfImWgm6iFLtP3BE7La4KIWwGGqSdIWZ7HM89hNk7yBUPELuJey8e3D+gUDfRUHTOXlfq59NfHRdwi7kQD4h3C/FITS5hmplw8izjHo/dtLOhyLezltrLWhwEsQ/DdbzJ/ZN5MandAxC3iThTaM+8nZUx6Qo/axoIFvQ1zH1tZ8KNAW+ZNWuQwyh0o4hZxxwlsSR2l15Dwub+UKm8bCYpHPW2P+Iw5Ta/NG4m4RdxR4gAt6pm6A25aw0h4AE93r+V/L7ITaEhLIGBxhyPNe/oZw2GmrohbxB0mMFL/h/kNKYcNy8cU9xKHgWqRrZxRBP7QU6/LZ+kXqcP4Zndncp+zr2W2E3GLuIvFrnotiP39weTTussibsxiXlfFAm1Zo0SzReMIUs4+X9Ev26/2Scz1TFtzNXNfEbeI2w+whob/9o9a1JsU88fc23q47hvnch+LVmtt54DIp/KPPJpvWFog7Ar+V8Qt4q4OMOxgG2sB826q5f40rM0s5DUuI/fDotOSwOvvXS5zETa4PL7LsjcRcScbPfSDHkcBxNhybxvAgt5YIO7ZtrhhBoSNP7iIe3wMoh7IXMjcyJzI3FnEnSxgdP4nKTfbfoE2PVE35hDmGcwtRZSBPdkjmWt/FfYKXn+P7BvxNRzn8oL5rshRXMQd8poHkUxjKcVbVdzrmmFJwLOHr5jv2coV1PS76sD8E9Hq84l2Q0QdEoWcFeH5/+2xPOgt4o62IzRiNq+8dstfK0brvpRicG+zWNDvO9b/5bZKs5Qm7KKFghf1VhE82TEe4u4l4o6ma7OI7duZP+sHP4nocoxanzNfIIP8v4sQMyz1bzHnMZ/n///e1XJvWWNT+C7DS3socy4p63qYT/poF2F/w2wo4o6mq19etQG+3kBUHy6jVgqF3Z6Fu6qK84z7ttwrKZ6wIGsMEotcTyo5RlhPHAa1OXrweEcMatF29/96TJ12SeUU3LKu8xDyUpdY7XMp3cAybAKphBFNEmzzEXHXsruP9xD39ikV980eSRdw/Ptf1tuWdaedwpmLC+D8glkaEkz2FHGnStzvXO0i7LeYuVSKm2jfgpjswqQLOzLrMDvakRibEgckhkB6q1NE3OkQNwwqC4jGjORuzw/QnsV8gJnq8EwW78CCafgcW0VZCZQ1HbnCh1HImWRF3OFiKKlooi5Z7MV6lEbihZxouhLgsPQ2qbx2m4i4zRM30vZ8y2wnfVngAiRuxDbou6Sy0oq4DRA3Rqm/k6qS0Uz6sLGzjr62ZT3CHGmHl7YK22P/YE6ieG0QIm4fQDTXv0gF90syfHOFfZHDGLgxZFfZa0nFibcWcdco7g4jtFN/1wjPC48g5AR/NSHrKEFt7QS53AqXrbwpIZ8aW2UzKJ4sOaaI++rniNatK9hyeoZZP6SO0BC+0cuIDipT6yd4WzUSiRgt7mYeDjhLIjj9laQCiNqIuKs2TQuin3922VO+JIRO0Jkb/LuKxp+Xy62eQ7SbyCNxYu2Q92sn2t333+Ry37qI+/8iumQk6UABxu1E3JWb5UgPT7DAG4Yb+12XDvBuQaeqZ1vW9XxsBnMW/3y3LdP1qIV9RaVEFZY1xvaXVHJv/vyCgnadFnFSSazBkQSzpYj712bp5iHukQF3mrrc4OtdxL2+oiIld6RhLhFQI0VykQl7V0cGmgoX2JN9/n0TZj/mUXY8xlHsuGCdXybi1uMp0efTHcJGpoyOQZ9pZS63zEXccwtG9jkuv19nS12qqMR9rkeI6QOG3AK2VJGlFvXR64i489jyFqKxCNR4mfkYMwzvsJMGEi2p4j9dkIHDI8RxvUzNIxP34R7BK1cadBuYBSKa7CERt0LY+9x7k3L+R16y3/FIMDpPosMqdS7Lutelcz0jsotM3NjSmuh4/vBxb2zYraBe2qfMS0J6Us2Jxg0nGjJLGaSzK25YMDH17uejc22WN+CotfkGWFvtBHms8bUcnPe+UmGXXVIq8EYYqfnZv5SvHkq0taG3AmMeUnEF7Ehj76GzplYsYRcze2ZR3HBSmVzsG5SfWH1b/W2SOv1gF1tA1ipbmgaIjoUYZJlhe4qL8fmzLIp7NKnUw6aPZrD0r3ZZMkwS/SQeZ5PaIgsgPbVdh7nBRdzlzAZZEjcK2n9BKfAX59Zr5eGB9aNoxwgg0ATekAFksMknUnTJeZ6dkRsVQOB6uGtK1qG5vFNNVXE/L7oxAhhgkDX3ogB6Qx/H6I2KJcdnRdywVOJN1idNvYNb8VC78j79DFvizk3CDszFpLawSu0NexJ9/QTRAwv4572TfuNBivvhNKyzPQTelHm89sCSvXfzcBpzGnPTAL4rc/vcmJ4gBG8L6UeRv3jaMK/OZ1El2keeiCeQO+BeVYDA7qwq2oi4awJyTaPQfW/pP5ELu6sjvzlKDg2SJ+OKZkTnLCXasEavmVcxz0yCuOHQgVIy2JxHHeqnqHonj6mq7Svxi5DEjWu5SfpODOKGs09VY99PsnRwfVrdico3Oqzd+H+PuMX9Gqka1Ej1isAO+NC+Uc3nMZKeQSoVTQW3DkHcx+q1jCRdiEfc0zx8wjvK06nytK70iIK8Jk5x76BH3sIG66iPtfX4m5/IX6K6UsQN54A5lL7Kk+Z0V8sa6yLu5fD0k6dT5WkN9BD3oDjFDd/spS7Hcay/y/EyLXzEQ8OJfiYpj7FWAYv7dlJhdlle8yLJxFAW1CfMj/j/l9jhhxwWnr+tI1x2nd8Y7Ay2VnOHb3iFf3iLOMXNb5x8Tm8ncMztrdNSr83vYyLpISyo7zHhB1s3IHF3JuWs0iLT3cWybok7TJLPtznzJOYAO56EgSa12M66RveXRC8v4lXlsFp8SaDiHuAh7unk3zLaRo/m+zmOP6un1pMLeFgN34XA+InMCzM+atfx8DufJyIK5XkfmJ8lEZ0WUBIP5PBDOHLTOMXdl7nM5fgK8hFO6fi8c8p2FxPZNroXsKYYXpwTeaPr+RIAUS/mMXZMeaZ1pNm2dtVZS6nf20iHpzrFvVJKAgU+Q7rd8Yw/toNxSHlQz3BjE/e2etTtVHAMMcXl+ndOdNWCrevyHT1LnJbX1zOGo310/jJHAgBUrjwxYmEPzp9XnX+ZrZxtgvv+XO51l9RET2ZkJG3H9/oQ3/Mr/O/f7JDK/PD37uaxI3B1AF/fUi8vO8UlbgCRLeP0RWA7bIKeUlfgPE2gub7g+0klTMD6eDzzQ6o6ohQrbpzjLZ9v2+s8RrWtIup8vTxK4+4U4Dm2cbzAXreLn+aZKOw2dtW8eJPCMCbyd57ikd8tqJfoVaR8SGITNwQxWk+tlzMfp8quns9pVmB3/QKo+PwYKn2fG44Rc8lnTmtugPEeb9xIPNm48Ud4nP/8EDrg9nb0yfHjnCbf4PFsDw7h2XYLceQGsLsEO8kecYk7LBQj7ov1DMJvBxjl0Sg7R9QBh3qc/3RZxZb8bB+N8tny+f7hONc3drBxDChu8HJWxY31FNxefcdp61zYqx1Tqcj2xfNrwlxulTOdMnJsizxLfrZnemSrbR/S+XLMPtx/hmPmFZAxrRBI7TWb/Dl+pU7cQxzTfr+N0oEb5H4kN9CNUifiTtiDz/0WtqeYL0gcdmDP1eJ2fdjhz36O4bc12OfMNFXihoX8B1JZVgSCQpF3QrleO7pSPmGioV57d86SuE9iviNdORKxbMk8gzkEhqSU32v9fG1v2GYs69aEGCShhYezIm5sncE//RiRXuidHdtLswumuhvsX7c303avyE33smPNvsh2j3+IEnDgwjZyqyyIG+6qcH2tI/ILucNXXsNWcHVYziExi7unx971iARcHgKuhmdB3EjEMESkF0GHh0ul+/ZSjxSKu7+HuMcm4PKw5sbOUP00ixtTEzjASF60aEbuJz22l5qkUNztXMs6JycYCYFR/dMs7ssppdlME9rhO1Spdkp0WYrvd7BD4C8lKMnE6aTcvVMpblRoQA7yXhkXXLd8UIRlXRuWg4bjfC2YV2jrce8MPN+2zFN15GCSIumQNgyGtXZpFDc8dVBnKZdhYfd3jCwwbu1v0PU35JfEMB3Y8ipyrsv8qCjAsHZtGsV9v56WZ1XY2KpZ4LIGnmzQGn60y5r2CNGsb0DE36ZN3IgDR/nTHTMs7sYeBQBXGHL9LTyu/03RrG9g1opCG93SJO5DmFOy3rIshJku4njbEHF39hD3p6LZonATVc3Jb7S4Ufb0ssyLm+hwR1QZPKi6G3Lt9Vyrk1rWLaLXooDcBbOosu3JWHHDSo4N/PbSrnmRtGaerZPyNTXs2n/Lgv6+QNyIiiuTVi16aj6TVDJF48W9J/NradPUvJwa5EUuL+tS8HdSqZiMFzdM/3dIewoEv+BIUnn/jRf3B8xDpT0FPmcGSBD5NHM+86uUprDCUmYV/Zry20hxI2vqSlJB6wJBTcKui31/F6Pd6KBzxCcAKLx5vMnixsW/Jt1W4FPcu3tst8FZ5oqU3e4lpDzWjBX33RRxnSuB0eLex1Pcudz3KbtdCPork8X9EVWtJSYQeIm7keteuuIPKbvdenrd3cxEcSPbxxpS0TACgV+BI331Cpd1920pvF1U7TnKRHEfyJwUQWfYxJaUTWkTOPzYJzA3Mn9mYf/TVtVp0obrtF6MEzcMIHeG2AHa52tJVeS45je7iDx1Im9o+6j8ajCO0qO3ceJG7qozQmp0CwELLhZVyc0mMAlIvbyUlBenUeKGy+nuIYm7o4fR5X3pLwLDgOwsR5skbqRxXUchGdPyqXTcxT1B+orAMCDk90KTxP0AhRgsoqflE12m5WdKXxEYBgSR3GOSuJ9hPh3mSVjIzW3LGsOiXplPgkD0J+knAgMxgFSpX2PEPUH/KxAIqsdBzM9MEjcu9DxpN4GgRrQllV/QGHHPZP5O2k0gqBEomrBBj95GiHsZs4O0m0DgC/Pp1yCSRAOZHTeS+JQLBH6BKfl0Ey70LuZP0l4CgW+8xTQipPVBUl43AoHAH55lGhHS+oQpFyoQJASP6HV34gHnlWnSXgKBbyA78MI4LwAZI8bqNwyKDDyljzkBb5vPpL0EAt8YRio6LDYg0SGKh+/C7Egqe+MbLp/DscnSXgKBbyD3wfK4Tr4D5d2586KuQEd9rK3js7D8fSDtJRD4xsXM2Kq99vOYNuBYf8cxhLBJXLVA4B8XkMrvHwsGknvRcBwb5Dg2kfmutJdA4Bvnk8qEGgsGeIh7uou4P9XHLy1gp9qclOf8rXRFzNNNq4gpEBSBs0llCo4FfUn5izuxQk/ZCzGVOZdUwoYK7lELYe9n53JLC5IxzOdjv5F+IEghzopT3NsqvVUagbswy/XvAp+Ws5inuaRRmij9QJBCDIxzWg68QGorDALHdtgEUm5zTpRsUOO3yOYeOdLWSD8wG9y2OWYf27Ku5X/72qoYfdYRq0EN2Io5Wk/FsSf3OHMLl88FshXGQp7nIu5PpB8YLWzkvnvB0aYv4njGH81fKMatsGIQiBMLN/ip3PDlBZ1gPR87XCRitLiP96jgeULGH02sTizF4BUKKGVM3qhmWbcyR9gh5UAXRChuTMXdxG1Z12X80QylmN1P/QKZT/8nXVng8rI+QUZuV9zOXGTChY5izpauLHARdx0W83iHuMdLjTd6mAwJ+XyIuVi6ssBD4HXzTkmWNVw7J9WVp5IPk55jwoXezVwr7SUQ+MabZEiapZtJObfUkzYTCHwBu0vfmXChSG2MDfmdpM0SNR0+jHkR8ziZCicOSEv2tSninsU8UNosIcK2rEed1VBF4IkBjInryaCiBLjQs6TdEjFi7+ex/SQVUZOB7UjtcRtTTugdUjW6BfGL+wIPx5F7UnBv8FPfxvCttP1JRVIaI24EmYwSaSVCAId7jNwXGH5fffk+5ur7WWKbW3jyj6Ss5caI+yGSDKjJGd1yuecd4v6Ij29i8D215Xv4yXFP5Xx8bwNv5zbmSJPEfQNzHakKhoL4xWAxT+ap+F3870BmQ8Pv53yPpcZNBt4OAq0uNUnc15Dat+si0hKEIO4/eiw1/mrg7aAGQB/TxP0c8+SQG7kOs51mHen2mRF3U5d4/5V21TTbScfWpHxCepom7uH657AauA036OSCxp2CY9L1MyPw9jrxA0Q+ztCQYPiCILHJb00T9xGkUjGF07jI4FF1avaSdHuBQbiKVJ0w48SNtEyrqcCoxm/XHbVBZwBz+xLe2nVdrKV2/pj4tAvMARKb9DNR3MBUfeEQZG8W36oCIa7gY/uUMHIvcBH3Qkm4JzAEsBEhtVIrU8X9IPMiLcYvXMQ4pYTR+zIXa+nl0mcEhqArc4b+2Uhxn8F8mkXXwJHssDDpYZ1aijuXn97ncpPyVNN9GbUFpgBVfEabLO7WpCqV1GMBznAR99fSxukAv1iPQKLDvC870ZbyRGoEtopPN1ncwOfMXtzgx7CY1zmMX4dJG6dA2MhQW/mlPYvbtpk8GU80IJWnvIXp4h5Byh0Vb/cu3BFu1EnyJvPPYyVlsfEj9k6uSy7LukWejid6Mz8q+L+x4saNfPxLZ8jlnnN0hJ9tbVEXGCnuozxKP70qT8cTN1cMeKaLG1MQuNi14I7QycPhf6y0t9Ej90YZuYsCIib3T4O4gSeZg7gj7Ovxlh8v7W30mnuoY2r+Fbd1E3kyrkARTeQor5MWcSPqZaKu3rnYZX/6KmnzRIzCVt7qTTTELrKmuq2MptcgfZPJseIhPt26zEuJZs8l+h/TPikt4kb8MLbEtuOGP5IFvaxA3C+bHl+cEmHX57Z4o1LyA8v6Wwbue7No/CPse5i2g+emQdzAE8wL9QPdKj9FJ9pVHE8S08kHe8RId0/p/e7N9/dNgSt0iGmn7DLmBhdxf5UWcf+eAqjbLQht3fyIh7jPTqGwG7vGJoTmd2E3dRE2OCct4kaebF5rUDeRUiI7/BUe4t4vhffax2PX5uHwzlo+2UXcD6VF3BW/u1eklMgOvyV38OmODv9CGpdNnnvzlnV/eGf942CiGWsLhD2B2SRN4t6RVAL2MpFTIjs9jEuDdCLF/rCep/Q+sWvzvUPcG2HxD/G0LxNtOpDPvjNzRyyECn6XCnEDbzNPESkJYhY40jS9zlzL/JL/f2yIp9uGlC+5V0BNasR9PPM/0r0EGcJ1pHIbUNrFjenINKrsficQpBWbMRcx22VB3MD5zBel3QUZwBDmMzV8JlXihnsikrF3krYXpBjY/p3J3CtL4gbg2nh/Ei7YJmqeL3VLtLP0R0GAgH3Jj+NW6sTdlLmYuVPMwj7FzuXWFOx1jkrrFpAg8lH7G+aRcYgbaXAQRz1fT5GfoupT40xVWqjEL0oQN3A989EYhd26krBT7HIpiBynMd/3+dnAxf0acxyp+NKOzNdJVRz0wjxS2UxbF3DrEsWNfb+F+vxxiLufh6fSv6RvCkoACnEgbfFBcYh7Bz3yFoqqoz7mVUztJ+a+AU7LKzCMVMRYHOLe30Pcd0v/FJQADIITi/h8oOJGCZOlLsdxrL/L8TItfBQJ/5SUBRA5l1sFIO4t9KygWwzirsNinugQ9yo+3kH6p6CWgFZmUXHVdAIV90Dmty7HcWyQy/GWem1+H6lKCbjw90jlgqrr+OzNWuBbFbC+jzfdZIrBkMVCLuOR+vq8yC3rCbGYC0oEdoEeL/Jvai3uk6myEQydd4CHuKd7iNsNbfT3OUMCX2auYS4pYL8avgui/pB5qvQN3y+lRvwyuoNfSouYP+ZTRtf8Eg37mpoy/5K/LqKjM9gs7XR/3yYqcWPa27mA6AB9SaU9cmKFDyE6P39yidPyCuylp+ebJ1xU6MAnME+03Zcl0VyHZT3oYiu4Lcbn0qZKAgTLui9j4kYVkaG1+LtAp+Xb6lG30EOsC7Nc/84JTMUfcEzBK76jZ0DiBuCmd0OChd2FO+38gg68NOQwQa/rqKejmZzGwCVxxV+zkO9xuZ5yu3qf6jQBlvHZzE3jFjfwAqmtMAgc22ETmM8W/P48TaC5nm7Ao2w7PQMYr6fSuQDFjfIq2HPvkUhxq0opzg48NSZxu9UnXxqbuFX4pFsml8MzIGwEh3zP/F0t/z5wccPQNVpPrZdrI8AWjinGcwX/312/ACo+P4ZK3+d2w1mkDHX1Eyju1R4dOPKlBI+UT7pMyx+McZlwk0cV19YZEPfdWg+UFHGHhVLFjZHnTUpgLnPurJ+5dOD5fkZLW73dgxy9m7CgRuvp+Zp8gsMY7RV87i1Qb73guWwIN6NoYoCdox9JuVOLuH2gk14GJGq/GVPM/GhUeU15ag1/c0BB+WJMm88JfIrOTMjzqZdPPog0yUS7ZkDYDfQss9Q2zZS4AcR8w3e9YcIEvjs82JBMr6bMoDbqo0HQVY1M+5IgDbiLVHwGibiLn54jocMIU1seo7qHe+vfRRfG4wjmnBKn45kVNwArPXKd9zZU3F6BKTeJNowGDMnw2Dw4oO/LpLiB40htMzQ3UNwwMv3gUo+8i+jDWFh6RhlkopHMiptxzGP8PBcSlb/HkhnO3Nwgge/Cgh6nvbem2LXfCxUkA8hk+i4Fa8TMqrjtrsx1jjIs7zmSuidFyLuiMibzdv75UNFB6oB6dz+QcrYiEXfpkhnpUUStZ8KE3auK1xjRYNFDaoB0YNjP3iuE786suJ/1EPcxiRJ3Lve2i/FsVVL2oAUlAZ6b2M++OKTvz6y4z68q7JXlRAP2TJi453q4pm4n2jAacFR5h9SeNom4gxU3r63tBwrEvZjoKpQ+ReHyxgkS90su4l4gmVSNBnwtkE8PMRaWiDsUa3lePtsy92FWhNRB4EjyWDcha24UlpvlmJIfKfowGlfo6XjYuzNZF3cVYC2LAJNRlJDRMV/jGnnQic6z3ePiBeYA0YlwoGoTwblE3C5oRCquHP69daQ/CgLCuVrYUSWaEHF7oEwbPB4hWd8KSsdJpBKGRBnVJuKuBnDe/5x5i/RNQQk4ilSJq6h9KETcNQClkD5iIuIqJ/1UUCSO1yP2XjGcW8Ttc4oOIxu2L8R5ROAX8CREcsOOMZ1fxO0T2CqDke1JEbjAB5AUBP7icdaKF3EXAVjRX9Oj+JbSfwUuwNINabRR+menmK9FxF0ksDWGrJSorCIlggSFQLLKV0iFbjZJwPWIuGsJeBnN0w9QIIDhFXWzX6DaFRAQcSfsmk5gLqKqpY8EASNfZimXe9q2rEf55z0Sdnm/YX7HvJ2S5fQk4i4R7ZmoDvJPvSYXBC1sVS3VWZTgkIRcHtIPIx47iQUKRdwBrbX+zZxCEooZ9Ihdli9EUDUybmLMl9ZIv9DxYk9q3XURd0DAdGw4qX3Ng0SWgYl7R9d49lxuRoyXtQOpOvIoK904wY9PxB3CA52m3+plIs+SxW2xkGe7pHF+LKZLQimjhWRGzXcRdwhorKfpH5OqXCooTeC9HRVWpvKxljG06ROkqtWYkkJaxB0iziBlTcd0vYHItCSBN0OSCuaBdvSloP5AqljAnWSW0VTEHTJQReIp5nQytMJJhoHEGK+SChz6jYHXL+KOCFijLdBv/yaim0QDKbbO12vrm8ncLU4Rd4RA0vnHdacZQBE5POg0TWcwhzC7iXarRW+9rv6A2cPwexFxxwCkT55EyuDWK2Rht3YkWNyIXGyiYdcpONJqIZLrFEpH7L6IOyag88BtFdFD45h7hyJuy/qHyx7xGjv8zJumAE5HDzKXkarXlabty6yJ296COZT5EpPXU3bLmO+rnl6Pw+A2kQI2urGQP/IoarBHxkXdVo/US5ncHxLtjCLi9tHVGzC/cFQamc9smoD724T5F1JGN0QW7RPQyP0vF3HDN7tpRkUN7zJU+VjOvJfZOsX3milx9/OoD3ZJgu4T4YLwgpqh1+WoI15rw5suarDQMWpfnkFRwziGLDpLmLeSsfnfbZ7p2RfrmefoGgpXZkrcF3iI+84E3i+2YxBS+iGpxBAXMZvXUuAtmFfwKH6rna29dji79GdOIGUow0t8C7NvyX7c0Xc3VCPwTIm7O3Oji7j7JPzeMUXHFtoK5jPMI0iKJVSH3Uhly0E64be0wOubf1v50ldug9NzIm71gHgEtNfrh1LOvMOg1oXR58/ML/VIhHzq3UXLebTRNgtsLy7Uz6Z9um7R7uEh7kki7l8fEncE+2imyY2/h+7AM5n/I2Xx3SNjIzraD85AqAwDqzci8VAAIKV+/HYj5o8u4r4+yeLurN+4P0Qj7tiBqKL/BPRd2C/fUwv9G1KBKqNJOWI0i+h+Po5olIQb6GGkCkR8q0doCProAAV9rZ4BJFXgh6hy078I+01mWVLFjRpKc5hjMiRuuIBODem7sX87iFQWzlWkgh6w/jyRwssSg/xhYSTe30rbFxBVN17fz3+1ANFxw6jhNoJ5dcJHcPhq9GJ2hSdDNR+MXdynMVsxTxdxhzLS7Utq6+tFParjGcMo91fmMaT2fZMgbgTTHMi8kPmYfj6rSaUJvlGPzlHszRsgbt9IzJpbxB0+8JbfWY/iN5EqsDBPr1eRlneUHiXPJLVlhlKzWwYkbkwdsbfcU5//SuZDpFxv0e7YCUDqIjiWnKPtB3FYuEXcMYkb0TpjU0BMmWcn6HrgEYctow/0S2emXs+uZP7MLGeuJeXVtUQTGT/n6/tYq18Sc/WxxfozeGmsYW7QXK1/B3/6r5iTmW+TykWWlGcB+8GHKelnb0YlbgRJ2AXcuUhxo65xP2EsRCaSs/WUGdP5YZq3ufCGgt9fTCoC7WR5hrHx4CjEDc+gzgWsX6S4BQKBoRBxCwQpA0IuEZnzZ71ua63ZUB6NQGA25jvW4hU8Vh6NQCAQCAQCgUAgEAgEAoFAUBr8RI8lBYjIgpcQDIbIi/YUVR+lBc8wp1Hxi5Tdo7SL6MUVfqPHkgL4bcNvehdSftevM9+o5vPYBkRtsdYF3Dpl9yjtInpxhd/osSRgB/2GLwym6KiPtfX4m59IRXCZgtrco7SL6KVamHCx8Old6nIcx/q7HC/THWwk81NSgRujdeOk5R6lXUQvqbjYgaQyhTiBY4NcjrfUa8D7mF1JJURE+ONnpLKgpuEepV1ELyVHjyXhegd4dKLpRXT8Nvr79ktopwniHk285qS3S6LFbVr0mNv19iVVe8qJFXpq6Bcr9MsjiQjqHk285iS3i0zLQ8a2+u3eqeAYkiKWk3ulC0z5HnBM9Sq+o2dK7lHaRfTiCdOix5DdZJzuSNh2mcB8tuD359GvpXRRSQSZS+4nlcQQMwAkBERGkCSXkq3pHqVdRC++YFr0GDJ3jtZTOKQnepwql7J5TrMCu+uOVvF57E8mfT+1pnuUdhG9CAQCgUAgEAgEAoFAIBAIBOnB/wNopWkMeaxRlQAAAABJRU5ErkJggg=="},{"partUri":"/media/image2.png","contentType":"image/png","content":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAASAAAADwCAYAAACgyrkbAAAXbElEQVR42u2dB7hU1fHAx4ogGgREUBFs2KKoaBKNBRVFYyxR/McSaxSNBVsM9kBiAf1bMPYWbLHFqEETW2zYe+9SooioCCigaIy5457jnbdv777dt7t3b/n9vm++T9nde8+dOWfeuWfOmREBAAAAAAAAAACAvNAxkAVRAwDETfdAPgzk5UB6oA4AiJNdAvnWyVjUAQBxMsI4IJUtUAkAxMXNRQ7orUAWQS0AEAevOMczyTih4agFABrNwoF85ZzOIYHc6f57biAroB4AaCRrmFnPZoGsFMgX7v/vQj0A0EhsBMyH4EeYf9sJFQFAo/DO5mPzbx0CecP9+3uBdEZNANAIbnKO5qGif9/SzILORE0A0AhedU7mwhKf3eg++zqQ/qgKAOqJjYAdWuLznoHMdJ8/Esh8qAwA6sXq5jVr84jvHGG+sx8qA4B6YSNgS0V8Z4FAnnPfmR7IkqgNAOrB76V1BKwU6wfyjfvuZagNAOqBX2R+uILvXuK++99ANkR1AFAr/gzYRRV8t2sg09z3NW/QQqgPANqLOhAfATuswt/sI+Ga0VGoEADay2pSff4fDcPf737zWSDLokYAaA9DjAPqWcXv1jAzp5tRIwC0h5MlDK1Xy2jjvLZFlQBQLTc4BzK+Hb/tFMhE9/t3pFBRAwCgYl52DuTidv5+RzMLGok6AaBStP7XPOc8htVwnb+7a+i1VkWtAFAJq0p9KmAsF8hsd517USsAVMLOxgH1qvFax5tr7YpqAaAtaomAFaMpPV5z15saSBfUCwDlqCUCVopNpHBGTK85BvUCQDlecs7ikjpe8xp3zf8Esi4qBoBS2AjY4XW8ruYT+tRd92kp5BECAGiBjYANqvO1DzbXPghVA0AxOxknsXSdrz1/II+7a+tsqAfqBgDLSVK/CFgpBkhhHUjvcRXqBgDL9RJWuWgU50uYPXFzVA4Anhedc7i0gfdYPJAp7j5vSqHSKgDkHI2AfekcwxENvtfuEq41HYvqAWAV4xS2jOF+97l7zQlkedQPkG9+YRzQMjHcb2Uz47od9QPkmxOdM5gR4z1PMU5ve0wAkF/+4hzBozHeU7MlvuvuOzmQRTEDQD7xEbC4q5sONrOg0zEDQP6wEbAjm3D/WyTMnrg65gDIF/3MLGSrJty/dyCfS1gKej5MApAfbASsWQUFf2vasBcmAcgPJ0j8EbBSr4EvuHZonfklMAtAPrjODfzHmtyOn0qYPfEizAKQD/zM47IEtOVK15ZvAtkA0wBkG81O+IU0LwJWTLdAPnbtedG9mgFARllZwsXfwQlp0/6mTcMwEUB2sWWUl01ImzQM/6hr0yyJ52waADQBXzxwZsLatWYgX7m2XY+ZALLJtW6QP57Atp1tZmfbYCqA7PG8G+CXJ7BtiwXynmvf24EsgrkAsoONgB2V0DYOMbOgkzAZQHawEbCtE9zOO1wb9cBsP8wGkA12MA6od4LbuaKZqd2N2QCygY2AJf0E+snGWQ7BdADpx0fAnkhBWxcO5HXX3g8C+QHmA0g3z7kBfUVK2jvIzILOwnwA6cVGwI5OUbtvcG3+OpC1MSNAOllJ0rnJr6cU8hZpu58MZH5MCZA+tjcOaLmUtX2YafuvMSVA+jjODeDPJH05mPX10a9fTQ9kScwJkC6uMa8xaWR9KSQtS9MiOgA4nnWD98oUP8PF7hk0jetATAqQDnThdo4bvL9N8XNo4vpp7jleCWQhTAuQfFaUcBH3Zyl/lr3NsxyNaQGSz3Zm0PZJ+bPoAvr97ll0VtcX8wIkm2MlvRGwUmg5Z5898RbMC5BsrpZ0R8BKMcrM6n6OiQGSi4+A/TlDz9QpkAnuuSYFsihmBkgeNgJ2TMaezeY3+iOmBkgeK0h2ImCluN0927xAVsPcAMnCRsD6ZvD59FzbbPd8D0o2FtkBMsNwyVYErBTHGie7GyYHSA5XuYH5VIafUbMnvuqe88NAumB2gGTwlBuYV2X8Obc0s6CzMXtT6RzI4EAOF9blco2+cn3uBuXwHDzv9RJmT+yP+WOjk/sDcGogjzr9+z8GN6Ke/NLXdITtcvC8vQKZ5Z5XBwLZExtDx0C2kMLWh/FSiEB+GyGkTskx25qOsEJOntlmT9yfLlAXtET2wEBGBvKwFApGlnI26vy1qOQD5t/WR3355XeuE8zO0WxAn/NxMyB60Q2qRksg7RLImECeCeQ/EQ7nk0BuCmRo0R+4FyU95Z+ggfgI2NM5e+51zaC5hm7QJou7V/RR7XQ4lg3M9/dBtfnmaclHBKwU55mBsDldoazD+TrC4Uw3DmeNCmfRY81vO6Lq/GLPgA3P6SCb4p7/zUA65LgvqCMYVIHD+bTI4SxQ5X26BTLXXeschmC+WV5IV7Gb0cFxOXruRSp0ODOcw9H9OgPa4XCKOUrCnN2rMATzTR4jYKW41+lgbob10ME4nEfMLKRYZtbZ4Vh0z9nb7j73M/wgjxGwUqwsYUnqf2bU4cyR6LD4OPcKXm+HU8wW5r67MPxgrOsMz6AK+YMZHDuksP0LO4czws3oohzOXPe5dzgLxtjGm10bpsR8X0go/gzY1ajiu0XYd5w+/i2Fs0ppcjizIxzOF010OJalJVxnOoXuBnk7A1YJg83AHZ2wtml9s42crdLgcIo50bVP9w/1pqtBX8nXGbBK+auEh1XXamI75ncOxDuczyMczpfG4WwkydxXs4CbVWp7x9HFQPmZEAErhR7LmOn0Ml7iS9BW7HA+q9DhdEqBTm3GzW3pYqAcI0TAojjaDJi9m+xw5rnPR0hhzadTCvV5p3uWCfQ18PzZdYpnUUUrdN3keQnPNnWvwzXnMw5HX0M+inA4WkxRw+ajnMNJeykh3ez6jXu24+la4HlSiICV40dm4Fwco8PpnDE9nmZmckvRrcAPDr+oeSzqiORypyN1RBtU6HB0B7HuJJ4W4XD0Ws9k2OFYOhjHewPdCTx9zIDYHnVE0tUMoJekdUhbD2MOdQ7nwwodzmI50t+uRg8D6U7g2dp0jJVRR1kOMrrS2eLugVwayFsSnWJUNzReFsgeUtiAl1fukzDTALXY4Hv8iWTdtLYA6ohE1yz0tPzcMs7GR3euDGRPYZOdp58UTrx/6/obwPdc4TrG86iiBT2kcEjyAglriJWSSVI4R7e3e52F1pwl4dmzbqgDLD4f8nU514OG13cO5E+BvGL+YhfLe1LYruDz2GxGFyqL7saeLvnNtAllsBGw43L27FojXheNdevBu2VmODr7ucTNhnq63y7mHJF+rjltFqErRbKfUPECygzCNKedqIbeNTicUuxsfvN7ulIkPssCaV6gFduYQbRSxp5t2UD2cs6knMOZYBxOtSV5xkl4JouUoq1ZT6i5BmXw55x0cTDt53K6OyfSlsOZaBxOrWFxXXT2qTDupju14jIJk9d3Qh1QTJojYN0qdDiT6uhwSnGSkFq0FEtImInxPNQBpXjCdZBrU9DWrhU6nMnuO/r6FUdqEc1G+Lq791QpVAkFkcMkjBSuhjqgGI2AzZLknkxewjgcXRj+JsLhaHIrXVgeKs3LZbSphGH7s+la3/Wt15w+HkIdUIreZhDvmECHE1Xu9/0EOJxS/EXCNKPr5LxvbWLstRtDDUrR7DNgXZzDGSPl64t/LG3XF08CelRjhmuzhp7znGzLO+Np7hUVoBU2AhbHGTBdG7H1xaMczicpcTjl1j1UDshpv9IjLPOcDkYxzCAKHwF7oUHXX7zI4USV+52eYodTjM56njTPtWQO+9VwCVOPLM8wgyjqfQZMz/xUUl/8U+Nw1pDsncBfz8zursxZn1IHPEGyVVUWGoCNgJ3QzmssUqHDmSGNqy+eVC6UMAQ9MEf9yq4rktwOImlPBKxShzMzhw6nGF3v+sDpQ0/WL5ST575Nws2f5JaCSGzVz34R3+lgHI4mSp9bxuGMk7D6Jh2vwJ5GR8fk4Hn7mFfPkzA/lMNnQbQRsOL64nMiHM6sIoezIOqM5F9OZ6rLvhl/1pESVvboiemhHJdKuIt4WCB/kzBpVLHoq5YuWJ/uZk6Lor6KWdO8qt6Y4efUP2I+P9KtmB1KoQvPGnU6VArnpaLOUuk0WkPJo6WQrqMzqquJ/ze63Sqjz7iDecbBmBwUXfjcSMJyv59LZfXFO6K6uqL69KHpyRmdQd4jVLwAKSTFOjCQ66VwMvvbNkRD8JzebjxZzp6om0f9QeHfYer8oJu+BpgZzmcRTmae+3yEFBaZdzSfkcUvPm4z9shSeoozJSzrRMWLDFNcX/zjCh1O8ZT/SPPqRQQrPnqb1+CHMvKqonvDPhEqXuTC4XwU4XA07Kn7dEZFOJxifJrMF1Fx7Aw3dtsjA89j9zptgHmz4XB0B7HuJJ4mldUXrzZK9Zi7zvWoPHYWdI5f9a/147uk/Hkedc/yLKZNv8P5sEKHs1iN9/RnwNit2hw00uizJ56f4udYx/TRgzBrOtB9OEOdw5kag8MpZhlzn50wR9MYa2z945Q+gz9wO0vYK5ZYVjAO54MIh/Nf43A0f04jc8hsZe67KuZpGhot8kEEtX3azs9prie/oH4h5kwOyyfM4RRzhBABSwoHmj5xSIrbvhambB59q3A4mvtYcyD3aGJ7/RmwlzBd09F9XI+Z15ilU9R2v5D+KGaMlz5SWX1x63CWSlD7fdTiBkyZCOxh1WtT0uafmH6+JyZsLL0rdDhaRsZX30xyKgJfteFkTJsYxph+tEUK2jvWtVX3pXXAfPVlWSlU1Gyr+ua7xuH0SsmzLW3avzOmTgwa5Xzf2eWthA9qXTz3ienOxHS1010qK/c7URpbXzwOthQiYEnll9LygHBS8YnsdPvAipitfR68EoczKQMOp5jDhQhYkvmHhFkqk1ieSDexvu3aeA/mqq/Dmey+s5ekvzZVFJe4Z32ZbpFIVpLCifKklrTZQtjE2iZdpWV98W8iHM6/c+BwinlEiIAlnZGmj/4iYW272bXrfWbQIUsUOZyocr+ar1YjWVmovtleiIAlH12AfsP8kUzKEQddhvDbBUbm2UBdnMPR0GW5+uK6zT0r5X7rQS+jmyGoI9HY4zJnJKRNJ0pYoGDZPBnjBxU6nE9wOGUZZHS1GupIPDeZAd+/yW1ZwM3GtD23ZV3xeshNz0eNasPhTMfhVIWPgGm2xIVQRypmrDOdzcZLc7MnbicZrnjRUSor9/upcTia0mJ++mhV+AjYK6giNRxp+v++TWzHnRJukkx9GtlK64sXOxzK/dbGeMl+cbysoX3+ObPE0L0JbdDsDpmpeKGOZ4ZE16Z6UAqJ1DcVzpjUG584fASqSBWbSJg98YIm3P+P5tV9qbQrc7S0LPf7sBRCepsJxfAaSU+j9/9DHaljrIQVa9eL8b66VuhTzGRi75h60DeNMgkHx8PmxgH9EHWkDk1YN93Z7+kYlyOGmH4zMCvK1GRcrxsnxF/kxjNMiIClnX2NMzgipns+IBlNXqc5diaaQbEt/auhXCxhziJIJ/MZh6AVcJdp8P3WMA7v0CwqdCXzfqmnfzejjzWMh52eb0IVqUadwlcx2fJcd5/ZUtgMnEl0PcJXBpgTyMb0sYbwqRABywpnmplJo94ctLKuj1ZfmnWF9jcDRHd+DqCP1RUiYNmik4TLF+9IYV9dvdnP9Jn+eVCqJrn2NYY+dlNNqA82AoZes8FOxqaNmNU+5a79ZN4Gik/GpCWKV6Gf1YXDnE517YAIWHYYJ+EG3nqOlQGSjOMfTWGwU6jPhdKXflYzFzl9voYqMoWWh5ot9U+PepmEh71zuTlYp5f+fJjmn+1FX6sJHwG7GVVkjhPMbOWXdbieJveb4653Tp4Vq2lS/QE4zV/cjb7WbvwZsJGoInMs7Ga2at+pUnu4/FAJK/nmvmrKwca7P++8M1THUkIELOvYw6rn1nCd+Ywzux+1FjjCDCCtod0ZlVQFZ8DywXUSHm1ap53X2Jg/VqWxVQLuk8bse8gqfkr9lZuuQ3Znun7joIbQ25OszzuxKUK0tBU2jcddQq6gSrlQiIDlhUPMGBla5W/1gPg899tTUGXp99OLjIJvEeoSVcJDTl9/RRWZR2c9T0iYTbRHFb8dbl7hlkOV0U7ocuOErhLyQrcFEbB8MUDCIg5jq3BcvmLwOFRYHk3EdINxQhdIBpJkN4geUt89IpAOLpAwlF5JhomtpfGHWzOFLpDdYZR2DiopyWZGR2uijtygZa2mSFgBpa0F5dvcdyfyRlE5ukX8fjPAKDXcGr8oqbvKiYDli1+ZsVGukkUf88p2PGqrDk1L8HCFis7zVPx1VJFL/iVhnq3lI77jt7hkouJFM9Ct58+Yd97foJLveVDCiCHkj34SHuwuVU7ZVrygVlwNaLG2V40T+jUq+Q6fafIPqCK3nGreELYr+szmFBqIqmpDp49vSLiXIe9RHxsB25XukVt0rdSH2CdLIdWq5x6hVHddsZU29OjBz3Osi4HGAa1F18g1Nsx+mnk98wdYD0dF9cNW2tCFtcE51YONgHFsBW41Y2J1CRPb6wJ1F9RTX9aVQnJ7n+R+vRzq4E/u+d+kO4AUomBzXZ/Q6Ng0CU8TQANYW1pW2sibE/IF7P5GVwDHMeZVzMuPUEvjsEnup7r33rzwoXvuU+kG4NDNqK8b5/MCKmk820qYZiAvSe67m062O10Aiv4o+8XnA1FHPOwi4XZzLeS2dMafdxPjgNbG/FCE7v85RTieEys7S8tKG1l2Qj6XtjpdskcCJIS9Jay0oZsWs3r25XwhAgaQSGzaSl2Ey2KlDZ8lgAgYQAI50jihxyV7lTb8Hg9y+wIklFOME9KNWVkpQUsEDCAlnGEG692SjSMLNgLWHxMDJBfNJX2xGbC6ZpL2Shu/ESJgAKlBc+Bea5zQ1ZLuvLj+DNhbmBYgHWiljRuNE7pC0ltpw0fAbsWsAOlBd4XaShvnpvQ5OAMGkFI0EvaAcUIjUtZ+GwHbA3MCpA9NWTneDOThKWo7Z8AAMoBW2nhWwiT3B6ek3QcJETCATLCktKy0sX8K2nyehIdtASDlLCNhRQGdVSS9uoQvRHcbpgPIBssFMknCShvbJbitU6Vl5QMAyAArm8Gt2RW3TmAbu0m4AP0rTAaQLdYM5BMJS5lskrD2bWwc0DqYCyB7aAWBz9wgnyXJqrRxoGuXJlzriKkAssmGgXzuBrvWXv9hQtrlI2DvYCKAbDNIwnI/mvxr1QS06T7XntsxD0D22UEKUTEd9O9JofJkM/GL5KdjGoB8MESSUe6nq4QL0HtiFoD8sI+ElTa0CkUzKm1sZBzQupgEIF8cKi0rbXSN+f42AtYJcwDkjxOME3oikMVivPcYIQIGkHtOM07oESmk9oiDe909/44JAPLNmcYJ3SPxVNr4QIiAAYAUcklfYpyQ5mZuZKUNImAA0AKtqnGdcQzXSOMqbfzU3GcAqgcARStt3GScw5XSmEobQ4UIGACUQCtt3Gmc0JgG3ONcd+13UTcAFKMn0x80Tmhkna9/j7vuOFQNAKVYPJCnjBM6to7XnuKuOQo1A0AUXSSstKFyVB2uaSNge6FiACiHVtp4TcJKGwfUeD0bAVsP9QJAW2iljQkSVtrYrYZrHSBEwACgSrTSxmQJK21s387rnOOuMQGVAkA19JOWlTa2acc1fATsDtQJANWyViDTJay0sWmVv3/f/XY0qgSA9vBjaVlpY/0Kf9dFwgXovVEjALQXjWbNds5khlRW12tD44DWR4UAUAtbBvKlcygfBbJaG9/fX8JwfmfUBwC1smMgX0tllTbOdt+biNoAoF5oXXef5L5cpY273XfuRGUAUE/2da9WvtJGzxLfec99fgbqAoB6M0zCReYXpWWlDRsB2wdVAUAjOElKV9rYQIiAAUAMnC6tK20QAQOA2DhLWlbauMD99yRUAwCNpjjJvQ/V/xPVAEAcaGmfW40T+tbNjAAAYkGLHN5lHNB+qAQA4kQTj10VyAPSMjQPAAAAAAAAAAAAAAAA9eR/oeqV3auxALsAAAAASUVORK5CYII="},{"partUri":"/media/image3.png","contentType":"image/png","content":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAUAAAACWCAYAAACvgFEsAAA0/klEQVR42u1dB3wUxfffcne5Sy69kt5IJUB6CIQmvRNqaCJVQAQBRUEEAcVKVaQJigXsBQsIP3sXAbErCgoooNj/CneXm/+8u1nYLNfr5m6+n8/7QG53p7y38503u2/fMAwFBQUFBQUFBQUFBQUFBQUFBQUFRfPDbCx/YEFY4vxQ/14sB7Cw1BSy1Ucxlo+wvINlP5bO1EwUzQVfYnmdSCMhuvfJ32ewRGE55EcCPIzlOBbewfNPYTHIUM+utMvSNc7qwxd4ltwfuVjKsbQOEP1TBAH2if7/D7mRM8nfO2VAgCos6gC40T1FgM7qwxc46Mf7gxIghVvgbBCgcEwgwFuxHMHyE5be5Ng2LHpyfD45doocG06WRA8Qz2WkqC5bxwRMw/K7aHAtxfIf+XsxuV7cloXkuJGQ9zg7dVlr+yPEGz6HZQOWk0QHNS722Vq7orG8gWULWdqC511roy9SfQD6k/qgzR9iWYZF6YCuYrHswvI2lk2kHQOs3CPW6gBcLWrT06T/YtjqoxT29G5Nv67cF460ay5pw3tYNpPyz2L514r97ZW5RnQdjKUfiDffh/T9H2KPWEpL/oGUABkJAQ7Fkkz+/7GF4wvJwICbpZTczA+Sc54kxi+yc4yxUnec5O9e5EaRtkU609ury1LbAcJzzzZYssjg+VE08J3tsyUPJAJLPfk/T875zo7XItZHCSl/Nzk2hRy7yQFdXUf+HkT+7oClwYL+7dXB2Fkh2OujFLb0bku/zt4X9tpVS8p4i/y9iPy9xYb9HemrcN1gLGXk/zDukkQ2mUOpSJ4E2IL8bSDPB6XHs0W/XU9++5rMum9i+YzM4LaOOUqASVbaIr3R7dVlqe3CQNSL/v6VnFfkYp8tkRlPBhB4Cw9h+ZNcr3aQAOeT/68nx3qSvw84oKvR5JiReEzQdo0F/durwx4B2usj44TebenX2fvCXrsWSPrdzwoBZjvZV+G6BCLw/y/IsVHk75WUiuRJgHGim+tXOwPgBvLbZtFvcDzMzjFHCdBaW6Q3ur26rA3ePyTlWCNAR/tsicwmk/NvJH9/Rf7WOkiAgsdwnx0CtKYrWPJuJ8s6OG+dBf3bq8MeAdrrI+OE3m3p19n7wl67Fkj63d8KAcY52VfxdXHk/5+RYyPJ36spFTV/Aqwk3oU4ZOM1LBl2jrlLgEdJ2SzxXuzVZYsAHVkCO9pnabvEy6opxHs4Ixkwlq6xtAR+RfTMFJFy7ekKwpu6kv+Xk/N22FgCW6vDHgHa66MzS2Bb+nX2vrDXrvbMxYgIljxntEeAjvSVEqBMX4Q8SW4aMMBLIhKcIbop7yIPhuFGOk+eA4mPP48lXlTu5Yw5zGYfeW4z2sFjjIW6tzvQFuGBNyw9XsByt526bLX9D9HDeHjI/QmWdg5cZ60uS+3KJDf/EbK8OyHxOqTXSPUheCbwrOsx8q/wgsKerhrI+dvJg3d4eF9g4yWIpTos6UI6edrroyUCtKZ3a/p15b5wpF1zSRs+xfKo6Lg1+9srU3zdHVhWkP//SfrxPvn7W4bGUlL4GcJApAhOvYNTMF709yxCTpdTE1EEOm7GomPMbxzXUHUEpd454hXvIHJY9BiCgoKCgoKCgoKCgoKCgoKCgoKCgoKCgoKCgoKCgoKCgoKCgoKCgoKCgoKCgoKCQp4oyIgPgWSN8VQVsoY6IVzRHoSRX7ZoiqaIJ2Mqn6pCvuA6Fkc+/dCcgr8OripHNw5P/7FdXsQMqhYZjqb4kJbDOsR/8vqtbfUgQ2rjDyUkhORQzcgPNQURMxeNSP8RxtT2awr+6lAc+RTTNCM7hRxQmh0+59N1FTrDMx2RIAuGpUM6IZqqW2boXxW75/xTdRfsdO7JOtS3PGY31YzsELtoRMYx8Zg6vKb8fNvssNlUNTJDv8qYHWJDgbx/Z6kxIULRk2pHXpgzIPVzqa1m4d+oZuSFpAi+94d3lRmltupdFrODakdmqM6PWH32kdomhtowPQ+SS2ZT7cgLozolvisdVA0dE96hmpEdcrfMzD8rttOvD9eiqrxwmv5ehogcWBX76S8PmUnwjVvbGuqKIx6gapEfWqdrh946Nus3/dMdEcjy0Zm/tcoIH0w1Iz/gMbT9rRVtG2FMncFja0BVHCR7jaCakSfGpsWFoPrcaBSnVUBK9FCqEnlCq+H3tW8d1VhbEtUYGsK/QjUiW4TFRqj0dTUlKBWPLfz3GKoSmYJjmJ1larX+cG4uYi9u40chP2h4nv33zptboTuWFCP4P2N5dzcK/wO2zESTV+5GiVlFeo7j6PM/2Q4qlv13aUIC+r2wELVWq3WceeNmCvlhEMsy6Mj+7iaB/+PfBlK1yNCp4LhHEzIK9XO2HUQdh89GLMfDZEXjNmWIAeD1fYK9PyDAxZgIFSwLO8aFUNXIbVAxD7dtFakznB6IQNoUR+jwbw9RzcgOakx4/1c3bBYCApx054swUYH0p6qR3/J3e6uQEB2QH8jHOTmCsfpS7cgKIQqe/eeWhUVIIMDlC4qQQsH+jY+pqHpkBRg7aOIdL5gIECQutaUez2APUtXICyrs7f29KD4eCQQIUoAJERPjNqoeWQE+qUJfvtftAgHC/8lk1ZuqR1ZuxbbYlFy9QH4g7etnII7j6WQlM/SCAbQfe31iAlyACRETI2xZqKAqks3yd2tRfviF5a8ghXnhsLPaFqoh+TgVQHS1g6cjMQFeseJZYbKiHxjICJvzVCqdmPxA3s/OFozVnapIFlDgpe4fi68rQFICvOnaAlgG08lKPoAxg8bf8nQTAgSJTsqEyWoTVZFcBhX28ubHxSEpAYLkYGLE52ygapIFusGgOvxm10sI8JM3ugiT1WVUTbLAhqjEdJ2U/EBqBk5BHM//TicreaArDJx3sbdniQDnYWLEBHkWn8NTVfkd9+Vmh12y/BUkNysMJqv1VE1+B8/xirM1AyYjSwQ4dunjwmTVharK/7g328LyV5A3s7IEY3WiqvLvoFIq2LMLZuchawR4w6w8WAbTycr/gLGCie4xiwQIEhmfBpPVPVRVfh5U2Lv7dY6V5a8g6UolGGstVZdf0REG1f59na0S4Ed7OwuTVR1Vl1+xNiI2WWeN/ECq+l6Bl8EKSDZCcwP6ER1gwLyOvTxbBDgrNhYpWfY0PpelKvMb1mSkaqwufwVJx+fgc1dTdfkNHCa205V9xiNbBDh68SPCZNWeqsx/WJWKvTtb5AeyLzNTMFY7qjK/gMXL39PzZuQiewQ4d3ouUiq4U3Sy8htgjKBRix6ySYAg4TGJMFnR1Fj+GlR4+XtqJvbu7BHgb1iSFAow1l1UbX5BDQyqd1/uaJcA33mpozBZVVO1+QV3aaPidXO2HrBLgOU9x8Iy+CSdrPyDKhgoe7F3Z48AQabFxCAVy56gxvIL7kxOVOv0pwbaJUA4Jwmfi6+5g6rN904FJrQTZT1G2yU/kJELHxAmq0qqOt/j9gTs1f3mAPmBvJyRIRirnKrOt4NKpeJOXD052y75CTITn4uvOU4nK58DxgYasWCrQwR4zdaPUWhEDExWt1HV+RjYmzs+BXt1vztIgGexxPI8GOtWqj2fogwG1evPdXCYAF97toMwWZVS9fkUt2rCo3VAbI4QIEjby0biZbDyGFWdb9EWBsgL2KtzlABBJkZHwzL4KFWfT3FLXIxKp/t5gMMECOfCNfja5VR9vgMmsqNtug53mPxAhs3fLExWrakGfYdl4M2ddYL8QJ5LTxeM1Yqq0Eeeuor7/srxWQ6TnyBTL8+EZfD3VIM+QwmMjaHXbXSKAK+5/2OkDouAyepmqkLfLX+/uwJ7c787SYC/Yok0L4OXUC36BDDRoFeerHWaAPc8UStMVsVUjT7BEiAyIDRnCBCkdad6xCmU31AV+gZFMDCext6cswQIMjYqCpbBX1E1+gSLoyKVuvM/DXCaAM+dHIDgWlzGTVSNPlj+KpRflXQc7DT5gQyZu16YrAqpJr2PReDF/VJQ4BIBPpmWJhirgKrS68vfL68YleE0+QkyviEdlsFfUE16HTAWUP2ce10iwNlbPkIqTRjswngjVaX3l7+fj8ZenCvkB3IGE6eW48BYC6g2vT+oXni0xmUC3PVIjTBZ5VN1ehULgMCAyFwhQJDi9v2N2Iv8jKrSu4CBgB7DXpyrBAgyIjLSiIn0MFWnV3FDuFah/+9Ef5cJEK4NN+/vfD1Vp1eXv4eLavsZXSU/kEGz1wiTVQ7VqPcwH7y3My4ufwV5NDVVMFY2VanXlr+fjB6aanSV/AQZNSTViMs6RDXqNcAYQAOvXo3cIcBZmz9ASpUaJqtrqUq9t/w9OAx7b+6QH8gpTKCh5mXwPKpV7w2qpx+sRu4S4FMPVNHJyruYB8QFBOYOAYLkV/eCZfDHgaScSBm1JQsGwkPYe3OXAEEGR0QYlSy7X0b9U2IJdfFaSCCqlVFf5oZqeP0/P/Z3mwChDA0uC5c5R0b90zKuJ20NJbaWA0KBsPKrehrdJT+Q/jPugonKyMgrdMl5DivNDGuY1D3p8IZpeT/N6JP8RVWOdqYMOnKNmmUNP7u5/BXkgZQUwVhp/nZsOxRFPrFiXNbRtZNyf+xdFv1GrIZJcfBatiY/fPOiERnf3Tu15YnBNfHv5SRq/R7krVLyHw0bkOL28leQoQOSG5VK9gN/9ysvIaxNPdYx6Bp03i4vciPj4PfKyTGatD4VMW+AjW8Zm3UU2/wxxk9bS7aIUmf0LY95c/WknJ/mDU5DXTq3Q+68ABFkxLQFqKFjonHN5Nw/J/dIPlyaGzbWX7aqzo2YfVXf5C+AwyZ0SzpclqMd4dCFGg2TetPIjGOGZzoiQe6ZmnsqJy7Ur0kEsLf2/sDw8EZPkB/IT/n5KAQTKi76an/2q7YgYs0PW6obBV3/+0Qd6lMWs9eRa6tyw+cfXF12TrhW/3RHNLhd3IeMf5MIwIRi3Lm5EnmKAHdsqkAsa5qsUvzYL3ZIbdx+0LGg7/0ry85V5oQ79BilT3nM//57su7Ctcc2Vze2K4xY5Y+O9K2Ife2cqC1HNlWj+nGT3SK/Cbc8jlZOr0Bi3lg2KuvHKDWT6fPnL7GhlfdNa3la3JaFwzPgE9gWdi9uk6WdeXRzdZOO/N/jHVDLJPXfmDBO+EvwiDaGc5whWaE45ylRsKyRY5h//NmvQVWx58W6BoEJCOYie7YaURf/ovTa9VfmNoaouJ/VIdwJf4hSyf0CjyrSkjUoOyPMI5KKy4IyoWx/9Qt0umF6y0apvoe1j3/RgTGpvrkh85j02v6Vced5heqEj+Xk/CHpBmlbJo3s5hYB1o+biv7a2aFJmSe31aDkWPWfvu5jbovQv8GRELflu01VqFWWAyvZnETNyIOry43ii396oB0anhqNFick+EW6hoUhTFTwVQCEQ8z3oGyFgQVJVf3Vt1FlcUh6M84ekPI148BeCwMqY5+WXrtkdCZaNDcf3XpjkV+kRaIalbeN8ni5ZW2iTGX7q183Yp0uG5N1ia3wUvJpR57Rzh2U+s0l19Zmo8o+V/hUKnpfjmYOzrukH+OG93SLAIdOvhYd31rTpMxP11agNjWdfd7HsopydHp7uyZtwd56Y1a8ZpgjHmTIgKrY/X/saH9hSda/IqYxQ6nUvWFn7w1vSW1oqJ5nmOe94C2HY5Y5vzIpyS/9WpGYiOK0ysaN01temJEPrio7364gwqHdtgqS1J1WXJ51Srj22w2VxqK0UFTfr0Xjr9/08dgS1FE5frgn4jgWbV9f7vGyH7y3DPE8i376vJfP+/XL131MOi1OD0Xfb66+4BxgQvw5L1nTwaFHHfkR6w+tLj9/0VNvaYgODzF2HnUt8sQLCGdk4IhR6NDai6u8jbPL0Oh5d7tV5lXr30Zj+leYVotQJniDw/tWo6s3vu/Tvo1ZsgNFxCXr+lfFGgUvELisf2XsR848c43rURrzwNguif/rVRazI1TJlClZ9n94ydgIXstZH5LENy1bgvcHz3+88kAVE+uLHUNDDb4kvu/y8lBvrdaAl/WNuAlLitPCGkbWxb8wtnPC7qp87UJnnuHlpagvG1ob99zYzol7awsjb8c/9YKlYnKSWv/mrjqfEsX6O9sgvFxEv3/X1+Nln/22D8LLUHTfXW182ifQYUqSWkeW9r06YB2DrkHnOUlqZ/bF5WryI24EGw/vEP9CcWrYKPzbGnhemlPayTD9ntd9RhKQ72/YpFlowqg+aPyo/mjMtas9Uu7UlbvR8JFD0YSGnmh4wwg0be2rPu1Tu4FTEcvxjSzPf6BRMB16l8XsBA7r3iZqG9ZzrNsPgbHMwmRk6BQWZvgKE5MvyGJ1ixZAgHovhuSMB4I9gknJF/15MSMDQSZr2M8E193ZS31KUCnYPdgba7xpXj5yJhefO9KtUzzq0z3Ra+X37paIundO8ElfIIED6A50qODZ3Vin8V6y1QA8aH8Li4zVDb/+fp97g4EgU1ftRWmFVQZMUeAoLcGi8OZLlioVy/4QzfP6J9z8JM0R6YLJFntpe7zYnygg2HWYaL3ZD/h6ZWpMDIKXObg/T0K93n5raZqwOFZf1y5W/8Ohnl4ljNNf9YaXFGjL6lKv1bF5VampjjNf9/ZqX0BXdTWxetAd6NAHb9UTWZbfw7KsEZ5heSIkJVhk8DX3QD5CbCsetr3t6qs3zRF4EO+AFwhTYmKMpz0UmyeV77FXxuNlN65nojc7g+vY1w0TrbfI72BuLipVq/W4nnO4unE+jgooVyrZ7yMjlLpnH6r2GmkA8QE5ARF6k2QVChbdv8Z7JAs6Al2Bzhjf7h9jmrBYltMnZhXpJ96xixKczU/vPkSl3UaCx4fw5PGiFz10mxiHB/W/JWq1bn9OjseJY31yMnhMBk+s3+1gChDtUS8sg7ekpMBndwbsNUNapyI/xbCF8zyzGW6WSWMyjH8d6+dx4oClLyyBvb00vaxjPOrbPcnj5YJOQDegI6KrcD/ZqpLj+aPKEI2+75W3UbKzIFfc9hyKT8/XsRwHDsUUP8e9MvlKlj2sYVn9BkxYniSPnlqtQcGyr/mgDwnwQmKTB9t/Mj8fDY+MBO8VnmGuhVgwGXxNMwx7UH/n52r1h17v4vEXFPASxNsEeO8dbUx1/XbEcy9aQBdYJzrQDdbRUBnYKZzhuIfh3ils17dx5oZ3KfER6TV5OVKo1HjJqwCHQjaf3IVguQ2eb8HXGsc84En9iAlEaV7+TvNFBxQM83ZfD31p8mZWFspSqeBFx5+46EGMvJChUrIfqlScYeXyEo8QCIS9QIjKyc+8H6ICdUBdD93nmVAb0AHoQmX+1C5DZrYaxnH831EJabqxN+8MauKbce+bqGVFN5NDQd6eaxgZojsmrV9TFArdKw5uWG5r6UjCRBJ91PaZeJlqOIGJ19U2/0Zi+yBcCBPqO4z/vzO2wffMEvi8bGDvJIO7LxUG9WmBOtbG+Sw0pa5dLBrct4VbZUCfoe/kEzuvvzl0A5k8r3if4xUGU8zg1gNBR34NCx9E2ugEHZ4MwKEYzMgcCZgAXoHQkvlxcS7HDPbH3hgm03d92O4U8GC3YeJ1pb1ft2yJOoeFGTjzM8v5jOvZQnyJrnjZdzopUa1zZt9esfx5tB/SqHm0+pYSnxHgKuy1hWp4U92uXA99hT5D37EOujQDOynMJM0as9p0NExf93pQEN/F2D7OyPKKt2XsUFiNGdR3CA3Vf+lkzODJi4kKZvmy0ZhwPxwUEWF0ZZ8RCAvCHuSPuJhqpnkhTsGzLwkxg85uYARJD7AXhbwdZiMNU4E6H9tS6WpsnxH3Gd4cxjUzW3XleP60Whupr597b0CT35RVr6DUggo9a+aB5uJQXIIKTCrHYBMjZ9LYP2jO2AxLk1Qft3euMym3Tl+M7UPYOjsZeeVPdHbCmsLz7Pn21TH6Ywd7OEwqwwemoJqKGJ9/nlZdHo1GDEpx+PxvP+qOKkqj9LiPsnhz6M6ExfL8C/AsrLRbQ0DGDA6avRaFhIbrMNn/0AwdikvfaGFP0PRGC2IGTzlALvX+S1aaCcT7iANJVz/OyUGt1WodJr7/GN/H9nkLxbCTW1ioQv/opgq7pPL3D/2QNkyB7lhS7HMCvH1xsaluR5KuQl+gT9A3JjD2GTZNWCzHn0vIKNRPuH1XgMT2fXAhto/jOIgzjmACCKaYwTyVSv9udrbNdPUajgO31y97C2Di/WS4nbT7G5OToY16fO6nTOBtr6nheWYd3ISwp4etmEFIeQ9L0SP7u/ucAKFOqPuZ7dU2n09CH8wDyhSKpAkwWxVDOAiktO895ZbmHdu34lkUl9ZSh0n93wByKC71sDBpfARvW+FtqSVy2WHetxdu2iw/tfGGMExulr5uOXExtg/aB6/iQ5jARb1Cwf7VMjtMd+DVzhYJZsywNFOaKl+TnyClrSPR2OFpFo9Bm3Nx23Ef4M1hfQDbSU3uRRIz+E6zIz94u80rlAZM5rBTY8BvgWoOwcAkAm96pV9fjMTeF/Gs/IWWcDNJv3Pem5mJUpVKHYT54OM9mOBAulLBvqdUmmMG9acuEsy/J/qjqEglWr6gyG8EuGxBoakN/4q23oQ2QluhzbjtEEWQHiS2GgxhIpHxqTpICdUciG/6PW+g3LKujQxrSmJwW4A7FJfgMkwmZ5IVCt1uEjMIXleYebe2hf5sGG7Xl2PI5usQxgPhPBDWg3/fx/guLlFeExbLNPbrkWQQvvWFDc9hovjyvW5+I8Av3r3M1IYXd7S78K0wtBHaysg7ts9bSINwEUgJBeEjU1a+jOrHTEBDRoxAY69f59+3une/hOpHj0f1I0aicQvuQyMXbEPaqHgdxysgzVg3JkgRzzPMy5hcGoFkHjcvf0Hy/NyuxRGYiD/NzUWQK5Ck4/JFVhA5ozP2qE4lxofoXn2mPZowKgO1KozwG/kJUlwQgSaOzkAvP9YOxUar9NBG3NZOQWwnyCQ+S8FzhlFdU9E/JEHp4zdWomETZvqF/CaueBLNGFqKhASlD8wpQSmxaiPL85BmLIEJcoDBFkLSA0w64GV9K4M2tQYi1nIcfIv8c5APKDESeZ59BWIGw0J5dOOcfL8T4MJr8hC0BdqE27YnCD10i+jYKuoN8cZHINc0VJiyqPiaAIcP7Yd0TzdNwT+hWxJ4fkpqqYtoj8WICec3/G93Py/5lgEBYlKGZ5Gx1DRNAF7wStAPfP724yc9/UZ+EBANbSCrhruD3ENvgtGdE/ZJ9/3YMLscTb7rZZ8T4KSGbpfsQbJ6Ug44FvHUUqJJC27kutBQc9gCw9zJ+H7f1AyV+S11Y6VGA8kYzjAObE4UhFjXIlGtz8vRorgYFfJmnkFrAuEvUDe0ISlBDY8o1lKzXESH4sj10o2Bxg+u9su3xENGNqCzj9Q2acvIuoQD1EpNcQ9kUBE2Ldea4+y+JstRX2AK9j7/bRUSov8oJwfBBlDEs6ijpmkCXqFgz94wK88UhDxzcrYpHg9CYlz9NtcZ+eP7vqa6oE6oG9pw/ayWSKkwvZ2nk9VFhHVqFXVo142t0MHVFejqhkp0+Q3++XwONkwa1KcOvbK8DH28qhyN65L0W0l2WDdqItFzQEw+v86Ojb0QdnI4NxdVaDQ62LWN8e43wVE8wzzFkuQNZ0Txf2lKJWzFuZqapwlgQkAf7b0YG/jUA1XwAgIVtNSi/fs6e438oOz8XK2pLgjCFn7/8JVOwmTVgZqnCdaGR8fpR12/EV1139t+T2ww7uYdqKS2F+J4xRn6uOLS53/oVcm2m79gMhJCUDBBPofPifF0vbjcE7E8r3smPf2SwOerMSHj46epsZpgdXqqRiclJ/h+GJ7HQbJS2EfXkxsxQVlQJpQNdVj6Vjk9RQOT1SpqnovPajHRnIZ9gGWV2mrRdmGyqqEmuoiVEGRs7dOz5zMyULxCAdlWfibPCt32OLEsgfCby8LC9N9YyVYDwc/EWNXUROZBpVRwp+ZOz7VKVCQI2ZQe/8Sn7idIhTKgLChTGowtljnTchEJgaGTlRlAMCbCkVUA9NYDKCwqDiarO6mJyKDCxHZyRkyMzeQDsG1lT622UdhLl3E9PU6KkmXfhKSl8Dneb3YSnCYqFGCsO6iZTKiCQfXOSx1tktZ7uzui7IwwlBAXciFI2RWBYGsoA8qCMm2d+/aLdcJkVUnNZMKdQDRyTKBa1n0ULINP0MnKDLhh0R4HM0jfl5yM4Fti+KaYcf574f6Y+P5IVyr1r0mW29YE0l3h+o5TY5lwe1JCiM6aFyYW2LNj5ODUCy8rxJ+r2RM4V3i5AmU4sv8HtAnaxpg/qaLLX0wwQDRy/ARuxA33C5NVGTUVvmHjeV7nTObo97OzUZ5KpeNZ9h98/UgH6gjFy91NoPTLo6IajzuR9h42MifGKg12Q6lU3HEgJme8uG3rykyByuVto9BX79v/bA7OgXPhGrjWmbqumpQN+3scp0PKRCwmopFrpmdNeBRMVrcEvaWwd3VscnS00ynoIW0W5BdkzDGDD+F/w6xUUYK9xa/gG2NX0t0DMcNLElzO8iA3VVvQ9WvPOp86/7O3L0NtWkWicK0CPXivdVJ74J4y0zltSyJN1zhbD3yiRyarNkFuq1uAYIBo5JoIoU2XYXgZrPw+2PkPblS0C3tZrm5E9FBqKgrH5IaJ9AgZpGLM4hnmfLlGo4OwGlfruAITNC4/2I21LC5GpXP17e7/HW8aMwjxfMKxX7/pg4YOSIbP2dB1M1uaznX1bXFstAomq6XBbCggFiAYTxDV1JUvmz5lmz+xK5rY0AONvm6NR8odeu1GYbJqFcy2WhqFvatf3dyG8jNMbtUajR6THdz8EDMYo2CYZ4XYvl8cTHFvTSBEhhirOFgNhZeW30+9PNPtt7qQtCApIQRlpoWaXlyAwP/ht92P17pd/pRxmbAM/jaIxxQQiolg3CUpSLc/vF8t+j+SUAFkxZQyNGHpo+4vg+//2JT2Hrd1cTAvfy+knXJXIIB5VmwsfL9r5Fn2L0ys+p1O7ENir+xwc5quG4PUVJD1Gr20s51H4vq++7g7qq2KQeoQziTtq2PR9wd6eKRseOtMJqv8ILXVIpUmTO+JvULG3rQd7V5a2uQTNv3THdGQEcM94gUWtx+AOIXqs2DlP0h5dUniUXf25F2akIAw+UGojC7ZA3sTi6XB/4la/YmF4VqF/r8T/T1CUqe+7I369kg05e0Dgf8LeQbdFXiDDG3FbV4QlMtfhfLz4g4DjJ4gqFELt6J9y0svSWQwZMQwD21+tEaYrHKC0VY3gld1xs3lKcjneAlcY14Cw2dzsOMXfN72OOzWBmEsnqhDlKsw6DwLvKT8YnxDukcI6vmHayCLs06lZI/hoiuwVML/4Tc45ok6Lh+ZbsRlfh6snnr9nHs9sEHRh6ht1xGoR2k0Ov/UxbRad12RbUzJyDZMvOMFjyyxVepQmKxuCMbl72cNdjYfckTuT0kxZZGG5bSFZ3SwEdN/bdVq/YGcHLeXwVrzMvj6IDNVvulF1SPukRO82LhqUrawQRHsEhguqiOc/GY6x9WXIGKSJZNVbpDZagEQirvL3wm3P4/iMwp0LMf9p1Fx8/tXxe6bOzD1q4ndkz7ITwldhL3ML5UhGn3fabe5TYKFtX2NuLxPgo3/4MZE7jyj+8m8QZFpQGFvD3Yxs7bjVyFeun4Gu7htSk52iwShPky0h4LMVte7u/yF1PUlRRE6nmft7fg1Ds5pjc+Fa9wh27BQHiar+UG2/D1cVNvPreVv7ynLkUKl1uOybO14qGE4DtKPwUZMxpkb3nW5voFXrxImq+xgstV1oZiQTrm4NH0rKwvlqFQ6Bcv+xTi245ewc5YRdnU74UQgtFgeNm/YHlTGwkvJQ6OGpBpdJSP4fhcvoQ0qJQ9f7mQ6UGUmPnc/XAPXulpvQ30qLIMPBtGYgudomFBWu0REsJtcQbu+zu54WA8bMUUlpOvGLn3MpXqv3vQBEC5sgzsnaCyFPbIDQyIijK686IDvd/H1BgXDvMc4v+PXQPgULkup1L3h4KdwFvYs1geRseBTQ1O6K2cJ6Jev+6DBfVvASw4YUEsY5zYoEjZiMtb3a9EIZTlb/5PbqoTtVTOCxFbzYH9g2Fzc6QwtCx9E2ugE2KDIlR0P03le8S6+1gDbW7ry7XFeZXcjr1B+GCz8BzekcTv2ppwhn+/y8lAvrdYgSobg6o5fibC7G7wthhjBs06S4MDwcHgbHCzGmqNR84a/f3Au0embu+pQSpJap1SysO/DZW7U3w3KgLKgTGfaAG1Wh3DgWVwTDIbCBLI/v6qn0dlP0mD3OJbjjJjE9jKu76dimrAYPKay23YyTF/3ulME2G/a7cJklRYMtpodgj24n5xYhr6QkYESFApY8kK6o84eaAMkNpjFMYwBdn37yko6LEuyNSUlaIyFyeeDoQOSGx0lnfM/DUA3zcs3bVCk4NmXGc/s+RCPy9oNZULZUIej7RnSP7kR9+H9IBhTsBIy9ptxp+NfeKzai9IKKg0sy3lyx8NO2BM8FRoRqxt+/RbHl8Eb30O8QgWT1VXBsPx9b0B4eKOjb14hjIU1v+h4Al8e6eHmVKlY9odontc7Go8IxA0EHgTGSoYl6I5NFQ5vUFRXE6vHRKVjPL+FqHnC4lg91AF1OdKmR3HbyRI8OcBtdTWvVBmASBwhnMHXrDN9hcHxymOM59OHRbI8/wTLssbSbg3I0TfSOaVdGrEX+xYT8IMKz1T3O5CU4GBuLoLwFbxUPcfYfnPoLiIwue6AZ12QXOG0Ay9m+mACx0Qe6Ma6Cl5EiL/ZtSawKVJkhCm2D76X9maKo3KoA+pyZCMmaDv0AV83I7CXv4p3csu6NjoS21fabSQJReIg7CjCi80ax3L8uYSMIv3EO3Y59PYZz3PQtqRAttV0yOVn7y3sFkyQoRxnwOd+ga8p9FHbIGbw3xK1WrffTszgxuRk0yd3gWwsvHR8a0CvJJvL37+O9UOjh6aaQ5F4m6FInoSG1IUmjckwQhtstbF/z6RGpYJ9M4DHVAsgjj5Tb7VJMFeseBbFp+fpMCnZC0XyJIo4hfILZYjG0OfKFTbbN2P9W5AkFZ7vTwlYSylY9o1eWq3V5e9Jc2yf6VU8Z97mUO3jJuZjz+6whmX1G2zEDP6A2wkZpQPYWAnwiZqt1FWHXu8CW1LqFQpTKNJQP7RxKNSdn6vVQVtspdiCvkCfAtRWVwJx2Nr0CN7O4uWlAV6UMM4nEXYXeAxfiBlstBUzmFnS3sDyin1MwA4qhmm8zwqxQFhKljm27w987iA/thPin24DD28gXuoey8uz2N5uWq0BtzVQjTUFk0vj2W/7WI/tU3KNeDn6AePfMJMMaIOtmEFItQV9wedODsjlL694LbtNncXl7/R73kAty7s1MiwrhCIp/djUQeaYwTTd2Jt3WiTAHhMW48nK9Hw9JhBtNQm8JimhQBjK4oQEk0elYBh4rpYqk/Z2x97gLykKhc5Suv51LVpAItaANBYmjH29uiYYpGRy5uveCC+LDcSjggGlkENzGXPMYCO0DdoobXePLgkG3Ke9ATim4iH0pNekpZeQycgF2yC2T4+9QwhFkssevGmYsN8Bj9VSzCCEz0BIDj5vfCAuf1/pGhZmEJPI1y1bok74N0Ik8OZQbhtbJ0C7OZJX8KwkLpEzPwcMNGPFcBxr2LyqtAmJQCZo2HNDYd55rYsM290Ft+00tPH155pmrd60qi3iAtOzmMjxfOOMe9+8NLaP5Rox2eyR4dIfNjObD55eenG1YerqfU1IMK2wCi+D+ZcCjf+igeTWYK9JIJAn09IQhJ+oWPZHxrzbmFwhxAzq24eG6r8UxQxCDCG2ZqAZ63KeY41CeipRbJ+R59kX8PE4Gbc9TsGzL0JbxTGD0BfoE+O7h/++uTExwWW0qjEI5DFl5R6Ukl+m93Bsn7dQzfHKH9XaSH393IvZa7qNWwjkDaFUkYFkq3HgLcHWlqdFe3lg8tjJePdVvCdRgZfExyJ5XickcViZlAReYEAZCwikc/s40/L324+6o4rSKD0mPghFmsI0j53xoI1ToM3QdugD9KVTbawB/7YrkJwK8KK6j19kzq03azUKCdXqyD4b5c2kD5HYg30cuABiBiFM58o1+4AAgR9GBYKRwrMTNaOiwxSHajQaw8c5OQjCTDDx/ddMZ+NwsgGTicQ/zc1FPMsYU2NVD8dpFZ2aua0y85M180OUnH7dba3RIxsrUFioQq9ScZbSjDUHFEPboQ/Ql7UrWiPoG/SRcSwpg3zdXHyv4XvuEY7njJPvfhmVdKoXYvs2M03TjDUfB4nj/0vIKNRPuH0XSsop0UdrVR8Dd+Bj2uZppLiQvFGdEj/7cn0l+nZDFRpfl4hiQng99qK+YZr/rl1TMBGeq2oRZlg5IQed3FaDdi9u9V+n4qjnGPk9x7SL2oLwOXdekf3zCdyPfctaox5VsYixnLevuSEC9+ER6Av06dXlbdDxrTXo9suzf2qXHzG7GfaH69QqateeJSX/wT1354RcVJiTbGA5Djz0Sc18TJViEjwSERGmH3NZuokzvri3Eo3skPBZfHxI88vr2LcyZrc4kyzsJdC9bfRJxrE0O80BHW4Zk6UXpwo/tLr8fJss7cxm1o+oG4alHxX34+jmalSQGvpEoKwVC9NCn/phS02TtO7X1ad934wev5jQNjNs9uF1FTpxP5aOzoJHMO0CxFTqHqXRJ8X9O/dkHepTEdP8nrPP6JtySLqXwJyBqUeayXMku8hLUY/7dF3FJfsldG0T/VBz6odGwXTAnpFe2o+B1bHPBAoBDsJ9kfbvlZtb63Hfa5tTPy7D95a0H4fXVCCyVAwEsPMGpR2R9nFGn+Tml4S4vibuVWlHGjomBExmDpWKKdo2K/93cf+Oba5uLM0KbW47xyXcNi6ryaz768O1qCYvYm2g2Kp9YcQ9Zx+pbXIvLhudeYKR9xvtS1CWrV38w/3VRnE/Ns/I+y0kJHD2qRnVMeEDKW8Mbhf3v+bnISVrOswZmHLsn8c7oP+wG3vzqIzjJZnhA5kAQl1R1LZ9y9r8A0Y6sqGqsXub6LeBG5vdWr4o4u6nbyj+Ax5THN9aY+xXEftxc1se2kFk/8rYA/CME/r4xPziP9oXRt7RDPsR0qNt9DvfbaxqNHmxS0v+6VAcuSWQxlSbDG39stEZJ4Az/n6sA5o9IOVoboqmeS7xW0SpM+oKI5a2L4q8NTkuJCB3U2uVET64c0nUqrIsLWQcUTbXfhSmh3TvUhJ1d3luOGS7Dg1AU4VW5oTPgz4Wpqu7NeN+qOBeg3uuKC10YCCOqYxYVWEd5gzgjmiN01nfKSgoKCgoKCgoKCgoKCgoKCgoKCgoKCgoKCgoKCgoKCgoKCgoKCgoKCgoKCgofIlrsEAqpiiqCmorClmjBAvsN3JlsCqgNxZIvf0EkXexxLpZZoUfBlUvUie1lbxtFQx2kgL2z7gTC+Q0nItlOyEeueBIMBMg4BxRQC65OWGznAIse4nhnmXM+318heVmLL8y5j00VjDmjZDisUCC1NssDKq+WBYQ40NZK8mxJVgOSNqxAQvsVHYTlq+x9LDQjtvI9ZAR+CiWPuTaFeT3eYx5O05xnYC2FtovbguycM3touPHsFyH5QXSDrnYappEPyorfbVnq9EW+m/NVo7YCdqx1IKtxHZSWrg/7NnpgIVr5GgnMSBpK2xkHkb+LsJShqUlltewXI9lH7GNJZ3lYYEMTDdgeYqUIdbBRiwfYvmc6PQlUr4lmzCi8QD3wa1Ez1daOV9ad18Hx5ajfZMNAT5Lbp4V5LetWL7AMpPcTJBZ5A1y7ARRWC8yqBjSSUsEmE8GzD3k3GpyLAdLsgXvwEiMCLPk6xbacQW5PhHLY1jus+BZSOuEGfgBC+2vELWli4VrhOMp5MaA4+3JbyEysZVUPyVW+mrPVtU2+i+1lSN2KhFdL7aV1AOU2uohO3ZKtnFPyclOYkBG7u8t/L6FkITghS2xorP7GfNWswJ5WtLbMCx/MubU8wusjGGx1wn3yJvk/yeJni2dL63b0bHlaN9k5QGCe34WSzRRxgfkOKSx2SVSmNDJy4gnwJBjlgjwUSzbsHQiv3WwseSyNrDE7egpun4nUTQj+p0ny0NxnWqJwaUDy1I71ZLju7HchaWG/KaVia2ekegn1kpf7dnqJTv9d9ZOsZLrBVuJ7aS0oPedduzE2Lmn5GInWx5ga8a8P4stkhDrTExCBVb0BlsenMHyIOk7Y8UmYgIU6v5RRIDS86V1Ozq2HO2bbJ4rPU6I7yhxdUEBB4kLuxBLKVnurBC5zNGkc3PIsX1EIcKDdXCXp5Nly3LRQIN/x1ohwEbiSn9FllZFknbMIdcPJOW+R9qRRW6AdVjWS+rsTtx0afuvEbVluoVrhOPDscCmQq+S8xAZyHKw1ZfEuxD0w1rpqz1bvWij/2NdsBMrul5sqzKRnVpZ0Pu1duzE2Lin5GInS88A7yKkBf1bRchQvEz8H1kmWtJZjWgZeqsVHYC97ib9Z0Qem9Qm4iXw54ScfycemaXz8yV1T3dwbDnat+jm8pwwXzSbnsYywksPyA3NuP3BYitqJ4qgsxnM8vcS5n+GeA+exgriWXRppu0PFltRO1FQm1FQUFBQUFBQUFBQUFBQNA/AWxpnd2qHfUvvcqGuKYz5NXigoRxLQoDYyBN2k7udvW0vX9jX27r3tA3lpnMTIMJ7lQPGgeBKiN2ZR35LdXFwFUiUKpTtjRvKWrnifrhahhQQ3rBWBjYS989VG9mym1h3tvRTYGfweMI+7tw/3rSXJ+3rSv/s6d7X5chR5xcAsTtTyf/BCLcw5gh7iAWCTaYhGn8GGUgQFQ5xZ93I4IJ4sQUWOgXBqBCcC/FDlRaU+iyZXTZjGSMapFdJyrRUjnAObFa+iDG/aboaSwZjjq3SSryfqyy0U9yPTNJniEuaLSn/C1HbILbwDnJ8kuQ8YeP0u0l5/rSRuH/jbdgI2g8R/hDTtZJcAzGfCiyzCAGBrfZIBoNYd4KOp5DzphGbcqLzpfq15J1KbeSofW4UlSO1j/Re2c+YP4tbQ871pr3cse8gSTvF97E1e0n1XyQiLmv6t2QzqQ4FG7IWxhl4c/eTciCOb6AN/YvvP7no/AKgsUNEXt1qxhxkCh2OIf8XPAmpByh8vvKYpMxM0umZpMNSAhS+G1wuIUBpmZbKgXM2kv/Dp00QlH0rmSnFga5Cmy21U9yPlWSAAFncIyk/R9S2lcToAAgcTpO0AwCBnp39bCOpB2jNRuJjL4nskSPRz04bHqBQp9SmnUXnS/XLWFhBSNvpqH1CROVI7cNKzp1OzitlLn4H6y17uWPfiZJ2WrqPpfaS6v9KEQFa078lm0l1WCgqRzrOVovO3WVD/9L7Ty46v4AGMssIeJKwOhhlE5YkkRHKyY3JSW7gnZIy4esLCIhMlCjdEQIUl2mpHEvLuheJB8RbIUBpO8X9AMO1J793sjAwxQRYRv6/lxCgtB2LyeD1p43E/Uu3YSNLellC7APR/PPJb89JCNDSPWCPAMX6tUaA4rY4ah9GQoBi+7CScxPIvbQDS1cv28sd+0rbaes+FuxljwAt6d8aAYp1WChZAovH2RoLBGhL/ztlpvMLAJf2ZtHfwrICDPGWaLmyhxjmQaKUEeQ3cLc/FClDcOPXEbaHT5EiLbjeMCNvId4L/B1tocxlFsoRzskQlXkdOZeRLLH2kJtB2s4EUT+yiHs+kdyI4vJHidoG591GyptkpR2bGO98cO+MjaJF/fvSho2E9peRY/Dv80QPoYw5w0cDOVZC7AbHs0W6E+psRx5rTCU25UXnl0j0y1hor9TuvR20j7icMol9GMm5c8k9ByuJPC/byx37StspvY8t2atAov9p5JjWgv6kj6LENpPe41NE5UjHWSFZAoNn+YqFaxkrHCEXnTfBdOItNFf0IF6LP9GW8e5nP3K2kacflgeCvXxpX1f078o14nGWRCZKhjyvbY46bwJlMyW/xcRDDAb9ydVGUk+hOUAZQG1yRf/OXiMdZ+DNrSIe65AA5xgKCgoKCgoKCgoKCgoKCgoKCgoKGeH/AS8924+Oc1GZAAAAAElFTkSuQmCC"}]}

The area of a polygon (or any plane shape) can be evaluated by <a href = "https://www.mathworks.com/matlabcentral/cody/problems/179">Monte Carlo integration</a>. The process involves 'random' sampling of (x,y) points in the xy plane, and testing whether they fall inside or outside the shape. That is illustrated schematically below for a circle.<img src = "https://upload.wikimedia.org/wikipedia/commons/thumb/2/20/MonteCarloIntegrationCircle.svg/247px-MonteCarloIntegrationCircle.svg.png">^Steps:<ol><li>Define a 2D region within which to sample points. In the present problem you must choose the rectangular box (aligned to the axes) that bounds the specified polygon.</li><li>By a 'random' process generate coordinates of a point within the bounding box. In the present problem a basic scheme using the uniform pseudorandom distribution available in MATLAB must be employed. Other schemes in use include quasi-random sampling (for greater uniformity of coverage) and stratified sampling (with more attention near the edges of the polygon).</li><li>Determine <a href = "https://www.mathworks.com/matlabcentral/cody/problems/198">whether the sampled point is inside or outside the polygon</a>. Due to working in double precision, it is extremely unlikely to sample a point falling exactly on the edge of the polygon, and consequently if it does happen it doesn't really matter whether it's counted as inside or outside or null.</li><li>Repeat steps 2–3 <tt>N</tt> times.</li><li>Based upon the proportion of sampled points falling within the polygon, report the approximate area of the polygon.</li></ol>Inputs to your function will be <tt>N</tt>, the number of points to sample, and <tt>polygonX</tt> & <tt>polygonY</tt>, which together constitute <a href = "https://www.mathworks.com/matlabcentral/cody/problems/194">an ordered list of polygon vertices</a>. <tt>N</tt> will always be a positive integer. <tt>polygonX</tt> & <tt>polygonY</tt> will always describe a single, continuous outline; however, the polygon may be self-intersecting.For polygons that are self-intersecting, you must find the area of the enclosed point set. In the case of the cross-quadrilateral, it is treated as two simple triangles (cf. three triangles in the first figure below), and the clockwise/counterclockwise ordering of points is irrelevant. A self-intersecting pentagram (left & middle of the second figure below) will therefore have the same area as a concave decagon (right).<img src = "https://upload.wikimedia.org/wikipedia/commons/thumb/0/0f/Complex_polygon.svg/288px-Complex_polygon.svg.png"><img src = "https://upload.wikimedia.org/wikipedia/commons/thumb/8/80/Pentagram_interpretations.svg/320px-Pentagram_interpretations.svg.png">^EXAMPLE<pre> % Input
 N = 1000
 polygonX = [1 0 -1 0]
 polygonY = [0 1 0 -1]
 % Output
 A = 2.036</pre>Explanation: the above polygon is a square of side-length √2 centred on the origin, but rotated by 45°. The exact area is hence 2. As the value of <tt>N</tt> is moderate, the estimate by Monte Carlo integration is moderately accurate (an error of 0.036 in this example, corresponding to 509 of the 1000 sampled points falling within the polygon). 
Of course, if the code were run again then a slightly different value of <tt>A</tt> would probably be returned, such as <tt>A</tt>=1.992 (corresponding to 498 of the 1000 sampled points falling within the polygon), due to the random qualities of the technique.

The area of a polygon (or any plane shape) can be evaluated by <https://www.mathworks.com/matlabcentral/cody/problems/179 Monte Carlo integration>. The process involves 'random' sampling of (x,y) points in the xy plane, and testing whether they fall inside or outside the shape. That is illustrated schematically below for a circle.

<<https://upload.wikimedia.org/wikipedia/commons/thumb/2/20/MonteCarloIntegrationCircle.svg/247px-MonteCarloIntegrationCircle.svg.png>>

*Steps:*

# Define a 2D region within which to sample points. In the present problem you must choose the rectangular box (aligned to the axes) that bounds the specified polygon. 
# By a 'random' process generate coordinates of a point within the bounding box. In the present problem a basic scheme using the uniform pseudorandom distribution available in MATLAB must be employed. _Other schemes in use include quasi-random sampling (for greater uniformity of coverage) and stratified sampling (with more attention near the edges of the polygon)._
# Determine <https://www.mathworks.com/matlabcentral/cody/problems/198 whether the sampled point is inside or outside the polygon>. _Due to working in double precision, it is extremely unlikely to sample a point falling exactly on the edge of the polygon, and consequently if it does happen it doesn't really matter whether it's counted as inside or outside or null._ 
# Repeat steps 2–3 |N| times. 
# Based upon the proportion of sampled points falling within the polygon, report the approximate area of the polygon.

For polygons that are self-intersecting, you must find the area of the enclosed point set. In the case of the cross-quadrilateral, it is treated as two simple triangles (cf. three triangles in the first figure below), and the clockwise/counterclockwise ordering of points is irrelevant.  A self-intersecting pentagram (left & middle of the second figure below) will therefore have the same area as a concave decagon (right).

<<https://upload.wikimedia.org/wikipedia/commons/thumb/0/0f/Complex_polygon.svg/288px-Complex_polygon.svg.png>>

<<https://upload.wikimedia.org/wikipedia/commons/thumb/8/80/Pentagram_interpretations.svg/320px-Pentagram_interpretations.svg.png>>

EXAMPLE

% Input
 N = 1000
 polygonX = [1 0 -1  0]
 polygonY = [0 1  0 -1]
 % Output
 A = 2.036

Explanation:  the above polygon is a square of side-length √2 centred on the origin, but rotated by 45°.  The exact area is hence 2.  As the value of |N| is moderate, the estimate by Monte Carlo integration is moderately accurate (an error of 0.036 in this example, corresponding to 509 of the 1000 sampled points falling within the polygon).  
Of course, if the code were run again then a slightly different value of |A| would probably be returned, such as |A|=1.992 (corresponding to 498 of the 1000 sampled points falling within the polygon), due to the random qualities of the technique.

Solve

Solution Stats

9 Solutions
3 Solvers

Last Solution submitted on Jan 15, 2021

Last 200 Solutions

Solution Comments

Solution 1576522

1 Comment

1 Comment

David Verrelli on 24 Jun 2019

An elegant solution.

Problem Recent Solvers3

Problem Tags

area integration monte carlo polygon polygrams random self-intersecting polygons

Community Treasure Hunt

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

Start Hunting!