function y = f(x)

% y = cos(x);
% for step function initial data:
if x < pi 
  y = 1;
else
  y = 0;
end