function output = build_parabola_se(s)
f = inline('x.^2','x');
xsize = 2*s+1;
ysize = f(s)+1;
offset = s;
temp = zeros(xsize,ysize);
for i = -s:1:s
    temp(int16(f(i))+1:ysize,int16(i+offset)+1) = 1;
end
output = temp;
end