7 October 2023 4:06:16.868 PM feynman_kac_2d(): FORTRAN77 version. Program parameters: The calculation takes place inside a 2D ellipse. A rectangular grid of points will be defined. The solution will be estimated for those grid points that lie inside the ellipse. Each solution will be estimated by computing 10000 trajectories from the point to the boundary. (X/A)^2 + (Y/B)^2 = 1 The ellipsoid parameters A, B are: A = 2.00000 B = 1.00000 Stepsize H = 0.100000E-03 X coordinate marked by 21 points Y coordinate marked by 11 points X Y W Approx W Exact Error Ave Steps -2.000 -1.000 1.000 1.000 0.000 0 -2.000 -0.9000 1.000 1.000 0.000 0 -2.000 -0.8000 1.000 1.000 0.000 0 -2.000 -0.7000 1.000 1.000 0.000 0 -2.000 -0.6000 1.000 1.000 0.000 0 -2.000 -0.5000 1.000 1.000 0.000 0 -2.000 -0.4000 1.000 1.000 0.000 0 -2.000 -0.3000 1.000 1.000 0.000 0 -2.000 -0.2000 1.000 1.000 0.000 0 -2.000 -0.1000 1.000 1.000 0.000 0 -2.000 0.000 0.9960 1.000 0.3971E-02 34 -1.800 -0.9000 1.000 1.000 0.000 0 -1.800 -0.8000 1.000 1.000 0.000 0 -1.800 -0.7000 1.000 1.000 0.000 0 -1.800 -0.6000 1.000 1.000 0.000 0 -1.800 -0.5000 1.000 1.000 0.000 0 -1.800 -0.4000 0.9626 0.9704 0.7855E-02 301 -1.800 -0.3000 0.8959 0.9048 0.8983E-02 884 -1.800 -0.2000 0.8585 0.8607 0.2251E-02 1207 -1.800 -0.1000 0.8283 0.8353 0.6941E-02 1549 -1.800 0.000 0.8225 0.8270 0.4450E-02 1591 -1.800 0.1000 0.8310 0.8353 0.4299E-02 1488 -1.600 -0.8000 1.000 1.000 0.000 0 -1.600 -0.7000 1.000 1.000 0.000 0 -1.600 -0.6000 0.9925 1.000 0.7491E-02 60 -1.600 -0.5000 0.8897 0.8958 0.6161E-02 941 -1.600 -0.4000 0.8116 0.8187 0.7151E-02 1671 -1.600 -0.3000 0.7581 0.7634 0.5327E-02 2200 -1.600 -0.2000 0.7196 0.7261 0.6582E-02 2607 -1.600 -0.1000 0.6994 0.7047 0.5281E-02 2874 -1.600 0.000 0.6912 0.6977 0.6509E-02 2977 -1.600 0.1000 0.7012 0.7047 0.3508E-02 2830 -1.600 0.2000 0.7199 0.7261 0.6214E-02 2597 -1.400 -0.7000 0.9692 0.9802 0.1103E-01 250 -1.400 -0.6000 0.8518 0.8607 0.8877E-02 1278 -1.400 -0.5000 0.7640 0.7711 0.7011E-02 2136 -1.400 -0.4000 0.7030 0.7047 0.1735E-02 2803 -1.400 -0.3000 0.6544 0.6570 0.2648E-02 3403 -1.400 -0.2000 0.6221 0.6250 0.2909E-02 3776 -1.400 -0.1000 0.6023 0.6065 0.4281E-02 4072 -1.400 0.000 0.5968 0.6005 0.3649E-02 4109 -1.400 0.1000 0.6019 0.6065 0.4674E-02 4028 -1.400 0.2000 0.6190 0.6250 0.5999E-02 3836 -1.400 0.3000 0.6469 0.6570 0.1010E-01 3512 -1.200 -0.6000 0.7516 0.7558 0.4194E-02 2265 -1.200 -0.5000 0.6737 0.6771 0.3357E-02 3162 -1.200 -0.4000 0.6140 0.6188 0.4755E-02 3865 -1.200 -0.3000 0.5717 0.5769 0.5215E-02 4513 -1.200 -0.2000 0.5466 0.5488 0.2214E-02 4824 -1.200 -0.1000 0.5313 0.5326 0.1302E-02 5001 -1.200 0.000 0.5198 0.5273 0.7470E-02 5283 -1.200 0.1000 0.5287 0.5326 0.3859E-02 5073 -1.200 0.2000 0.5432 0.5488 0.5632E-02 4891 -1.200 0.3000 0.5741 0.5769 0.2818E-02 4448 -1.200 0.4000 0.6137 0.6188 0.5053E-02 3917 -1.000 -0.5000 0.6045 0.6065 0.1995E-02 3999 -1.000 -0.4000 0.5495 0.5543 0.4787E-02 4813 -1.000 -0.3000 0.5127 0.5169 0.4146E-02 5311 -1.000 -0.2000 0.4882 0.4916 0.3464E-02 5748 -1.000 -0.1000 0.4750 0.4771 0.2139E-02 5936 -1.000 0.000 0.4687 0.4724 0.3655E-02 6082 -1.000 0.1000 0.4734 0.4771 0.3724E-02 5935 -1.000 0.2000 0.4872 0.4916 0.4441E-02 5782 -1.000 0.3000 0.5153 0.5169 0.1557E-02 5282 -1.000 0.4000 0.5495 0.5543 0.4782E-02 4773 -1.000 0.5000 0.6047 0.6065 0.1820E-02 4027 -0.8000 -0.4000 0.5024 0.5066 0.4266E-02 5476 -0.8000 -0.3000 0.4698 0.4724 0.2587E-02 6019 -0.8000 -0.2000 0.4482 0.4493 0.1148E-02 6415 -0.8000 -0.1000 0.4324 0.4360 0.3617E-02 6716 -0.8000 0.000 0.4288 0.4317 0.2945E-02 6692 -0.8000 0.1000 0.4334 0.4360 0.2661E-02 6659 -0.8000 0.2000 0.4468 0.4493 0.2501E-02 6427 -0.8000 0.3000 0.4653 0.4724 0.7029E-02 6120 -0.8000 0.4000 0.5056 0.5066 0.1021E-02 5472 -0.8000 0.5000 0.5448 0.5543 0.9529E-02 4786 -0.8000 0.6000 0.6108 0.6188 0.7981E-02 3948 -0.6000 -0.3000 0.4338 0.4404 0.6655E-02 6665 -0.6000 -0.2000 0.4133 0.4190 0.5638E-02 7054 -0.6000 -0.1000 0.4083 0.4066 0.1731E-02 7144 -0.6000 0.000 0.3998 0.4025 0.2759E-02 7341 -0.6000 0.1000 0.4027 0.4066 0.3879E-02 7245 -0.6000 0.2000 0.4124 0.4190 0.6535E-02 7055 -0.6000 0.3000 0.4368 0.4404 0.3681E-02 6618 -0.6000 0.4000 0.4674 0.4724 0.4974E-02 6072 -0.6000 0.5000 0.5068 0.5169 0.1005E-01 5434 -0.6000 0.6000 0.5716 0.5769 0.5330E-02 4474 -0.6000 0.7000 0.6533 0.6570 0.3729E-02 3377 -0.4000 -0.2000 0.3974 0.3985 0.1106E-02 7293 -0.4000 -0.1000 0.3826 0.3867 0.4116E-02 7671 -0.4000 0.000 0.3792 0.3829 0.3678E-02 7785 -0.4000 0.1000 0.3810 0.3867 0.5723E-02 7721 -0.4000 0.2000 0.3956 0.3985 0.2890E-02 7382 -0.4000 0.3000 0.4179 0.4190 0.1080E-02 6916 -0.4000 0.4000 0.4428 0.4493 0.6571E-02 6518 -0.4000 0.5000 0.4860 0.4916 0.5665E-02 5762 -0.4000 0.6000 0.5411 0.5488 0.7663E-02 4955 -0.4000 0.7000 0.6173 0.6250 0.7738E-02 3906 -0.4000 0.8000 0.7185 0.7261 0.7608E-02 2619 -0.2000 -0.1000 0.3706 0.3753 0.4691E-02 7930 -0.2000 0.000 0.3718 0.3716 0.2213E-03 7892 -0.2000 0.1000 0.3701 0.3753 0.5239E-02 7917 -0.2000 0.2000 0.3870 0.3867 0.3011E-03 7627 -0.2000 0.3000 0.4023 0.4066 0.4252E-02 7290 -0.2000 0.4000 0.4297 0.4360 0.6313E-02 6694 -0.2000 0.5000 0.4729 0.4771 0.4195E-02 5990 -0.2000 0.6000 0.5256 0.5326 0.6983E-02 5137 -0.2000 0.7000 0.5984 0.6065 0.8139E-02 4071 -0.2000 0.8000 0.6917 0.7047 0.1297E-01 2944 -0.2000 0.9000 0.8244 0.8353 0.1090E-01 1567 0.000 0.000 0.3644 0.3679 0.3486E-02 8044 0.000 0.1000 0.3694 0.3716 0.2191E-02 7938 0.000 0.2000 0.3795 0.3829 0.3360E-02 7790 0.000 0.3000 0.3975 0.4025 0.5024E-02 7385 0.000 0.4000 0.4283 0.4317 0.3455E-02 6755 0.000 0.5000 0.4636 0.4724 0.8769E-02 6179 0.000 0.6000 0.5287 0.5273 0.1415E-02 4997 0.000 0.7000 0.5909 0.6005 0.9566E-02 4203 0.000 0.8000 0.6843 0.6977 0.1335E-01 3024 0.000 0.9000 0.8101 0.8270 0.1681E-01 1702 0.000 1.000 0.9928 1.000 0.7180E-02 59 0.2000 0.1000 0.3723 0.3753 0.2991E-02 7896 0.2000 0.2000 0.3840 0.3867 0.2754E-02 7654 0.2000 0.3000 0.4016 0.4066 0.4929E-02 7240 0.2000 0.4000 0.4322 0.4360 0.3817E-02 6706 0.2000 0.5000 0.4752 0.4771 0.1907E-02 5968 0.2000 0.6000 0.5233 0.5326 0.9291E-02 5180 0.2000 0.7000 0.5961 0.6065 0.1046E-01 4140 0.2000 0.8000 0.6944 0.7047 0.1026E-01 2971 0.2000 0.9000 0.8274 0.8353 0.7893E-02 1496 0.2000 1.000 1.000 1.000 0.000 0 0.2000 1.100 1.000 1.000 0.000 0 0.4000 0.2000 0.3945 0.3985 0.4049E-02 7463 0.4000 0.3000 0.4150 0.4190 0.3988E-02 7009 0.4000 0.4000 0.4436 0.4493 0.5774E-02 6517 0.4000 0.5000 0.4867 0.4916 0.4971E-02 5738 0.4000 0.6000 0.5399 0.5488 0.8883E-02 4928 0.4000 0.7000 0.6128 0.6250 0.1218E-01 3926 0.4000 0.8000 0.7126 0.7261 0.1353E-01 2711 0.4000 0.9000 0.8472 0.8607 0.1349E-01 1329 0.4000 1.000 1.000 1.000 0.000 0 0.4000 1.100 1.000 1.000 0.000 0 0.4000 1.200 1.000 1.000 0.000 0 0.6000 0.3000 0.4379 0.4404 0.2526E-02 6560 0.6000 0.4000 0.4674 0.4724 0.4949E-02 6064 0.6000 0.5000 0.5168 0.5169 0.7646E-04 5248 0.6000 0.6000 0.5704 0.5769 0.6510E-02 4517 0.6000 0.7000 0.6523 0.6570 0.4751E-02 3402 0.6000 0.8000 0.7498 0.7634 0.1359E-01 2324 0.6000 0.9000 0.8927 0.9048 0.1214E-01 932 0.6000 1.000 1.000 1.000 0.000 0 0.6000 1.100 1.000 1.000 0.000 0 0.6000 1.200 1.000 1.000 0.000 0 0.6000 1.300 1.000 1.000 0.000 0 0.8000 0.4000 0.5052 0.5066 0.1369E-02 5454 0.8000 0.5000 0.5497 0.5543 0.4639E-02 4776 0.8000 0.6000 0.6137 0.6188 0.5118E-02 3879 0.8000 0.7000 0.7015 0.7047 0.3183E-02 2808 0.8000 0.8000 0.8095 0.8187 0.9281E-02 1664 0.8000 0.9000 0.9589 0.9704 0.1150E-01 335 0.8000 1.000 1.000 1.000 0.000 0 0.8000 1.100 1.000 1.000 0.000 0 0.8000 1.200 1.000 1.000 0.000 0 0.8000 1.300 1.000 1.000 0.000 0 0.8000 1.400 1.000 1.000 0.000 0 1.000 0.5000 0.5993 0.6065 0.7196E-02 4104 1.000 0.6000 0.6687 0.6771 0.8322E-02 3232 1.000 0.7000 0.7593 0.7711 0.1173E-01 2225 1.000 0.8000 0.8843 0.8958 0.1149E-01 995 1.000 0.9000 1.000 1.000 0.000 0 1.000 1.000 1.000 1.000 0.000 0 1.000 1.100 1.000 1.000 0.000 0 1.000 1.200 1.000 1.000 0.000 0 1.000 1.300 1.000 1.000 0.000 0 1.000 1.400 1.000 1.000 0.000 0 1.000 1.500 1.000 1.000 0.000 0 1.200 0.6000 0.7493 0.7558 0.6461E-02 2285 1.200 0.7000 0.8534 0.8607 0.7320E-02 1240 1.200 0.8000 0.9899 1.000 0.1010E-01 82 1.200 0.9000 1.000 1.000 0.000 0 1.200 1.000 1.000 1.000 0.000 0 1.200 1.100 1.000 1.000 0.000 0 1.200 1.200 1.000 1.000 0.000 0 1.200 1.300 1.000 1.000 0.000 0 1.200 1.400 1.000 1.000 0.000 0 1.200 1.500 1.000 1.000 0.000 0 1.200 1.600 1.000 1.000 0.000 0 1.400 0.7000 0.9681 0.9802 0.1213E-01 272 1.400 0.8000 1.000 1.000 0.000 0 1.400 0.9000 1.000 1.000 0.000 0 1.400 1.000 1.000 1.000 0.000 0 1.400 1.100 1.000 1.000 0.000 0 1.400 1.200 1.000 1.000 0.000 0 1.400 1.300 1.000 1.000 0.000 0 1.400 1.400 1.000 1.000 0.000 0 1.400 1.500 1.000 1.000 0.000 0 1.400 1.600 1.000 1.000 0.000 0 1.400 1.700 1.000 1.000 0.000 0 1.600 0.8000 1.000 1.000 0.000 0 1.600 0.9000 1.000 1.000 0.000 0 1.600 1.000 1.000 1.000 0.000 0 1.600 1.100 1.000 1.000 0.000 0 1.600 1.200 1.000 1.000 0.000 0 1.600 1.300 1.000 1.000 0.000 0 1.600 1.400 1.000 1.000 0.000 0 1.600 1.500 1.000 1.000 0.000 0 1.600 1.600 1.000 1.000 0.000 0 1.600 1.700 1.000 1.000 0.000 0 1.600 1.800 1.000 1.000 0.000 0 1.800 0.9000 1.000 1.000 0.000 0 1.800 1.000 1.000 1.000 0.000 0 1.800 1.100 1.000 1.000 0.000 0 1.800 1.200 1.000 1.000 0.000 0 1.800 1.300 1.000 1.000 0.000 0 1.800 1.400 1.000 1.000 0.000 0 1.800 1.500 1.000 1.000 0.000 0 1.800 1.600 1.000 1.000 0.000 0 1.800 1.700 1.000 1.000 0.000 0 1.800 1.800 1.000 1.000 0.000 0 1.800 1.900 1.000 1.000 0.000 0 2.000 1.000 1.000 1.000 0.000 0 2.000 1.100 1.000 1.000 0.000 0 2.000 1.200 1.000 1.000 0.000 0 2.000 1.300 1.000 1.000 0.000 0 2.000 1.400 1.000 1.000 0.000 0 2.000 1.500 1.000 1.000 0.000 0 2.000 1.600 1.000 1.000 0.000 0 2.000 1.700 1.000 1.000 0.000 0 2.000 1.800 1.000 1.000 0.000 0 2.000 1.900 1.000 1.000 0.000 0 2.000 2.000 1.000 1.000 0.000 0 RMS absolute error in solution = 0.657801E-02 feynman_kac_2d(): Normal end of execution. 7 October 2023 4:14:13.321 PM