function ai = gweights(n) % this generates the appropriate % gaussian weihgts on the interval [-1,1] switch n case 1 % this means n=1 ai = [2]; case 2 % this means n=2 ai = [1,1]; case 3 % this means n=3 ai = [5/9, 8/9, 5/9]; case 4 % this means n=4 ai = [0.347854845137454, 0.652145154862546,... 0.652145154862546, 0.347854845137454]; case 5 % this means n=5 ai = [0.236926885056189, 0.478628670499366,... 0.568888888888889, 0.478628670499366, ... 0.236926885056189]; case 6 % this means n=6 ai = [0.17132449237917, 0.360761573048139,... 0.467913934572691, 0.467913934572691, ... 0.360761573048139, 0.17132449237917]; end