function w = mysimpweights(n)
if rem(n,2) ~= 0
    error('n must be even')
end
w = ones(n+1,1);
for i = 2:n
    if rem(i,2)==0
        w(i)=4;
    else
        w(i)=2;
    end
end