p = [ 0.000 0.000
1.000 0.000
1.000 1.000
0.000 1.000
0.500 0.000
1.000 0.500
0.500 1.000
0.000 0.500
0.500 0.500
0.250 0.000
1.000 0.250
0.750 1.000
0.000 0.750
0.750 0.000
1.000 0.750
0.250 1.000
0.000 0.250
0.750 0.750
0.250 0.250
0.750 0.250
0.250 0.750
0.500 0.250
0.750 0.500
0.500 0.750
0.250 0.500
0.125 0.000
1.000 0.125
0.875 1.000
0.000 0.875
0.625 0.000
1.000 0.625
0.375 1.000
0.000 0.375
0.375 0.000
1.000 0.375
0.625 1.000
0.000 0.625
0.875 0.000
1.000 0.875
0.125 1.000
0.000 0.125
0.875 0.875
0.125 0.125
0.875 0.625
0.625 0.875
0.750 0.875
0.875 0.750
0.375 0.125
0.125 0.375
0.250 0.125
0.125 0.250
0.875 0.125
0.625 0.375
0.875 0.250
0.750 0.125
0.125 0.875
0.375 0.625
0.125 0.750
0.250 0.875
0.500 0.125
0.500 0.375
0.625 0.125
0.375 0.250
0.625 0.250
0.875 0.500
0.625 0.500
0.875 0.375
0.750 0.625
0.750 0.375
0.500 0.875
0.500 0.625
0.375 0.875
0.625 0.750
0.375 0.750
0.625 0.625
0.125 0.500
0.375 0.500
0.125 0.625
0.250 0.375
0.250 0.625
0.375 0.375]';
e = [ 1 26 0.000 0.125 1.000 1 0
2 27 0.000 0.125 2.000 1 0
3 28 0.000 0.125 3.000 1 0
4 29 0.000 0.125 4.000 1 0
5 30 0.500 0.625 1.000 1 0
6 31 0.500 0.625 2.000 1 0
7 32 0.500 0.625 3.000 1 0
8 33 0.500 0.625 4.000 1 0
10 34 0.250 0.375 1.000 1 0
11 35 0.250 0.375 2.000 1 0
12 36 0.250 0.375 3.000 1 0
13 37 0.250 0.375 4.000 1 0
14 38 0.750 0.875 1.000 1 0
15 39 0.750 0.875 2.000 1 0
16 40 0.750 0.875 3.000 1 0
17 41 0.750 0.875 4.000 1 0
26 10 0.125 0.250 1.000 1 0
27 11 0.125 0.250 2.000 1 0
28 12 0.125 0.250 3.000 1 0
29 13 0.125 0.250 4.000 1 0
30 14 0.625 0.750 1.000 1 0
31 15 0.625 0.750 2.000 1 0
32 16 0.625 0.750 3.000 1 0
33 17 0.625 0.750 4.000 1 0
34 5 0.375 0.500 1.000 1 0
35 6 0.375 0.500 2.000 1 0
36 7 0.375 0.500 3.000 1 0
37 8 0.375 0.500 4.000 1 0
38 2 0.875 1.000 1.000 1 0
39 3 0.875 1.000 2.000 1 0
40 4 0.875 1.000 3.000 1 0
41 1 0.875 1.000 4.000 1 0]';
t = [ 2 52 38 1
4 56 40 1
9 57 77 1
9 53 66 1
5 48 34 1
7 45 36 1
8 49 76 1
6 44 65 1
20 53 64 1
21 57 74 1
21 56 58 1
20 52 54 1
19 49 51 1
18 44 47 1
19 48 63 1
18 45 73 1
14 62 30 1
16 72 32 1
25 78 76 1
23 67 65 1
10 43 26 1
12 42 28 1
25 81 77 1
23 75 66 1
22 62 64 1
24 72 74 1
13 78 58 1
11 67 54 1
17 43 51 1
15 42 47 1
22 81 63 1
24 75 73 1
20 55 52 1
21 59 56 1
21 80 57 1
20 69 53 1
19 50 48 1
18 46 45 1
19 79 49 1
18 68 44 1
9 61 53 1
9 71 57 1
4 29 56 1
2 27 52 1
8 33 49 1
6 31 44 1
5 60 48 1
7 70 45 1
22 60 62 1
24 70 72 1
13 37 78 1
11 35 67 1
17 41 43 1
15 39 42 1
22 61 81 1
24 71 75 1
14 55 62 1
16 59 72 1
25 80 78 1
23 69 67 1
10 50 43 1
12 46 42 1
25 79 81 1
23 68 75 1
14 38 55 1
16 40 59 1
25 77 80 1
23 66 69 1
10 34 50 1
12 36 46 1
25 76 79 1
23 65 68 1
22 64 61 1
24 74 71 1
13 58 29 1
11 54 27 1
17 51 33 1
15 47 31 1
22 63 60 1
24 73 70 1
5 30 60 1
7 32 70 1
8 76 37 1
6 65 35 1
1 26 41 1
3 28 39 1
9 77 61 1
9 66 71 1
20 64 55 1
21 74 59 1
21 58 80 1
20 54 69 1
19 51 50 1
18 47 46 1
19 63 79 1
18 73 68 1
55 38 52 1
59 40 56 1
80 77 57 1
69 66 53 1
50 34 48 1
46 36 45 1
79 76 49 1
68 65 44 1
61 64 53 1
71 74 57 1
29 58 56 1
27 54 52 1
33 51 49 1
31 47 44 1
60 63 48 1
70 73 45 1
60 30 62 1
70 32 72 1
37 76 78 1
35 65 67 1
41 26 43 1
39 28 42 1
61 77 81 1
71 66 75 1
55 64 62 1
59 74 72 1
80 58 78 1
69 54 67 1
50 51 43 1
46 47 42 1
79 63 81 1
68 73 75 1]';

pdemesh(p,e,t)
axis equal
