function y=prvaf(x);
       y=[];
       for xi=x,
       if xi<0, yi=pi*exp(xi);
       elseif xi<1/2, yi=pi+sin(pi*xi);
       elseif xi<1/2+sqrt(pi+1), yi=pi+1-(xi-1/2)^2;
       else yi=0; end;
       y=[y,yi];
       end;
