matlab-validation.m

af_pts_0 = [1.       , 0.99891  , 0.99571  , 0.99034  , 0.98291  , 0.97344  ,       0.96192  , 0.94848  , 0.93299  , 0.91571  , 0.89644  , 0.8759   ,       0.8535   , 0.8297   , 0.80435  , 0.77773  , 0.7499501, 0.72115  ,       0.69133  , 0.66074  , 0.62938  , 0.5975   , 0.56525  , 0.53274  ,       0.49997  , 0.46733  , 0.43469  , 0.40243  , 0.37056  , 0.33933  ,       0.30861  , 0.27891  , 0.24998  , 0.22221  , 0.19558  , 0.17037  ,       0.14643  , 0.12403  , 0.1033   , 0.08422  , 0.06694  , 0.05158  ,       0.03802  , 0.0265   , 0.01702  , 0.0096   , 0.00422  , 0.00102  ,       0.       , 0.00102  , 0.00422  , 0.0096   , 0.01702  , 0.0265   ,       0.03802  , 0.05158  , 0.06694  , 0.08422  , 0.1033   , 0.12403  ,       0.14643  , 0.17037  , 0.19558  , 0.22221  , 0.24998  , 0.27891  ,       0.30861  , 0.33933  , 0.37056  , 0.40243  , 0.43469  , 0.46733  ,       0.49997  , 0.53274  , 0.56525  , 0.5975   , 0.62938  , 0.66074  ,       0.69133  , 0.72115  , 0.7499501, 0.77773  , 0.80435  , 0.8297   ,       0.8535   , 0.8759   , 0.89644  , 0.91571  , 0.93299  , 0.94848  ,       0.96192  , 0.97344  , 0.98291  , 0.99034  , 0.99571  , 0.99891  ,       1.       ];

af_pts_1 = [ 0.     ,  0.00026,  0.00107,  0.00235,  0.00415,  0.00643,        0.00918,  0.01244,  0.01627,  0.0206 ,  0.02551,  0.03094,        0.03695,  0.04338,  0.05032,  0.05754,  0.06512,  0.07274,        0.08058,  0.08828,  0.09602,  0.10334,  0.11046,  0.11695,        0.12294,  0.12776,  0.13176,  0.13385,  0.13444,  0.13306,        0.13099,  0.12732,  0.12347,  0.11794,  0.11262,  0.10557,        0.09914,  0.09093,  0.08384,  0.07483,  0.06739,  0.05781,        0.05044,  0.04047,  0.03369,  0.026  ,  0.01791,  0.0085 ,        0.     , -0.0035 , -0.008  , -0.01293, -0.01603, -0.02164,       -0.02484, -0.02993, -0.03315, -0.0379 , -0.04104, -0.0454 ,       -0.04826, -0.05208, -0.05447, -0.05747, -0.05906, -0.06103,       -0.06162, -0.06225, -0.06143, -0.06026, -0.05699, -0.05322,       -0.04826, -0.04337, -0.03793, -0.03278, -0.02744, -0.02259,       -0.01783, -0.01362, -0.0097 , -0.00633, -0.00333, -0.00093,        0.00111,  0.00258,  0.00363,  0.00425,  0.0045 ,  0.00438,        0.00395,  0.00324,  0.00236,  0.00147,  0.0007 ,  0.00019,        0.     ];

out_0_6 = [1.00000000e+00, 9.79894292e-01, 9.59798684e-01, 9.39715244e-01,       9.19646043e-01, 8.99593149e-01, 8.79558631e-01, 8.59544558e-01,       8.39552496e-01, 8.19578815e-01, 7.99616827e-01, 7.79659804e-01,       7.59701018e-01, 7.39733740e-01, 7.19751244e-01, 6.99746800e-01,       6.79713680e-01, 6.59645157e-01, 6.39534502e-01, 6.19375144e-01,       5.99164139e-01, 5.78902202e-01, 5.58590205e-01, 5.38229022e-01,       5.17819528e-01, 4.97362594e-01, 4.76859096e-01, 4.56309907e-01,       4.35715900e-01, 4.15077948e-01, 3.94396927e-01, 3.73673708e-01,       3.52911833e-01, 3.32131594e-01, 3.11359773e-01, 2.90623168e-01,       2.69948577e-01, 2.49362796e-01, 2.28892623e-01, 2.08564854e-01,       1.88406288e-01, 1.68443722e-01, 1.48703952e-01, 1.29214081e-01,       1.10002932e-01, 9.10999375e-02, 7.25345281e-02, 5.43361357e-02,       3.65341917e-02, 1.91581278e-02, 3.71665063e-03, 3.25760752e-04,       1.54611300e-02, 3.52478243e-02, 5.51436758e-02, 7.51440904e-02,       9.52463170e-02, 1.15447605e-01, 1.35745203e-01, 1.56136193e-01,       1.76613309e-01, 1.97164637e-01, 2.17778037e-01, 2.38441370e-01,       2.59142497e-01, 2.79869276e-01, 3.00609568e-01, 3.21351234e-01,       3.42082134e-01, 3.62790127e-01, 3.83463695e-01, 4.04099797e-01,       4.24701543e-01, 4.45272179e-01, 4.65814949e-01, 4.86333099e-01,       5.06829875e-01, 5.27308523e-01, 5.47772288e-01, 5.68224416e-01,       5.88668152e-01, 6.09106742e-01, 6.29543431e-01, 6.49981606e-01,       6.70425201e-01, 6.90878286e-01, 7.11344931e-01, 7.31829205e-01,       7.52335177e-01, 7.72866918e-01, 7.93428498e-01, 8.14023984e-01,       8.34657449e-01, 8.55332957e-01, 8.76048898e-01, 8.96786077e-01,       9.17521906e-01, 9.38233800e-01, 9.58899172e-01, 9.79495434e-01,       1.00000000e+00];

out_1_6 = [ 6.58362934e-18,  4.88680971e-03,  9.82961593e-03,  1.48278433e-02,        1.98809165e-02,  2.49882601e-02,  3.01492989e-02,  3.53634574e-02,        4.06294368e-02,  4.59384803e-02,  5.12774365e-02,  5.66330976e-02,        6.19922554e-02,  6.73417021e-02,  7.26682297e-02,  7.79586303e-02,        8.31996958e-02,  8.83782183e-02,  9.34809899e-02,  9.84940183e-02,        1.03385054e-01,  1.08103468e-01,  1.12597821e-01,  1.16816676e-01,        1.20708594e-01,  1.24222136e-01,  1.27305866e-01,  1.29908344e-01,        1.31978133e-01,  1.33463793e-01,  1.34313887e-01,  1.34476977e-01,        1.33912973e-01,  1.32653082e-01,  1.30756135e-01,  1.28281024e-01,        1.25286641e-01,  1.21831878e-01,  1.17975626e-01,  1.13776777e-01,        1.09294223e-01,  1.04586854e-01,  9.97135637e-02,  9.46694985e-02,        8.90886778e-02,  8.24780949e-02,  7.43446027e-02,  6.41950542e-02,        5.15363023e-02,  3.58752001e-02,  1.72874212e-02, -8.41460340e-04,       -1.50685973e-02, -2.57553004e-02, -3.36003310e-02, -3.91438379e-02,       -4.29259054e-02, -4.54866176e-02, -4.73660589e-02, -4.90935273e-02,       -5.09179849e-02, -5.27888727e-02, -5.46410037e-02, -5.64091908e-02,       -5.80282468e-02, -5.94329845e-02, -6.05582170e-02, -6.13387570e-02,       -6.17094174e-02, -6.16050112e-02, -6.09651647e-02, -5.97953797e-02,       -5.81488958e-02, -5.60800208e-02, -5.36430624e-02, -5.08923282e-02,       -4.78821261e-02, -4.46667636e-02, -4.13005485e-02, -3.78377884e-02,       -3.43327912e-02, -3.08398644e-02, -2.74133158e-02, -2.40991186e-02,       -2.09105573e-02, -1.78529031e-02, -1.49314273e-02, -1.21514012e-02,       -9.51809614e-03, -7.03678328e-03, -4.71273396e-03, -2.55121946e-03,       -5.57511040e-04,  1.26310442e-03,  2.87401767e-03,  4.14164198e-03,        4.91367935e-03,  5.03783182e-03,  4.36180141e-03,  2.73329013e-03,        0.00000000e+00];

out_0_4 = [1.00000000e+00 9.79901587e-01 9.59876998e-01 9.39913730e-01 9.19999280e-01 9.00121146e-01 8.80266824e-01 8.60423811e-01 8.40579604e-01 8.20721702e-01 8.00837600e-01 7.80914795e-01 7.60940786e-01 7.40903069e-01 7.20789140e-01 7.00586498e-01 6.80283359e-01 6.59881522e-01 6.39394892e-01 6.18837798e-01 5.98224571e-01 5.77569540e-01 5.56887034e-01 5.36191382e-01 5.15496914e-01 4.94817960e-01 4.74168849e-01 4.53563910e-01 4.33017473e-01 4.12543868e-01 3.92157424e-01 3.71872470e-01 3.51703336e-01 3.31664351e-01 3.11769844e-01 2.92029223e-01 2.72420122e-01 2.52907569e-01 2.33456557e-01 2.14032080e-01 1.94599133e-01 1.75122709e-01 1.55567802e-01 1.35899406e-01 1.16082515e-01 9.60821237e-02 7.58632253e-02 5.53908139e-02 3.46341220e-02 1.50559168e-02 1.80582963e-03 5.89944135e-04 1.40450265e-02 3.57382629e-02 5.85598625e-02 8.07636941e-02 1.02411353e-01 1.23565816e-01 1.44290059e-01 1.64647057e-01 1.84699788e-01 2.04511227e-01 2.24144351e-01 2.43662135e-01 2.63127557e-01 2.82603592e-01 3.02153216e-01 3.21835736e-01 3.41668113e-01 3.61640128e-01 3.81741116e-01 4.01960413e-01 4.22287355e-01 4.42711275e-01 4.63221511e-01 4.83807397e-01 5.04458269e-01 5.25163461e-01 5.45912310e-01 5.66694150e-01 5.87498318e-01 6.08314147e-01 6.29130975e-01 6.49938135e-01 6.70724964e-01 6.91480797e-01 7.12197430e-01 7.32876309e-01 7.53521246e-01 7.74136051e-01 7.94724534e-01 8.15290507e-01 8.35837778e-01 8.56370160e-01 8.76891462e-01 8.97405496e-01 9.17916071e-01 9.38426999e-01 9.58942089e-01 9.79465152e-01 1.00000000e+00];

out_1_4 = [-6.78422293e-18  4.75910166e-03  9.69141151e-03  1.47723234e-02  1.99772310e-02  2.52815281e-02  3.06606086e-02  3.60898662e-02  4.15446947e-02  4.70004879e-02  5.24326395e-02  5.78165433e-02  6.31275932e-02  6.83411829e-02  7.34327061e-02  7.83775567e-02  8.31516137e-02  8.77399083e-02  9.21356306e-02  9.63322562e-02  1.00323261e-01  1.04102121e-01  1.07662311e-01  1.10997307e-01  1.14100586e-01  1.16965622e-01  1.19585892e-01  1.21954871e-01  1.24066035e-01  1.25912860e-01  1.27488821e-01  1.28787395e-01  1.29802056e-01  1.30526281e-01  1.30953545e-01  1.31065066e-01  1.30762974e-01  1.29918001e-01  1.28400804e-01  1.26082036e-01  1.22832351e-01  1.18522405e-01  1.13022850e-01  1.06204343e-01  9.79375361e-02  8.80930851e-02  7.65416441e-02  6.31538673e-02  4.78017956e-02  3.08459856e-02  1.38529197e-02 -1.42857923e-03 -1.38351824e-02 -2.39336467e-02 -3.25118436e-02 -3.98327180e-02 -4.59822976e-02 -5.10464753e-02 -5.51111435e-02 -5.82621951e-02 -6.05855227e-02 -6.21670189e-02 -6.30925765e-02 -6.34480880e-02 -6.33194463e-02 -6.27925439e-02 -6.19532735e-02 -6.08837104e-02 -5.96218687e-02 -5.81774847e-02 -5.65598304e-02 -5.47781776e-02 -5.28417981e-02 -5.07599638e-02 -4.85419466e-02 -4.61970183e-02 -4.37344508e-02 -4.11635159e-02 -3.84934855e-02 -3.57336314e-02 -3.28932255e-02 -2.99815396e-02 -2.70078457e-02 -2.39814155e-02 -2.09115209e-02 -1.78074337e-02 -1.46864418e-02 -1.15972577e-02 -8.59629385e-03 -5.73996252e-03 -3.08467609e-03 -6.86846925e-04  1.39711262e-03  3.11079019e-03  4.39777343e-03  5.20164999e-03  5.46600750e-03  5.13443361e-03  4.15051597e-03  2.45784222e-03  0.00000000e+00];

difference_area_left_a = subtract(polyshape([af_pts_0; af_pts_1]', 'Simplify', false), polyshape([out_0_4; out_1_4]', 'Simplify', false));
difference_area_left_b = subtract(polyshape([out_0_4; out_1_4]', 'Simplify', false), polyshape([af_pts_0; af_pts_1]', 'Simplify', false));
difference_area_left = union(difference_area_left_a, difference_area_left_b);

difference_area_right_a = subtract(polyshape([af_pts_0; af_pts_1]', 'Simplify', false), polyshape([out_0_6; out_1_6]', 'Simplify', false));
difference_area_right_b = subtract(polyshape([out_0_6; out_1_6]', 'Simplify', false), polyshape([af_pts_0; af_pts_1]', 'Simplify', false));
difference_area_right = union(difference_area_right_a, difference_area_right_b);

sprintf('Area difference on the left plot: %f', area(difference_area_left))

sprintf('Area difference on the right plot: %f', area(difference_area_right))