7 October 2025 1:04:55.622 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 2000 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.9969 1.000 0.3108E-02 21 -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.9593 0.9704 0.1115E-01 355 -1.800 -0.3000 0.9012 0.9048 0.3631E-02 862 -1.800 -0.2000 0.8556 0.8607 0.5150E-02 1266 -1.800 -0.1000 0.8296 0.8353 0.5637E-02 1477 -1.800 0.000 0.8184 0.8270 0.8561E-02 1586 -1.800 0.1000 0.8370 0.8353 0.1768E-02 1416 -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.9908 1.000 0.9163E-02 66 -1.600 -0.5000 0.8873 0.8958 0.8574E-02 1003 -1.600 -0.4000 0.8054 0.8187 0.1335E-01 1784 -1.600 -0.3000 0.7611 0.7634 0.2272E-02 2093 -1.600 -0.2000 0.7136 0.7261 0.1250E-01 2699 -1.600 -0.1000 0.6991 0.7047 0.5630E-02 2894 -1.600 0.000 0.7035 0.6977 0.5857E-02 2794 -1.600 0.1000 0.6974 0.7047 0.7319E-02 2890 -1.600 0.2000 0.7167 0.7261 0.9401E-02 2695 -1.400 -0.7000 0.9687 0.9802 0.1149E-01 267 -1.400 -0.6000 0.8475 0.8607 0.1321E-01 1299 -1.400 -0.5000 0.7557 0.7711 0.1537E-01 2205 -1.400 -0.4000 0.6957 0.7047 0.8981E-02 2902 -1.400 -0.3000 0.6479 0.6570 0.9131E-02 3482 -1.400 -0.2000 0.6224 0.6250 0.2614E-02 3800 -1.400 -0.1000 0.5997 0.6065 0.6874E-02 4080 -1.400 0.000 0.5926 0.6005 0.7938E-02 4160 -1.400 0.1000 0.6013 0.6065 0.5225E-02 4063 -1.400 0.2000 0.6231 0.6250 0.1910E-02 3747 -1.400 0.3000 0.6473 0.6570 0.9734E-02 3484 -1.200 -0.6000 0.7502 0.7558 0.5625E-02 2261 -1.200 -0.5000 0.6780 0.6771 0.9194E-03 3097 -1.200 -0.4000 0.6215 0.6188 0.2729E-02 3734 -1.200 -0.3000 0.5649 0.5769 0.1210E-01 4540 -1.200 -0.2000 0.5426 0.5488 0.6206E-02 4936 -1.200 -0.1000 0.5373 0.5326 0.4722E-02 5048 -1.200 0.000 0.5301 0.5273 0.2837E-02 5032 -1.200 0.1000 0.5412 0.5326 0.8625E-02 4828 -1.200 0.2000 0.5433 0.5488 0.5524E-02 4871 -1.200 0.3000 0.5799 0.5769 0.2982E-02 4316 -1.200 0.4000 0.6133 0.6188 0.5448E-02 3873 -1.000 -0.5000 0.6003 0.6065 0.6254E-02 4084 -1.000 -0.4000 0.5502 0.5543 0.4107E-02 4739 -1.000 -0.3000 0.4997 0.5169 0.1715E-01 5590 -1.000 -0.2000 0.4787 0.4916 0.1297E-01 5991 -1.000 -0.1000 0.4690 0.4771 0.8143E-02 6111 -1.000 0.000 0.4657 0.4724 0.6707E-02 6278 -1.000 0.1000 0.4705 0.4771 0.6566E-02 6019 -1.000 0.2000 0.4886 0.4916 0.3044E-02 5780 -1.000 0.3000 0.5001 0.5169 0.1673E-01 5644 -1.000 0.4000 0.5363 0.5543 0.1800E-01 5070 -1.000 0.5000 0.6024 0.6065 0.4116E-02 4077 -0.8000 -0.4000 0.5029 0.5066 0.3749E-02 5542 -0.8000 -0.3000 0.4757 0.4724 0.3297E-02 5966 -0.8000 -0.2000 0.4463 0.4493 0.3004E-02 6386 -0.8000 -0.1000 0.4375 0.4360 0.1405E-02 6576 -0.8000 0.000 0.4237 0.4317 0.7983E-02 6837 -0.8000 0.1000 0.4259 0.4360 0.1013E-01 6788 -0.8000 0.2000 0.4379 0.4493 0.1143E-01 6525 -0.8000 0.3000 0.4694 0.4724 0.2973E-02 5984 -0.8000 0.4000 0.4979 0.5066 0.8733E-02 5582 -0.8000 0.5000 0.5504 0.5543 0.3928E-02 4741 -0.8000 0.6000 0.5979 0.6188 0.2087E-01 4192 -0.6000 -0.3000 0.4345 0.4404 0.5935E-02 6664 -0.6000 -0.2000 0.4152 0.4190 0.3764E-02 7000 -0.6000 -0.1000 0.4020 0.4066 0.4596E-02 7275 -0.6000 0.000 0.3922 0.4025 0.1030E-01 7496 -0.6000 0.1000 0.4037 0.4066 0.2915E-02 7238 -0.6000 0.2000 0.4065 0.4190 0.1242E-01 7201 -0.6000 0.3000 0.4437 0.4404 0.3281E-02 6491 -0.6000 0.4000 0.4604 0.4724 0.1196E-01 6244 -0.6000 0.5000 0.5251 0.5169 0.8259E-02 5122 -0.6000 0.6000 0.5713 0.5769 0.5614E-02 4368 -0.6000 0.7000 0.6531 0.6570 0.3914E-02 3367 -0.4000 -0.2000 0.3920 0.3985 0.6547E-02 7435 -0.4000 -0.1000 0.3863 0.3867 0.4285E-03 7601 -0.4000 0.000 0.3800 0.3829 0.2862E-02 7783 -0.4000 0.1000 0.3790 0.3867 0.7758E-02 7848 -0.4000 0.2000 0.4025 0.3985 0.3945E-02 7152 -0.4000 0.3000 0.4124 0.4190 0.6550E-02 7007 -0.4000 0.4000 0.4428 0.4493 0.6545E-02 6668 -0.4000 0.5000 0.4822 0.4916 0.9403E-02 5799 -0.4000 0.6000 0.5483 0.5488 0.5514E-03 4811 -0.4000 0.7000 0.6263 0.6250 0.1263E-02 3658 -0.4000 0.8000 0.7258 0.7261 0.3280E-03 2505 -0.2000 -0.1000 0.3684 0.3753 0.6957E-02 8030 -0.2000 0.000 0.3707 0.3716 0.8804E-03 8033 -0.2000 0.1000 0.3751 0.3753 0.1937E-03 7838 -0.2000 0.2000 0.3847 0.3867 0.2080E-02 7678 -0.2000 0.3000 0.4065 0.4066 0.2863E-04 7194 -0.2000 0.4000 0.4244 0.4360 0.1167E-01 6877 -0.2000 0.5000 0.4783 0.4771 0.1227E-02 5805 -0.2000 0.6000 0.5230 0.5326 0.9634E-02 5176 -0.2000 0.7000 0.5881 0.6065 0.1838E-01 4343 -0.2000 0.8000 0.6911 0.7047 0.1356E-01 3078 -0.2000 0.9000 0.8243 0.8353 0.1097E-01 1550 0.000 0.000 0.3621 0.3679 0.5756E-02 8122 0.000 0.1000 0.3735 0.3716 0.1924E-02 7854 0.000 0.2000 0.3812 0.3829 0.1663E-02 7728 0.000 0.3000 0.3959 0.4025 0.6593E-02 7352 0.000 0.4000 0.4280 0.4317 0.3712E-02 6801 0.000 0.5000 0.4618 0.4724 0.1052E-01 6095 0.000 0.6000 0.5132 0.5273 0.1414E-01 5382 0.000 0.7000 0.5974 0.6005 0.3101E-02 4084 0.000 0.8000 0.6944 0.6977 0.3297E-02 2857 0.000 0.9000 0.8173 0.8270 0.9704E-02 1570 0.000 1.000 0.9921 1.000 0.7900E-02 62 0.2000 0.1000 0.3675 0.3753 0.7861E-02 8144 0.2000 0.2000 0.3785 0.3867 0.8224E-02 7689 0.2000 0.3000 0.3977 0.4066 0.8892E-02 7386 0.2000 0.4000 0.4262 0.4360 0.9864E-02 6874 0.2000 0.5000 0.4724 0.4771 0.4749E-02 6025 0.2000 0.6000 0.5195 0.5326 0.1307E-01 5239 0.2000 0.7000 0.5968 0.6065 0.9691E-02 4087 0.2000 0.8000 0.6998 0.7047 0.4908E-02 2800 0.2000 0.9000 0.8178 0.8353 0.1747E-01 1618 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.3983 0.3985 0.2683E-03 7322 0.4000 0.3000 0.4150 0.4190 0.3997E-02 6973 0.4000 0.4000 0.4455 0.4493 0.3850E-02 6421 0.4000 0.5000 0.4867 0.4916 0.4979E-02 5800 0.4000 0.6000 0.5516 0.5488 0.2785E-02 4644 0.4000 0.7000 0.6186 0.6250 0.6441E-02 3817 0.4000 0.8000 0.7304 0.7261 0.4219E-02 2358 0.4000 0.9000 0.8519 0.8607 0.8832E-02 1339 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.4351 0.4404 0.5376E-02 6637 0.6000 0.4000 0.4629 0.4724 0.9446E-02 6273 0.6000 0.5000 0.5178 0.5169 0.9117E-03 5214 0.6000 0.6000 0.5617 0.5769 0.1522E-01 4566 0.6000 0.7000 0.6468 0.6570 0.1029E-01 3588 0.6000 0.8000 0.7567 0.7634 0.6705E-02 2224 0.6000 0.9000 0.8988 0.9048 0.6052E-02 822 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.5088 0.5066 0.2227E-02 5355 0.8000 0.5000 0.5477 0.5543 0.6606E-02 4877 0.8000 0.6000 0.6150 0.6188 0.3805E-02 3934 0.8000 0.7000 0.6988 0.7047 0.5847E-02 2840 0.8000 0.8000 0.8090 0.8187 0.9761E-02 1677 0.8000 0.9000 0.9598 0.9704 0.1061E-01 339 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.5958 0.6065 0.1073E-01 4083 1.000 0.6000 0.6791 0.6771 0.2065E-02 3095 1.000 0.7000 0.7713 0.7711 0.2061E-03 2067 1.000 0.8000 0.8842 0.8958 0.1163E-01 1042 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.7482 0.7558 0.7615E-02 2276 1.200 0.7000 0.8494 0.8607 0.1133E-01 1306 1.200 0.8000 0.9885 1.000 0.1149E-01 104 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.9735 0.9802 0.6719E-02 218 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.818503E-02 feynman_kac_2d(): Normal end of execution. 7 October 2025 1:05:59.244 PM