function y=Sierpfun(x,a)
%We compute y=1/2(x-a)+a
% x and a have to be two dimensional vectors for this to make sense
y=1/2*(x-a)+a;