function mesh = mesh_generation_basel()
    mesh.P = [-0.35324073384468713, -0.23148942261705546, -0.14478546688957589, -0.05805347168807790, 0.06366984803808229, 0.03137686415021390, -0.14546033508680190, -0.32094769238976228, -0.30191978697502653, -0.14478542851153825, 0.01234896593136026, -0.00809983115822205, -0.14478542851153825, -0.28140580238990115, -0.34606941133439217, -0.25003713434934871, 0.05649855191268827, -0.03953372267372784, -0.14478542851153825, -0.24131718486252185, -0.16904740745442054, -0.06214653977256362, -0.12277424711824714, -0.22939688242247736, -0.26106175815403176, -0.29730256364571783, -0.31152915587913999, -0.21714677591742637, -0.33875633474485356, -0.44450816637867546, -0.39103265728741010, -0.34581694381246608, -0.39322094420958287, -0.44101019774020039, -0.50000000000000000, -0.48122210591278192, -0.43552248136272542, -0.40088552764440810, -0.43676728980366780, -0.39915990475626456, -0.37749023752671201, -0.34691260560466625, -0.37611246117851616, -0.32733115936141921, -0.29169487272535471, -0.35306384707062000, -0.29115520315881949, -0.31864371247452450, -0.35520744516659919, -0.35799974243539495, -0.37518036418407058;
              -0.49810745903306392, -0.40315542066634780, -0.50000000000000000, -0.40315542066634780, -0.49810745903306392, -0.28345542466485962, -0.25454102733353401, -0.28345542466485962, -0.16071542035068842, -0.09541987623994347, -0.16071542035068842, -0.03095560282290616, 0.05642238733893190, -0.03095564120094380, 0.06389719881341827, 0.14229482819243922, 0.06389719881341827, 0.14229482819243922, 0.21389584071061163, 0.25013958691943117, 0.32776999303270715, 0.30867847841922591, 0.41964564790294845, 0.39309315493878216, 0.30988431636158165, 0.34857136915400244, 0.42594458724672390, 0.48496307708965236, 0.50000000000000000, 0.44536233916633938, 0.40178484499746236, 0.35715005987813386, 0.32326254048553160, 0.32782643273429779, 0.34197025073643261, 0.23101655913980046, 0.25525950977189982, 0.28029188704355301, 0.17960062422837664, 0.20914938893409107, 0.24768536418275133, 0.21481168700474795, 0.16884856124694658, 0.18238314948934065, 0.22752468541342663, 0.25844000808721201, 0.27989100685285428, 0.31228299168604579, 0.28433473986442626, 0.31223451542726355, 0.28769483540271312] + [.2; 0];

    %         1  2  3  4  5  6  7  8  9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51
    mesh.B = [2, 2, 2, 2, 0, 0, 1, 2, 2, 1, 0, 0, 1, 2, 0, 0, 0, 2, 1, 0, 1, 2, 2, 1, 0, 0, 1, 2, 2, 2, 1, 0, 0, 1, 2, 2, 1, 0, 2, 1, 0, 1, 2, 2, 2, 0, 2, 2, 0, 2, 0];

    mesh.F = [1, 2, 4, 2, 2, 4, 8,  7,  7,  6,  9, 10, 10, 11, 14, 12, 13, 13, 13, 13, 16, 18, 19, 19, 20, 21, 25, 21, 23, 24, 24, 24, 27, 26, 27, 27, 29, 32, 31, 30, 30, 33, 34, 34, 35, 37, 37, 36, 37, 39, 40, 40, 42, 41, 42, 42, 45, 46, 47, 48, 49;
              2, 3, 5, 7, 4, 6, 7, 10, 11, 11, 10, 13, 12, 12, 13, 17, 16, 19, 17, 18, 19, 22, 21, 22, 21, 22, 21, 23, 28, 26, 27, 28, 28, 27, 31, 29, 30, 31, 34, 34, 35, 34, 37, 35, 36, 41, 40, 39, 39, 43, 42, 43, 43, 42, 44, 45, 47, 47, 48, 50, 50;
              8, 4, 6, 8, 7, 7, 9,  9, 10,  7, 14, 14, 13, 10, 15, 13, 15, 16, 18, 19, 20, 19, 20, 21, 25, 23, 24, 24, 24, 25, 26, 27, 29, 32, 32, 31, 31, 33, 33, 31, 34, 38, 38, 37, 37, 38, 41, 37, 40, 40, 41, 42, 44, 46, 45, 46, 46, 49, 49, 49, 51];
end