Gm, Gn = 46, 46; 
Gcolptr = [1, 4, 7, 11, 15, 18, 21, 25, 30, 34, 37, 41, 46, 50, 56, 60, 64, 69, 73, 77, 81, 86, 90, 93, 98, 102, 106, 110, 114, 117, 121, 126, 130, 134, 139, 143, 147, 152, 157, 161, 165, 169, 172, 175, 179, 186, 189]; 
Growval = [1, 2, 6, 1, 2, 3, 2, 3, 4, 7, 3, 4, 5, 8, 4, 5, 9, 1, 6, 10, 3, 7, 8, 11, 4, 7, 8, 9, 12, 5, 8, 9, 13, 6, 10, 14, 7, 11, 12, 16, 8, 11, 12, 13, 17, 9, 12, 13, 18, 10, 14, 15, 19, 24, 30, 14, 15, 16, 19, 11, 15, 16, 17, 12, 16, 17, 18, 20, 13, 17, 18, 21, 14, 15, 19, 24, 17, 20, 21, 26, 18, 20, 21, 22, 27, 21, 22, 23, 28, 22, 23, 29, 14, 19, 24, 25, 31, 24, 25, 26, 32, 20, 25, 26, 27, 21, 26, 27, 28, 22, 27, 28, 29, 23, 28, 29, 14, 30, 31, 33, 24, 30, 31, 32, 34, 25, 31, 32, 35, 30, 33, 34, 36, 31, 33, 34, 35, 37, 32, 34, 35, 38, 33, 36, 37, 43, 34, 36, 37, 38, 44, 35, 37, 38, 39, 45, 38, 39, 40, 45, 39, 40, 41, 45, 40, 41, 42, 45, 41, 42, 46, 36, 43, 44, 37, 43, 44, 45, 38, 39, 40, 41, 44, 45, 46, 42, 45, 46]; 
Gnzval = [1.3333333333333333, -0.6666666666666666, -0.6666666666666666, -0.6666666666666666, 1.0666666666666667, -0.4, -0.4, 1.0857142857142859, -0.2857142857142857, -0.4, -0.2857142857142857, 0.8412698412698412, -0.2222222222222222, -0.3333333333333333, -0.2222222222222222, 0.5079365079365079, -0.2857142857142857, -0.6666666666666666, 1.0666666666666667, -0.4, -0.4, 1.5666666666666667, -0.5, -0.6666666666666666, -0.3333333333333333, -0.5, 2.0, -0.5, -0.6666666666666666, -0.2857142857142857, -0.5, 1.4523809523809523, -0.6666666666666666, -0.4, 0.6857142857142857, -0.2857142857142857, -0.6666666666666666, 2.6666666666666665, -1.0, -1.0, -0.6666666666666666, -1.0, 3.6666666666666665, -1.0, -1.0, -0.6666666666666666, -1.0, 2.6666666666666665, -1.0, -0.2857142857142857, 1.4228937728937727, -0.4, -0.3333333333333333, -0.25, -0.15384615384615385, -0.4, 2.4, -1.0, -1.0, -1.0, -1.0, 3.0, -1.0, -1.0, -1.0, 4.0, -1.0, -1.0, -1.0, -1.0, 3.0, -1.0, -0.3333333333333333, -1.0, 2.0, -0.6666666666666666, -1.0, 2.6666666666666665, -1.0, -0.6666666666666666, -1.0, -1.0, 2.9166666666666665, -0.25, -0.6666666666666666, -0.25, 0.7222222222222222, -0.25, -0.2222222222222222, -0.25, 0.9166666666666666, -0.6666666666666666, -0.25, -0.6666666666666666, 2.083333333333333, -0.5, -0.6666666666666666, -0.5, 1.6666666666666665, -0.5, -0.6666666666666666, -0.6666666666666666, -0.5, 1.6666666666666665, -0.5, -0.6666666666666666, -0.5, 1.6666666666666665, -0.5, -0.2222222222222222, -0.5, 1.2222222222222223, -0.5, -0.6666666666666666, -0.5, 1.1666666666666665, -0.15384615384615385, 0.5371794871794872, -0.25, -0.13333333333333333, -0.6666666666666666, -0.25, 2.9166666666666665, -1.0, -1.0, -0.6666666666666666, -1.0, 2.6666666666666665, -1.0, -0.13333333333333333, 0.47320261437908495, -0.2222222222222222, -0.11764705882352941, -1.0, -0.2222222222222222, 3.2222222222222223, -1.0, -1.0, -1.0, -1.0, 3.0, -1.0, -0.11764705882352941, 0.42291021671826623, -0.2, -0.10526315789473684, -1.0, -0.2, 3.2, -1.0, -1.0, -1.0, -1.0, 4.0, -1.0, -1.0, -1.0, 2.6666666666666665, -1.0, -0.6666666666666666, -1.0, 2.5, -1.0, -0.5, -1.0, 2.4, -1.0, -0.4, -1.0, 1.3333333333333333, -0.3333333333333333, -0.10526315789473684, 0.28708133971291866, -0.18181818181818182, -1.0, -0.18181818181818182, 2.1818181818181817, -1.0, -1.0, -0.6666666666666666, -0.5, -0.4, -1.0, 3.7888888888888888, -0.2222222222222222, -0.3333333333333333, -0.2222222222222222, 0.5555555555555556]; 
G = SparseMatrixCSC(Gm, Gn, Gcolptr, Growval, Gnzval); 
